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ABSTRACT 

This  thesis  develops  project  planning  models  that  allow  the 
possibility  of  specifying  alternate  ways  of  performing  any  of  the  jobs 
in  the  project.  The  "job  alternatives"  for  any  task  may  have  different 
times,  costs,  resource  requirements  and  possibly  different  precedence 
relations  with  other  jobs  in  the  project  network.  The  problem  is  to 
select  the  particular  way  in  which  each  job  will  be  performed  and  schedule 
the  resulting  jobs  so  as  to  minimize  the  cost  of  the  jobs  plus  the  cost 
associated  with  the  completion  date  of  the  project. 

The  problem  of  selecting  the  optimal  job  alternatives  in  networks 
wltii  no  resource  constraints  is  formulated  as  an  integer  programming 
probJem.  One  constraint  in  required  for  eacii  set  of  job  alternatives 
and  one  for  each  possible  path  in  the  original  project  network.  Arguments 
are  developed  to  show  that  a  substantial  number  of  the  precedence  con¬ 
straints  are  x'edundant  and  may  be  climinaccd.  To  accomplish  this  re¬ 
duction  in  problem  size  an  algorithm  related  to  the  critical  path 
algorithm  is  developed  to  reduce  each  network  to  an  equivalent  network 
containing  only  job  alternatives  and  maximal  distances  between  them. 

Jobs  with  no  alternative  are  eliminated. 

Two  branch  and  bound  routines  are  then  developed  to  solve  the 
problem.  One  of  these  is  tested  on  a  series  of  problems  and  is  sbo^’n 
to  he  efficient.  An  integer  programming  algorithm  is  developed  to 
serve  as  a  sub-routine  in  the  branch  and  bound  algorithms.  It  is  fast 
in  that  it  uses  the  critical  path  algorithm  to  solve  problems. 

When  resource  requirements  are  added  to  the  tasks  of  the  project, 
and  the  total  availability  of  resource  per  period  is  constrained,  the 
problem  of  scheduling  the  joba  so  as  to  minimize  tomplction  date  becomes 
extremely  difficult.  Nine  heuristic  routln-s  for  the  loading  problem 


are  developed  and  tested.  Of  these  a  serial  loading  rule,  operating 
on  a  job  list  ordered  by  late  start,  with  no  job  bumping,  proves 
superior.  Three  methods  of  generating  combinations  of  job  alternatives 
to  be  loaded  were  examir  o-f.  These  were  complete  enumeration,  pairwise 
interchange  and  multiple  pairs  interchange.  None  of  the  methods  pro¬ 
vided  good  solutions  in  reasonable  amounts  of  time. 

To  show  the  generality  of  the  planning  modal  developed  the  Integer 
programming  formulation  of  the  project  problem  was  adapted  to  the  m  x  n 
job-shop  scheduling  problem,  the  single  product  assembly-line  balancing 
problem  and  the  problem  of  planning  projects  under  incentive  contracts. 
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Chapter  I 


DECISION  NETWORK  PLANNING  MODELS 

The  growth  of  Interest  In  quantitative  solutions  to  management 
problems  has  resulted  in  a  rapid  development  of  planning  models,  based 
on  network  representation  of  the  activities  to  be  performed.  Process 
charts  have  been  used  to  show  basic  work  elements  in  a  single  task  and 
the  order  in  which  they  must  be  performed.  Networks  have  been  used  to 
show  the  required  Job  ordering  in  large  construction  projects  and  net¬ 
work  based  algorithms  have  been  de 'eloped  to  find  the  total  time 
required  to  complete  surh  a  project. 

Although  the  applications  of  the  models  are  at  different 
levels  of  detail,  they  have  many  common  characteristics.  In  each  case 
there  may  be  constraints  In  the  problem  that  effect  the  "time"  at 
which  the  individual  planning  units,  either  work  elements  or  Jobs  may 
be  performed.  These  may  take  the  form  of  an  explicit  restriction  that 
a  Job  must  start  on  a  particular  day  or  that  the  Job  cannot  be  started 
before  a  given  day.  Alternately,  e  Job  may  be  constrained  not  to 
start  until  some  prior  Job  is  finished.  For  example,  the  cellar  walls 
of  a  house  cannot  be  constructed  until  the  footings  are  laid.  This 
second  type  of  time  constraints  will  be  called  "precedence"  constraints. 
The  graph  of  Figure  i-l  shows  a  series  of  tasks  to  related  by 
precedence  constraints  which  are  graphically  illustrated  by  directed 
line  segrant.  In  a  particular  case,  say  Q — we  Imply 
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that  S]  Is  a  predecessor  of  $2  and  conversely  $2  Is  a  successor 
of  S|.  The  nodes  representing  planning  units  and  the  directed  line 
segments  make  up  the  network  to  which  we  heve  referred. 


In  addition  to  the  precedence  relations,  planning  units  may 
be  related  by  a  mutual  dependence  on  a  limited  resource.  In  a  process 
chart  of  the  man-machine  variety,  both  the  man  and  the  machine  are 
considered  to  be  resources  and  they  may  be  physically  restricted  to 
perfonr  only  one  task  at  a  time.  Thus,  if  the  man  Is  required  to 
perform  both  and  S]^,  even  though  no  technological  constraint 
exists  between  them,  the  Jobs  must  be  performed  serially.  That  Is,  he 
must  perform  then  of  then  S^,  but  not  both  together. 
For  problems  of  practical  Interest,  the  number  of  feasible  sequences 
may  be  large  and  the  problem  of  find  the  "best"  sequence  on  all 
resources  Is  a  difficult  combinatorial  problem. 

Finally,  the  planning  units  may  be  related  technically  by  the 
nature  of  the  project  that  Is  being  performad.  It  Is  conceivable  that 
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In  a  construction  project  there  may  be  two  methods  for  performing  a 
particular  job^wlth  different  costs  and  different  performance  times. 

For  example,  If  wooden  partitions  are  required,  they  may  either  be 
purchased  In  an  assembled  form  and  be  quickly  Installed  or  they  could 
be  fabriceted  on  the  site  by  carpenters.  We  call  this  kind  of  mutually 
exclusive  alternative  a  "job  alternative"  Interdependency  because  we 
must  choose  between  the  two  methods  of  performing  the  task.  Many 
additional  types  of  Interdependency  could  exist  between  planning  units. 
Perhaps  if  the  partitions  are  pre>ass«mbled,  then  a  particular  design 
for  the  electrical  system  Is  required.  We  wi 1 1  term  all  such  relations 
which  are  not  the  "Job  alternative"  interdependency  described  above, 
"other"  interdependency. 

This  thesis  will  develop  network  models  that  include 
precedence  constraints  and  the  possibility  of  resource  constraints. 

Each  model  that  is  discussed  will  include  sets  of  mutually  exclusive 
tasks,  that  is  "Job  alternative"  interdependencies  or  as  we  shall  term 
them,  "decision"  nodes.  In  addition,  '^ther"  interdependencies  may  be 
Imposed  on  the  sets  of  mutually  exclusive  tasks.  The  functional  setting 
the  models  will  be  the  (construction)  project  scheduling  problem 
but.  In  fact,  the  theory  developed  would  be  applicable  to  a  wide  range 
of  planning  models. 

in  one  chapter  of  tlie  thesis,  we  will  examine  the  problem  of 
selecting  from  the  sets  of  "Job  alternatives"  the  particular  Jobs  we 
wish  to  perform.  In  terms  of  our  original  example,  this  might  be  the 
choice  of  pre-assembled  partitions.  If  there  is  a  large  number  of 
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decision  nodes  in  the  network,  there  will  be  very  many  possible  combine* 
tions  of  decision  Jobs  that  we  can  select.  To  evaluate  a  particular 
solution,  that  is,  the  choice  of  a  particular  set  of  decision  jobs, 
one  from  each  mutual Jiy  exclusive  set,  we  must  evaluate  the  cost  of  the 
Jobs  plus  the  effect  that  the  choice  of  these  Jobs  have  on  the  comple¬ 
tion  date  and  thus  the  completion  cost  of  the  project.  We  will  call 
the  choice  of  a  set  of  decision  Jobs  a  "design"  problem  and  the  calcu¬ 
lation  of  the  minimum  time  for  completion  of  the  project,  given  a  choice 
of  decision  Jobs,  an  "operating"  problem.  Note  that  it  is  necessary 
to  solve  an  "operating"  problem  to  properly  evaluate  any  "design"  and 
that  an  optimal  "design"  Is  one  which  minimizes  the  sum  of  Job  cost  and 
completion  date  cost. 

The  interaction  of  "design"  and  "operating"  prob|)'ms  can  be 
seen  In  many  areas  of  planning,  if  we  wish  to  establish  warehouses  In 
a  manufacturer's  distribution  system,  the  "design"  decision  is  the 
selection  of  the  quantity,  size  and  location  of  the  warehouses.  To 
evaluate  such  a  "design"  we  must  find  the  total  cost  of  establishing 
the  warehouses  plus  the  minimum  cost  for  'V}peratlng"  the  warehouses. 

The  'Operating"  problem  is  the  optimal  allocation  of  customer  demands 
to  warehouses  and  vtarehouse  demands  to  factories  so  as  to  minimize 
production,  shipping  and  Inventory  costs. 

The  problem  may  be  Illustrated  graphically  with  the  following 
design  problem.  If  we  have  two  design  variables,  x  and  y,  each 
having  feasible  levels  I,  2,  3.  ell  possible  designs  are  represented 
by  peths  In  the  tree  of  Figure  1-2.  For  each  path  or  each  possible 
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design,  It  may  be  required  that  we  solve  an  operating  problem.  In  a 
facilities  problem,  for  example,  plant  layout,  this  Implies  that  we 
determine  the  best  method  of  scheduling  production  for  each  possible 


layout.  These  illustrations  suggest  that  if  problems  have  many 
discrete  design  variables  or  If  the  operating  problem  vm  must  solve  is 
a  complex  one,  the  determination  of  an  overall  optimum  solution  may  be 
difficult.  If  methods  can  be  found  to  reduce  the  number  of  designs 
that  it  Is  necessary  to  evaluate  or  if  efficient  methods  can  be  found 
for  solving  the  operating  problems,  then  it  will  not  be  necessary  to 
solve  design  problems  sequentially,  in  this  thesis,  methods  for  the 
elimination  of  some  designs  and  methods  for  solving  operating  problems 
In  the  area  of  project  planning  will  be  developed. 

Chapter  II  contains  a  review  of  a  wide  variety  of  planning 
literature.  The  work  is  categorized  by  the  particular  sets  of  con¬ 
straints  that  arc  to  be  found  in  the  models.  That  Is,  we  consider  most 
combinations  of  precedence,  resource  and  both  "Job  alternative"  and 
'Mother"  types  of  Interdependencies. 


Chapter  III  develops  a  model  for  project  planning  that 
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contains  precedence  constraints  and  both  "Job  alternative"  and  "other" 
interdependency  constraints.  It  is  assumed  that  in  the  project  there 
are  a  number  of  competing  methods  for  performing  some  of  the  jobs,  each 
method  having  a  different  cost,  a  different  time,  possibly  different 
precedence  relations  with  other  jobs  and  different  interdependencies 
with  other  jobs.  All  possible  jobs  are  considered  in  the  project  graph 
and  then  in  the  scheduling  phase  the  job  alternatives  that  minimize 
total  cost  are  selected.  A  numerical  problem  is  introduced  here  that 
will  be  used  to  illustrate  the  material  of  Chapters  III,  IV  and  VI.  In 
Chapter  IV  methods  are  developed  to  reduce  the  original  decision  network 
so  that  the  problem  may  reasonably  be  solved  with  standard  integer 
programming  techniques. 

In  the  solution  of  a  decision  network,  it  may  be  necessary 
to  solve  sub-problems  that  minimize  the  wost  of  the  "design"  selected 
with  no  regard  to  the  cost  of  the  "operating"  problem,  that  is,  the 
cost  of  the  minimun  completion  date.  Essentially,  the  sub-problem  is 
to  find  the  set  of  decision  jobs  that  meets  the  'job  alternative"  and 
"other"  interdependencies  with  minimum  job  cost.  Chapter  V  develops 
an  integer  programming  algorithm  specifically  for  this  problem.  Chapter 
VI develops  two  branch  and  bound  algorithms  which  solve  for  the  best 
"design"  given  the  cost  of  the  decision  jobs,  the  cost  of  the  completion 
date  and  the  interdependency  constraints.  They  solve  the  "design"  and 
the  operating  "problem"  simultaneously.  Computational  results  are 
given  for  one  of  these  methods. 

In  Chapter  VII  we  consider  the  full  model,  that  Is,  project 
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planninq  problems  with  precedence  resource  and  interdependency 
constraints.  To  solve  the  "operating"  problem  In  models  having  no 
limit  on  resource  usage  Is  relatively  straight  forward.  It  Is  simply 
a  matter  of  calculating  the  Itjngch  of  the  critical  path  and  evaluating 
the  cost  of  that  finish  date.  When  we  add  resource  constraints,  the 
problem  of  calculating  minimum  project  length  for  any  given  design  may 
be  a  very  complex  combinational  problem.  In  fact,  for  large  projects, 
given  current  techniques  and  reasonable  limits  computer  time,  it  Is  not 
possible  to  find  optimum  or  minimum  length  projects.  For  this  reason, 
we  develop  end  experimentally  test  several  heuristic  loading  techniques. 
The  best  of  these  heuristics  is  then  used  as  our  "operating"  rule  to 
evaluate  various  designs.  The  designs  to  be  tested  are  generated  first 
by  complete  enumeration,  then  by  pairwise  Interchange  and  finally  by 
multiple  pairs  interchange. 

As  we  have  stated  above,  many  nlannlnn  nrnhlenxi  can  be 
represented  by  the  combination  of  constraints  we  have  discussed.  To 
illustrate  this  point,  in  Chapter  V(il  I  the  basic  integer  programming 
formulation  of  our  decision  network  planning  problem  is  used  to  formu¬ 
late  the  Job-shop  scheduling  problem  and  the  assembly-line  balancing 
problem.  These  formulations  prove  to  be  substantially  more  compact  than 
competitive  formulations  of  the  problem.;.  Finally,  the  model  is  adapted 
to  projects  with  more  complex  criterion  functions,  specifically  the  cost 
structure  of  incentive  contracts.  Chapter  ix  contains  a  summary  of  the 
work,  conclusions  and  recommendations  for  further  research. 

Several  terms  that  are  common  in  the  literature  of  project 
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schedullngt  such  as  aarly  start,  will  be  used  frequently  In  the 
following  chapters.  These  terms  are  defined  rigorously  elsewhere 
&5,  4^  so  that  we  will  review  them  only  briefly  here. 

A  "path"  through  a  project  network  Is  connected  sequence  of 
nodes  (jobs)  and  directed  line  segments  extending  from  one  node  to  some 
other.  In  Figure  l-l  we  have  a  path  from  S|  through  $2  and  to 
S^.  The  nodes  S|,  Sj|^  and  do  not  lie  on  a  common  path.  The 
length  of  any  path  Is  simply  the  sun  of  the  job  times  for  all  Jobs  on 
the  path.  Irfe  now  define  the  "early  finish"  of  a  job  to  be  the  longest 
path  from  the  first  Job  in  the  network  to  the  Job  under  consideration. 
"Early  start"  time  is  simply  early  finish  time  of  a  Job  less  its  Job 
time.  The  "critical  path"  of  a  project  Is  the  longest  path  from  the 
first  Job  in  the  network  to  the  final  Job  in  the  network  and  is 
equivalent  to  the  minimum  numbe**  of  days  required  tc  complete  the 
project. 

The  "latest  start"  date  for  a  Job  Is  defined  as  the  day  on 
which  the  Job  must  start  if  the  project  is  to  finish  exactly  on  its 
due  date.  We  may  calculate  the  value  of  late  start  for  a  Job  by  sub¬ 
tracting  the  length  of  the  longest  path  from  the  Job  to  tl«e  end  of  the 
project  from  the  due  date  of  the  project.  Thus  a  Job  can  begin  no 
earlier  than  the  early  start  time  because  Its  predecessors  must  first 
be  con^leted  and  no  later  than  late  start  or  ! t  wl 1 1  delay  the  finish 
of  the  project  beyond  the  due  date.  The  difference  between  late  start 
and  early  start  time  is  defined  as  Job  "slack"  time,  the  measure  of 
permissible  delay  for  a  Job.  Other  terms  more  uniquely  related  to  the 
models  to  be  discussed  will  be  defined  as  required. 


Chapter  II 

A  REVIEW  OF  SELECTED  PLANNING  LITERATURE 


The  management  planning  literature,  like  many  other  direas  of 
management  activ'ty,  is  susceptible  to  many  possible  categorizations. 
Perhaps  the  most  obvious  breakdown  follows  the  functional  area  of  an 
industrial  concern  and  within  this  is  a  sub-grouping  by  problem  area. 

Thus  under  the  production  heading  we  have  extensive  and  largely  indepen¬ 
dent  literature  growing  up  around  the  "Job-shop  problem"  or  the 
"assembly-line  balancing  problem".  In  finance  we  have  the  "capital 
budgeting  of  interrelated  projects". 

A  second  categorization  might  be  by  solution  technique.  Here 
is  a  list  of  functional  irea  problems  best  solved  by  linear  programming; 
this  list  requires  integer  linear  programming  and  so  on.  This  approach 
is  closely  related  to  a  third  categorization,  the  one  to  be  used  as  the 
framework  for  our  discussion.  The  third  structure  divides  planning 
problems  by  the  type  of  constraint  found  in  the  problem.  To  be  explicit, 
we  have  defined  in  Chapter  I  time  constraints,  resource  constraints  and 
Interdependency  constraints  as  possible  dimensions  of  planning  problems. 
For  simplicity  we  wi 1 1  not  discuss  the  dimension  "uncertainty",  nor  will 
we  consider  motivational  and  social  problems  of  planning. 

The  possible  combinations  of  the  three  dimensions  and  therefore 
cur  subheadings  will  be  simple  time  constraints,  simple  resource  limits, 
simple  interdependency,  problems  with  time  and  interdependency,  with 
time  sr,d  "ttrurtc  cc*?  ,  v?' th  precedersce  end  interdependency  and 

finally  models  with  all  three  characteristics.  To  be  included  in  any 
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two  or  throo  dimension  category  the  model  must  emphasize  some  Interesting 
Interaction  of  the  relevant  dimensions. 

The  following  symbols  will  be  used  throughout  this  literature 
review  and  the  rest  of  the  thesis: 


Si 

d| 

''I 

C| 

t| 

U| 

Li 

•I 

t? 

Wi 

Wf 


•n  individual  Job  or  planning  unit 

a 

I  If  task  S|  Is  to  be  performed 

0  otherwise 

I 

I  if  task  S{  is  on  the  critical  path 

0  otherwise 

the  cost  or  revenue  of  task  S| 
the  time  required  to  perform  task  S) 
the  maximum  length  of  task  S|  If  tf  is  variable 

the  minimum  length  of  task  S]  if  t{  is  variable 

the  reduction  In  cost  C]  per  unit  Increase  In  t| , 
the  time  to  perform  Si 
the  early  start  time  of  task  S| 

the  early  start  time  of  Sf,  an  artificial  FINISH  Job 
that  Is  constrained  to  start  after  all  other  Jobs  In  a 
project  are  finished 

the  desired  con^letion  date  or  due  date  of  a  project 
Wf  -  0  ,  Wf  -  0>0 

0  otherwise 


{: 

{ 


f  ' 


df  >0 


•»! 


^0  otherwise 
the  usage  per  time  period  of  resource  r  by  Job  S| 
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K**  the  avatlabllity  per  time  period  of  resource  r 
P|  the  probability  of  job  S}  occur! ng 

Sj  an  artificial  START  Job  that  must  be  completed  before 

any  other  task  can  be  started 


Simple  Time  Constraints 

The  basic  PERT  or  CPh  model  |k3,  44,  4^  Is  taken  as  the 
example  of  simple  time  constraints.  If  lob  S|  Is  an  Immediate  prede¬ 
cessor  of  Job  Sni  this  relation  will  Pe  shown  symbolically  as  S;  S,n, 
gr^iphically  as 


and  In  mathematical  programming  notation  as  W|  +  t|  Kn  >  An 
alternate  graphical  notation  in  common  use  represents  the  Jobs  as  lines 
rather  than  as  nodes  <nd  represents  immediate  predecessor  relations  by 
intersections  of  Jobs  rather  than  by  directed  line  segments. 


The  length  of  the  critical  path  In  any  network,  say  that  of 
Figure  I  -  I,  may  be  determined  as  the  value  of  Wf  in  a  solution  to 
the  following  problem. 

HInIfflIze  Wf 

Subject  to  Wi  t|  4-  '^2 
•^1  4-  tj  4  Wj 
«|  -h  «l  4  «4 

*2  4 

*')  +  *)  4  ^6 
Wi4  ^4  4 

“5  -h  *5  4 


12 

*^6  4*  ^6  4*  ^7 

'^7  +  ^7  4 

Here  v#e  minimize  the  length  of  the  project,  subject  to  a  set 
of  time  constraints,  one  for  each  Immediate  predecessor  relation  In  the 
graph. 

The  dual  problem  formulated  by  Charnes  and  Cooper  [l^  defines 
a  variable  V|  for  each  Job  in  the  network.  Then  all  Jobs  that  have  no 
predecessors  are  included  in  an  equation  of  the  form 

VI  =r  I 

to  initiate  an  artificial  flow  of  one  unit  Into  tne  network.  For  all 
other  jobs,  not  including  Sf  ,  they  constrain  the  flow  in  from  the 
predecessors  uf  the  JoL  -'xirnurn  of  I  unit)  to  equal  the  flow  out  to 

Immediate  successois. 

-  V3  -  Vi.  -f  Vfi  =  0 

ror  the  final  node  they  establish  a  sink  for  the  one  unit  flow. 

-  V5  -  Vg  =  -  I 

Those  constraints  guarantee  that  e  set  of  Jobs  will  be  chosen  that  will 
form  a  path  through  the  network.  Finally  the  criterion  function  is 
N 

Maximize  ^  t;  V|  ,  Ns  t>tai  tiumber  of  Jobs  in  the  project 

i^ 

Thus  we  select  the  longest  path  in  the  netw..rk. 

A  third  possible  formulation  would  establish  a  cciistraint  for 
each  path  In  the  network  and  constrain  to  be  longer  than  all 

paths  .  Since  this  approach  Implies  complete  knowledge  of  ail  paths, 
it  would  be  a  simple  matter  to  select  the  longest  path  directly. 
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Simple  Interdependency. 

"Job  al terpatl ve"  Interdependency  has  been  defined  as  a  set  of 
mutually  exclusive  alternative  planning  units.  Originally  we  considered 
the  units  to  be  alternative  methods  of  performing  some  Job,  but  we  may 
also  consider  them  different  results  of  some  stochastic  process.  A 
natural  source  of  such  outcomes  would  be  a  research  or  development 
project. 

Eisner  [loj  proposes  a  network  with  decision  nodes  to 
represent  such  situations. 


occur  wi  th  probabi 1 i ty  Pj j  or  Sj^2  will  occur  wi th  probabi 1 i ty 
P;^2  where  Pj  j  -f-  Pj ^2  —  * •  Essentially  he  constructs  a  decision 
tree  with  time  value  and  Job  name  labelling.  Then  using  standard 
probability  calculations  he  obtains  the  probability  of  various  network 
outcomes  and  attaches  to  them  the  sum  of  the  times  from  the  relevant 
path.  We  categorize  this  as  simple  interdependency  because  there  is  no 
Interaction  t<etween  times  and  probabilities. 

Other  work  that  can  be  classed  as  pure ! nterdependency 
emphasi'.:es  solution  techniques  for  0,1  combinatorial  problems  rather 
ihan  the  application  of  techniques  to  planning  problems.  Examples  of 
such  articles  are  references  [62,  68]. 

Simple  Resource  Constraints. 


The  general  knapsack  problsm  may  be  interpreted  as  an  example 
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of  •  simple  resource  constraint.  We  have  a  number  of  available 
"planning  units"  which  may  be  selected,  each  using  some  amount  of  a 
scarce  resource,  or  resources.  The  object  is  to  select  the  set  of 
these  units  which  will  optimize  a  linear  criterion  function,  subject 
to  constraints  in  the  amount  of  each  available  resource.  In  our  nota¬ 
tion  this  may  be  written 

Max  ^  C|  d| 

m 

Subject  to  Xj  ^  K*”  r  s!  1 ,  2^..R 

i'sl 

0  ^  ^  I  integer 

Here  again  d|  =;  I  implies  that  a  planning  unit  is  selected  (for  the 
knapsack)  and  d|=  0  implies  that  is  is  rejected.  C|  is  the  value 
to  us  of  uni t  S{ . 

In  his  discussion  of  this  problem,  Dantzig  [24]  points  out 
that  linear  programming  solutions  of  this  problem  give  values  of  the 
d{  which  will  not  all  be  0-1  but  instead  have  fractional  values.  If 
an  Integer  solution  is  required,  he  suggests  that  rounding  is  usually 
"good  enough"  for  most  practical  problems.  For  exact  integer  solutions 
in  problems  of  one  constraint,  he  sug^sts  the  dynamic  programming 
approach  of  Bellman  [^s].  If  problems  have  two  or  more  constraints,  he 
suggests  the  use  of  linear  programming  with  constraints  added  to  elim¬ 
inate  fractional  extreme  points.  More  recently  Glover  [4^  and 
Weingartnar  and  Ness[743  have  developed  truncated  enumeration  methods 
for  the  solution  of  this  problem.  More  gene. ally  the  large  number  of 
integer  programming  routines  now  available  may  be  applied  to  this  problem. 
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Weingartner  [75i  7^  shows  that  the  Lorie-Savags  [52^  capital 
budgeting  problem  Is  essentially  the  knapsack  problem  as  described 
above.  The  cutting  stock  problem  has  been  formulated  by  Gilmore  and 
Gomory  [37)  and  Pierce  C6I1  as  a  form  of  the  knap..ack  problem.  To 
begin  we  wM I  define  a  resource  as  one  of  the  possible  customer  order 
widths.  Thus  given  customer  orders  for  K**  units  of  item  r  we  must 
schedule  our  cuts  so  as  to  produce  at  least  K*'  units.  The  planning  unit, 
S|,  is  one  pattern  of  cuts  that  may  be  made  from  a  stock  roll.  Thus  we 
must  enumerate  all  possible  combinatioits  of  order  sizes  that  might  be 
taken  from  each  stock  roll  and  define  each  different  combination  as  a 
separate  unit  S;  —  for,  say,  a  total  of  m  units.  Then,  if  we  wish 
to  minimize  the  stock  rolls  used  and  yet  meet  customer  requirements,  we 
can  solve  the  problem, 
m 

Min.  T  d, 

iTi  ' 

m 

Subject  to  dj  ^  K'”  r  ss  1 . ,R 

i  I 

0  ^  di  4  1  i nteger 

Gilmore  and  Gomory,  and  Pierce  have  developed  special  techniques  for 
the  solution  of  large  cutting  stock  problems. 

Time  and  Interdependency  Models. 

The  articles  to  be  discussed  in  this  section  are  typical  of  the 
general  literature  In  the  class  of  interdependency  included.  In  terms 
of  our  definitions,  "job  alternative"  interdependency  Is  more  common 
than  "other"  Interdepends.icy.  It  is  true,  however,  that  in  the  models 
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with  job  alternatives  generated  by  probnbi I i Stic  research  or  production 
outcomes,  a  Job  whose  sole  predecessor  may  not  be  performed  may  be  con> 
sidered  to  be  contingent  on  that  predecessor.  This,  for  example,  is 
true  of  the  article  by  Eisner  ^281  discussed  above. 

Elmaghraby  C291  introduces  a  series  of  logical  relations  to 
standard  network  formulations.  Converting  his  work  to  our  notation,  we 
have  a  planning  unit  or  task  Sj  with  a  probability  of  occurrence  and 
a  vector  of  parameters  such  as  time,  cost,  etc.,  attached  to  it.  We 
now  define  logical  relations  that  may  exist  between  the  planning  units. 

1.  "and"  or  logical  intersection  of  two  events,  propositions  or 
activities.  For  example,  unit  will  occur  If  both  events  Sj 
and  S2  occur 

d3  4  bd|  +  (l-b)d2 

0  <  b  <  I 

d|  0-1  integer 

2.  "Inclusive-or"  i.e.,  the  union  of  two  or  more  propositions.  Node 
Sj  will  occur  if  S|  or  $2  or  both  occur 

dj  4  d]  +  d2 

3.  "exclusive  or",  often  referred  to  as  the  ring  sum.  S3  occurs  If 
either  $]  or  $2  but  not  both  occur. 

d|  4-  ^2  +<13  =  2 

4.  "decision"  node,  i.e.,  a  node  at  which  the  system  may  transfer 
along  one  path  or  the  other  with  known  probabilities,  Elmaghraby 
uses  this  node  solely  for  non-determini  Stic  branches. 

*^2  <^3  4  * 
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Note  thot  relations  I,  2,  and  3  would  be  Included  In  our 
definition  of  "other"  interdependency  and  4  would  be  considered  as  "Job 
alternative"  intf.rdependency.  A  graphical  symbol  is  defined  for  each 
interdependency  relation  so  that  a  set  of  tasks  and  the  rsistions 
between  them  may  be  expressed  graphically.  However,  if  the  nodes  are 

given  performance  times  and  the  . *’*'®^*®*’  **  Interpreted 

as  an  immediate  predecessor  relationship,  then  it  is  only  possible  to 
have  interdependency  relations  between  units  that  also  have  immediate 
predecessor- successor  relations  unless  new  symbols  are  defined.  This 
is  an  unnecessary  restriction  introduced  by  his  attempt  to  show  all 
relations  graphically.  As  we  snail  see,  a  programming  formulation  has 
no  such  restriction. 

Given  the  model  as  described  above,  Elmaghraby  suggests  a 
complete  enumeration  of  paths  and  shows  algebraically  that  for  each 
such  path  a  time  and  probability  of  occurrence  may  be  determined.  He 
concludes  by  combining  path  time  and  probability  information  for  an 
overall  expected  value  for  project  completion. 

The  problem  of  determining  a  project  cost  function,  that  is, 
total  cost  at  various  completion  dates,  for  the  deterministic  case  is 
discussed  by  Kelley  and  Walker  C431 .  The  length  of  the  project  may 
vary  because  for  each  job  there  is  a  series  of  "job  alternatives" 
(actually  a  continuous  linear  function)  with  increasing  cost.  C) ,  and 
decreasing  time,  tj .  The  operation  time  for  job  $|  is  co'istrained 
to  be  within  the  u^-per  time  limit  U|  and  the  lower  limit  L|  ,  that  is 
0  ^  L|  ^  t|  ^  Uj 
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•nd  cost  C)  s  b|  •  a| t| 

where  b|  Is  the  cost  of  Job  S{  performed  in  t!ir«  L; 
and  a|  is  the  cost  of  decreasing  tj  by  one  time  unit. 

Then  given  an  absolute  due  date,  0  ,  the  objective  function  is 

m 

Min  S  bj  -  a|t| 
isl 

Subject  to  precedence  constraints,  one  for  each  link  in  the  graph 

“l  +  'i  <  0 - -0 

and  finally  ^  0 

Kelley  refers  to  a  form  of  the  Ford-Fuikerson  algorithm  C331  for 
an  efficient  solution  to  the  problem.  Fulkerson  C351  presents  a  similar 
algorithm  in  full  detail.  Essentially,  he  interprets  the  problem  as 
one  of  network  rlow  and  solves  the  dual  of  this  problem.  The  algorithm 
begins  by  setting  all  Jobs  at  their  cheapest  (longest)  value  and  calcu> 
iating  the  project  length  that  results,  setting  D,  (due  date)  equal  to 
this  value,  and  criculating  the  cost,  P(0).  Then  D  is  decreased  and 
a  new  value  of  P(D)  as  well  as  W| ,  t{  for  all  jobs  is  calculated. 
The  process  is  continued  until  a  shortest  feasible  length  for  D  is 
obtained.  It  is  then  possible  to  plot  the  project  time-cost  curves. 


t 
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Fulkerson  points  out  that  the  breakpoints  of  this  piecewise  linear 
function  occur  at  integer  values  of  time  if  the  bounds  on  Job  lengths 
are  integers.  The  algorithm  has  been  shown  to  be  efficient  in  practical 
applications  to  large  problems  In  this  formulation  it  may  be 

assumed  that  the  variation  In  job  time  and  cost  is  a  result  of  the  appli¬ 
cation  of  more  or  less  resources  to  the  Job.  Since  there  Is  no  attempt 
to  restrict  the  amount  of  resource  used  at  a  particular  point  in  time, 
by  all  Jobs,  the  model  is  not  considered  to  have  resou'ce  constraints. 

In  some  instances  it  may  be  unrealistic  to  assume  a  continuous 
linear  relation  between  time  and  cost  for  a  task  in  a  project.  For 
example,  if  a  Job  may  be  only  performed  by  an  eight  man  crew  on  regular 
time,  cr  by  an  eight  man  crew  on  regular  time  plus  two  hours  overtime, 
there  are  two  discrete  ways  to  perform  the  Job  (two  "Job  alternatives") 
and  linear  combinations  of  these  methods  may  be  technologically  or 
contractually  infeasible.  The  problem  of  discrete  "Job  alternatives" 

In  the  time/cost  problem  is  discussed  In  references  [l 9,  26,  57i  5^. 

Moder  and  Philips  (38]  summarize  some  work  by  Meyer  end  Shaffer  (37^  on 
this  problem.  Their  formulation  is  as  follows: 

Min.  21  Cj j  djj 

'  j 

S.t.  precedence 

Wi  +  t]  ^Wj, 

or  if  a  Job  alternative  situation  is  involved 

W|  -f  til  d|i  4-  t|2  d|2 - +  l|k(l)  8ik(i)  <  Wj 

and  Wp  ^  0 
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'U 


■f  ■' 

^0  oth 


where 

Job  S(j  is  performed 
otherwi se 

This  essumes  that  each  Job  alternative 
Sij,  J*1 

has  identical  precedence  and  successor  relations. 

Finally,  we  require  an  interdependence  constraint 


HD 

E  dj  j  S  I 
j  « 


0  d;j  <  I  i 


Integer 


A  more  general  and  more  efficient  integer  programming  formula¬ 
tion  is  presented  by  Crowston  and  Thompson  [191,  1967.  This  model  will 
be  presented  in  detail  later,  but  a  short  summary  Is  included  here. 

The  Job  alternatives  are  again  represented  by  0-1  variables,  d|j, 
which  for  any  particular  Job,  arc  constrained 


In  addition,  however,  any  set  of  alternative  interdependence  constraints 
may  be  written  on  these  variables. 

‘^ij  +  ‘‘mn  ^  ’ 

^  I J  ^  ‘^mn 

<*ij  ^‘^mn 

etc. 

Note  that  although  resource  constraints  are  not  considered  specifically 
here,  alternative  Interdependency  constraints  could  be  written  to 
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constrain  the  total  usage  of  a  consumable  resource  over  the  life  of  the 
project. 

The  length  of  the  critica*  path  as  defined  in  Chapter  •  con¬ 
strained  by  a  set  of  equations  which  represent  paths  in  the  network. 
Paths  which  cannot  become  critical  may  be  dropped  from  the  problem. 
Rather  than  generate  a  time/cost  curve,  they  establish  a  due  date,  D, 
with  overtime  penalty  and  undertime  premium  and  solve  directly  for  the 
optimum  set  of  Jobs  to  be  performed.  It  would  be  possible,  however,  to 
solve  a  series  of  problems,  setting  progressively  tighter  upper  limits 
on  the  critical  path  (Wp^O)  and  thus  generate  a  time/cost  curve. 

The  article  also  gives  an  outline  of  a  heuristic  technique  for 
solving  problems  with  time  constraints  and  the  "Job  alternative"  type 
of  interdependency.  It  is  assumed  that  all  alternatives  for  a  given 
Job  have  identical  predecessor  -  successor  relations  and  that  the 
following  inequalities  hold: 

^i,l  <  ^i^  ••••  <^iyt(i) 

^  ijl  ^  i)2  •  •  •  ^  ijk  ( i ) 

The  routine  is  described  as  follows:* 

1.  Technologically  order  the  Jobs 

2.  Sat  each  decision  node  to  the  alternative  having  lowest 

cost. 

3.  Calculate  the  critical  path. 

4.  Reorder  by  Early  Start 


*  Crowston  and  Thompson  Cl9i,  pp.  20-21. 
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5.  Go  to  7 

6.  Recalculate  the  critical  path  starting  at  the  position  in  the 
ordered  Job  list  held  by  the  decision  node  of  step  (lO) 

7.  Identify  all  decision  nodes  in  the  critical  path 

8.  For  ail  the -nodes  of  step  (7)  calculate  the  net  reduction  in  total 
project  cost  achieved  by  substituting  the  more  costly  alternatives 

9.  If  no  alternative  reduces  overall  cost,  go  to  (12) 

10.  Find  the  alternative  that  gives  the  maximum  cost  reduction  and 
switch  the  relevant  decision  node  to  that  alternative 

11.  Go  to  step  (6) 

12.  Review  all  decision  nodes  that  were  previously  chosen  to  see  if 
sufficient  slack  has  been  generated  to  allow  the  reintroductlon  of 
a  longer  hut  cheaper  alternative.  If  no  such  opportunity  exists, 
go  to  (14) 

13.  Introduce  the  cheaper  alternative  found  in  step  (12).  Go  to  Step 

(12) 

14.  HALT 

The  two  small  problems  tested  by  this  routine  gave  the  optimum 
solution  altiiough,  as  the  authors  state,  it  will  not  always  do  so. 

A  routine  originally  published  in  the  0.0. D.  and  N.A.S.A,  Guide 
PERT  Cost  (26]  and  in  Alpert  and  Orkand  [2],  1962,  and  extended  in 
Hoder  and  Phillips  follows  a  somewhat  similar  routine.*  At  step  (8), 
however,  the  replacement  Job  chosen  Is  the  one  with  minimum  Incremental 
cost  per  day.  If  Job  Sjj  were  originally  chosen,  the  measure  would  be 

*  Moder  and  Philips  (561#  pp.  I09'I22. 
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C|m  -  Cjj 

where  Sj^  ?s  the  alterna- 
time  bein9  considered.  However,  if  this  would  cause  the  critical  path 
to  shift,  then,  rather  than  use  the  incremental  cost  criterion,  they 
choose  the  replacenisnt  Job  so  as  to  minimize  lC{|„  -  ^ij^' 
switch  is  followed  by  a  review  of  previously  selected  jobs  to  see  if 
Sufficient  slack  has  been  generated  to  allow  a  switch  back  to  an  original, 
cheaper  Job.  This  is  similar  to  steps  (12)  and  (13)  of  Crowston- 
Thompson.  As  the  process  continues,  the  cost  of  Jobs  chosen  increases 
and  the  project  length  decreases,  mapping  out  a  time/cost  curve,  but 
not  necessarily  the  optimum  one. 

Resource  and  interdependency  Constraints 

Models  in  this  category  are  essentially  knapsack  problems  with 
interdependency  constraints  added.  For  example,  Weingartner  [75,  76] 
adds  both  "Job  alternative"  and  "other** types  of  interdependencies  to 
the  Lorle-Savage  capital  budgeting  problem.  The  resource  constraints 
are  budget  limits  on  the  capital  expenditure  by  period.  Thus  in  each 
period  we  sum  the  capital  requirements  of  the  projects  to  be  operating 
in  that  period  and  constrain  the  total  amount  to  be  less  than  the 
budget  limit.  In  cur  notation  his  model  is 
m 

Maximize  ^  Cj  d{ 
isl 


where 


Cf  is  the  r>et  present  vaiue  of  project  S| 
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Subject  to 


2  kf  d|  =Ck'‘ 

txl 


r  s  ) . m 


where  K**  is  the  budget  limit  in  period  r  and  kf  is  the  cash  used 
by  project  S|  in  period  r.  The  d|  are  again  0-1  variables.  Any 
linear  constraints  in  the  variables  d{  may  be  written  to  express 
Interdependency  conditions.  The  model  may  then  be  solved  by  integer 
linear  programming. 

In  the  iatc;  paper  (75J  Welngartner  adds  to  this  modei  a 
time  dimension  by  allowing  a  given  project  to  be  represented  by  a 
mutually  exclusive  set  of  projects  (Job  alternatives),  one  beginning 
at  each  feasible  starting  day  of  the  project 

. This  follows  the  practice  of  Marglin  [563.  Even  though  time 

is  introduced  invo  the  model,  no  explicit  provision  is  made  for  time 
precedence  constraints,  a  natural  dimension  of  the  capital  budgeting 
problem.  In  this  article  a  new  and  reportedly  efficient  algorithm 
based  on  the  dynamic  programming  solution  to  the  knapsack  problem  is 
presented.  Unfortunately  this  method  will  not  handle  the  full  range 
of  possible  interdependencies  due  to  a  rescriction  on  the  inclusion  of 
negative  variables. 

A  vary  similar  problem  is  discussed  by  Root  [633  and 
te..ned  the  "selection  problem".  Explicitly,  the  problem  he  wishes  to 
solve  is 

m 

Minimize  ^  Cj  dj 
i»l 


m 

Subject  to  2]  ^  r»l....R 

i  I 
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where  the  K**  are  assumed  to  be  integer  and  a  set  of  linear  Interde¬ 
pendency  constraints  are  written  In  the  variables  d|  l■l,..._,m.  For 
example,  the  fact  that  a  Job  may  be  performed  by  several  resources  or 
resources  in  combination  may  be  written  as 


I 


where  the  kjj  would  differ  for  the  various  alternatives.  Root  solves 
the  problem  by  applying  a  theorem  from  symbolic  logic  to  reduce  the 
total  set  of  possible  solutions.  Then  by  costing  each  remaining 
solution,  he  can  select  the  one  with  minimum  cost. 


Time  and  Resource  Problems 

Many  in  ortant  scheduling  problems  may  be  described  as 
time  and  resource  constraint  problems.  Among  these  are  forms  of  the 
project  scheduling  problem,  the  Job-shop  scheduling  problem,  and  the 
assembly-line  balancing  problem.  The  project  scheduling  problem,  of 
course,  appears  in  contei.ts  such  as  marketing  and  economic  planning 
(311  as  well  as  production.  All  of  these  problems  have  been  formulated 
as  integer  linear  programming  problems  but  because  of  the  high  number 
of  constraints  involved,  these  are  not  suggested  as  possible  solution 
methods  for  real  problems.  For  example,  Wiest  180)  estimates  that  a 
project  with  55  Jobs  in  4  shops  with  a  time  span  of  30  days  would  have 
some  5,275  equations  and  1,650  varTables,  not  including  slack  variables 
or  constraints  added  to  assure  an  integer  solution.  As  a  result, 
heuristic  solution  methods  have  been  developed  for  these  problems. 

The  essential  problem  is  that  the  level  of  resources  is  constrained  by 


26 


period  and  the  Jobs,  given  the  usual  time  constraints,  may  shift 
through  time.  Thus, in  programming  formulations,  it  is  necessary  to 
include  the  possibility  that  the  job  may  shift  through  time,  and  this 
requires  many  variables  and  many  constraints.  Heuristic  solution 
techniques  can  handle  this  problem  with  concise  bookkeeping  techniques. 
The  solutions,  however,  are  not  necessarily  optimal. 

We  will  examine  several  heuristic  approaches  to  the 
resource- level  I ing  problem  in  project  scheduling.  It  is  assumed  that 
in  this  problem  the  resources  are  not  fixed  but  that  the  criterion 
function  is  so  related  to  period  by  period  resource  levels  that  we  are 
motivated  to  smooth  the  daily  resource  usage.  Burgess  and  Killlbrew 
Cl 2]  describe  an  iterative  procedure  that  attempts  to  minimize  the  sum 
of  squares  of  daily  resource  usages.  This  criterion,  while  minimizing 
the  standard  deviation  from  the  project  mean,  still  might  allow  high 
peak  in  any  given  time  period. 

The  routine  first  technologically  orders  the  Jobs  by 
Early  Start,  and  if  jobs  are  tied,  ranks  the  shortest  first.  In  the 
second  stage  the  Jobs  are  loaded,  beginning  at  the  bottom  of  the 
technological  list.  Each  Job  is  scheduled  as  late  as  possible,  subject 
to  the  condition  that  the  daily  usage  should  not  be  too  far  above  or 
below  the  predetermined  average  daily  resource  usage.  The  late  start 
of  each  Job  is  strictly  set  by  the  assigned  start  time  of  its  successors. 
The  cycle  is  repeated,  each  time  attempting  to  reduce  deviations  from 
the  mean,  until  no  further  improvements  are  made.  The  best  schedule 


is  then  chosen. 
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The  method  of  Levy,  Thompson  and  Wlest  [50j  concentrates 
on  reducing  the  peak  usage  of  the  project  resource.  Heuristics  which 
drive  the  solution  to  this  goal  would  also  work  to  meet  the  Burgess 
criterion.  The  jobs  from  all  projects  to  be  simultaneously  scheduled 
are  scheduled  at  Early  Start  and  then  the  daily  demand  for  ea'h  resource 
is  plotted.  For  the  first  resource  an  initial  trigger  level  is  set, 
one  unit  below  the  maximum  usage.  It  Is  assumed  that  the  resources 
are  ordered  based  on  some  priority  system,  perhaps  daily  cost  per  unit. 
Then  all  the  Jobs  contributing  to  the  particular  peak  and  in  addition 
having  enough  slack  so  that  they  may  be  scheduled  beyond  the  peak  are 
listed.  Then  one  of  the  Jobs  is  chosen  probabilistically,  the  proba¬ 
bilistic  weight  being  proportional  to  the  job's  slack  time,  and  the  Job 
is  shifted  a  random  amount  within  the  slack,  forward.  The  resource 
peak  again  is  calculated  and  a  lower  trigger  level  set.  This  process 
continues  for  the  first  resource  until  no  further  improvement  is 
realized.  Then  freezing  the  lowest  feasible  limit  for  the  first 
resource,  the  procedure  is  repeated  for  the  second  resource  and  so  on 
through  all  resources.  Now  the  jobs  from  each  project  are  segregated 
and  again  an  attempt  is  made  to  shift  them  and  reduce  the  aggregate 
trigger  level  for  all  projects.  Finally,  when  no  further  i nprovement 
is  possible,  the  final  schedule  and  trigger  levels  ate  stored  and  the 
process  is  completely  repeated.  Due  to  the  probabilistic  element  In 
the  decision  rule,  a  new  schedule  will  result.  After  several  repeti¬ 
tions,  the  best  schedule  of  those  generated  can  be  chosen. 

The  next  problem  class  to  be  discussed  will  be  scheduling 
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to  meet  stated  resource  constraints.  The  early  solutions  to  problems 
of  this  type  were  obtained  with  Gantt  charts  and  their  use  In  Job-shop 
scheduling  continues.  The  resource  constraint  in  the  job-shop  problem 
will  be  the  limit  on  machine  availability  and  the  time  constraints  are 
provided  by  technological  ordering  of  operations  on  a  particular  work 
order.  The  criterion  function  may  be  to  minimize  completion  time  of  the 
whole  job  file  as  in  the  integer  programming  formulations  of  Bowman  181, 
Manne  fS^'I  ,  and  Wagner  (73l .  Alternately,  it  might  be  the  minimization 
of  idle  machine  time  (identical  to  minimim  final  completion  time  for 
the  fixed  job  file  case)  as  in  Conway  and  Maxwell  Ll8l  or  a  complex 
function  of  order  delay  cost  as  in  Carroll  rj4]. 

Many  heuristic  approaches  have  been  developed  for  this 
problem  (18,  20,  381  and  in  many  of  these  approaches  similar  decision 
rules  are  used  singly  or  in  probabilistic  combination  to  dispatch  jobs 
from  a  queue  to  the  machine.  We  will  quote  from  a  description  of 
several  such  rules:  * 

1.  SIO  (shortest  imminent  operation) 

When  a  facility  is  available,  select  that  item  in  the  queue 
which  has  the  shortest  machining  time  on  the  facility. 

2.  LRT  (longest  remaining  time) 

Select  the  I  tern  which  has  the  most  total  machining  time 
remaini ng. 

3.  J.S.  (job  slack  per  operation  remaining) 


*  Crows  ton.  Glover,  Thompson  and  Trawl ck  tZOl ,  p.  2. 
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Subtract  the  total  remaining  machine  tire  for  a  given  Job 
from  an  arbitrary  finish  date,  DO.  Divide  this  slack  by 
the  number  of  operations  remaining.  Select  the  I  tern  with 
minimum  Job  slack  per  operation. 

4.  Lie  (longest  imminent  operation) 

5.  FIFO  (first  in,  f i rst  out) 

Select  the  item  that  arrived  first  in  the  queue. 

6.  MS  (machine  slack) 

For  each  machine  calculate  the  total  machining  time  remaining. 
Select  the  Job  in  the  queue  that  goes  to  the  most  heavily 
laden  machine  next.  Break  ties  with  S.I.O. 

Many  other  rules  have  been  attempted  to  specifically  meet 
more  complicated  criterion  than  the  minimization  of  overall  completion 
time.  It  will  be  interesting  to  note  below  that  attempts  have  been 
made  to  apply  several  of  the  simple  rules  to  the  project  scheduling 
problem,  but  that  more  complex  rules  have  not  as  yet  been  translated 
to  the  project  problem. 

Two  main  approaches  have  been  su  :ed  for  solving  the 
fixed  resource  problem  in  the  project  scheduling  context.  Neither 
approach  guarantees  an  optimum  solution,  nor  consistently  outperforms 
the  other.  Representative  of  the  first  approach  is  an  article  by 
Kelley  i.44l  .  The  routin'*  may  be  summarized  as  follows: 

I.  Technologically  order  the  Jobs,  calculate  early  start, 
late  start,  and  slack.  Within  the  technological  ordering,  reorder  by 
slack,  lowest  slack  first. 
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2.  Start  with  the  first  activity  and  continue  down  the 

list. 

(a)  Find  the  eariy  start  of  the  activity. 

(b)  Schedule  the  activity  at  early  start  if  sufficient 
resources  are  available.  If  resources  are  not  avail¬ 
able,  two  alternative  routines  are  suggested  cl,  c2. 

(cl)  Serial  Method:  Begin  the  Jab  at  the  earliest  time 

that  resources  are  available  to  work  it  for  one  day. 
The  job  may  be  split  if  necessary. 

(c2)  Parallel  Method:  Find  the  set  of  all  jobs  causing 
the  resource  violation  and  rank  them  by  total  slack. 
Delay  a  sufficient  number  of  jobs  with  the  highest 
slack  so  that  those  remaining  on  the  list  can  be 
scheduled. 

3.  Repeat  step  2  for  all  jobs. 

4.  Repeat  the  whole  routine  with  various  orderings  on  the 
technological  list. 

In  addition  to  the  possibility  of  allowing  job-splits, 
Kelley  suggests  that  we  also  allow  resource  limits  to  be  violated  by 
small  amounts.  The  basic  heuristic  in  the  original  job  ordering  and  In 
the  parallel  loading  technique  is  the  use  of  job-slack  (J3)  as  a 
criterion  for  shifting  jobs  forward.  This  is  a  reasonable  approach 
similar  to  that  used  by  Levy  rsQLl. 

The  second  approach,  detailed  in  Moder  and  Philips  (S8I 
uses  Maximum  Remaining  Path  Length,  or  equivalently  Late  Start,  as  a 
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basis  for  delaying  jobs  (LRT).  Note  that  If  two  jobs  have  a  conmon 
early  start-time,  both  measures  are  equivalent.  That  is,  the  job 
delayed  because  of  a  higher  slack  would  In  the  same  way  be  delayed 
because  of  a  lower  Late  Start  time.  Some  other  features  of  the  routine 
are  different  than  that  suggested  by  Kelley.  A  main  feature  Is  a  list 
of  unscheduled  jobs  whose  predecessors  have  all  been  scheduled,  ordered 
by  Late  Start,  and  an  ordered  list  of  finish  times  of  scheduled  jobs. 

Thus  the  routine  steps  through  time  to  only  those  days  when  a  change 
of  resource  usage  is  possible.  On  such  a  day  it  examines  the  list  of 
available  jobs  and  ends  either  when  the  resources  are  exhausted  or  the 
job  file  is  completed.  At  this  point  it  again  jumps  ahead.  As  we 
implied  at  the  beginning  of  this  section,  examples  can  be  constructed 
CO  favor  or  penalize  either  of  these  approaches. 

Finally,  we  wi 1 1  show  that  Salveson's  [sO  formulation  of 
the  assembly-line  balancing  problem  has  the  structure  of  a  time- re sources 
constraint  problem.  The  resources,  in  this  instance,  are  the  work 
stations  and  the  work  times  at  the  station  will  be  the  resource  level. 
Then,  if  a  precedence  ordering  exists  between  jobs,  as  would  usually 
be  the  case,  a  job  may  not  be  assigned  to  a  work  station  unless  all  its 
predecessors  have  previously  been  assioned  to  that  work  station  or 
earlier  work  stations.  In  addition  to  the  above  statement,  the  problem 
may  be  complicated  by  any  of  the  interdependencies  we  have  discussed. 

For  example.  It  may  not  be  feasible  to  do  a  particular  pair  of  jobs  at 
one  station,  and  therefore,  we  must  constrain  the.  problem  to  prohibit 
this.  Also,  it  should  be  pointed  out  that  "alternative  interdependency" 
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relations  would  probably  be  the  most  efficient  method  of  actually 
formulating  the  problem.  Finally,  as  in  the  other  problems  of  this 
section,  it  must  be  noted  that  practical  problems  must  be  solved  by 
heuristic  methods  (n). 

Models  with  Resource.  Time  and  Interdependency  Constraints. 

Problems  in  this  category  abound  in  the  industrial  world. 
Wherever  there  is  an  element  of  physical  design  connected  with  a 
scheduling  problem,  then  interdependency  is  an  implicit  part  of  the 
problem.  For  example,  we  may  consider  the  problem  of  an  industrial 
engineer  attempting  to  design  a  job  with  the  graphical  technique  of 
man-machine  analysis  [4].  There  are  several  alternate  ways  for  an 
operator  to  perform  each  necessary  ^unction  and  operations  when  combined 
may  require  less  time  than  the  sum  of  the  same  operations  performed 
singly.  Finally,  many  possible  orderings  may  be  possible,  although 
sc<me  technological  ordering  constraints  are  involved  in  many  problems. 
Each  production  alternative  may  have  an  important  effect  on  the  sched> 
uling  problem.  Similarly,  design  of  construction  projects  will  have 
important  scheduling  Implications. 

With  the  exception  of  the  Weingartner  article  I75J  «vhich 
does  not  fully  Include  the  possibility  of  time  ordering,  no  model  was 
found  to  include  resource,  time  and  both  types  of  interdependency  con¬ 
straints.  It  Is  true,  however,  that  the  job-shop  models  of  Bowman  181 
and  Wagner  C7J)  could  be  easily  generalized  to  Include  this  possibility. 
The  use  of  job  alternative  Interdependency  in  heuristic  programs  Is  not 


uncommon. 
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Two  master's  theses  written  in  the  Sloan  School  of  Manage¬ 
ment,  M.I«T.,  expand  the  usual  statement  of  the  Job-shop  scheduling 
problem  to  include  the  possibility  of  Job  alternatives,  with  two  differ¬ 
ent  practical  interpretations  of  what  these  alternatives  might  be. 

First,  a  thesis  by  Russo  [.65I  reports  on  a  simulation  of  a  Job-shop 
with  the  possibility  of  alternate  Job  routing.  In  the  usual  simulation 
we  assume  for  a  given  Job  that  the  order  of  operation  isstrictly  deter¬ 
mined,  that  is 

0 - -0 - <5) - ® 

He  suggests  that  in  some  instances  the  ordering  would  be 


and  interprets  and  solves  this  as  a  case  of 


-  ‘^3.1 
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It  Is  clear  that  this  s!tu^itio^  can  be  categorized  as  a  time  and  resource 
model,  however,  since  his  heuristics  specifically  detailed  each  alterna¬ 
tive  route,  then  calculated  Information  for  local  dispatching  rules 
based  on  the  alternatives  and  then  decided  between  them,  it  is 
included  here.  One  of  the  most  effective  approaclies  he  discusses,  allows 
all  Jobs  in  an  alternate  chain  to  enter  queues  for  their  respective 
machines  when  the  predecessor  of  the  alternate  chain  is  completed.  Then 
in  each  queue  the  Jobs  from  the  alternate  chain  are  ranked  by  a  standard 
dispatching  rule.  Me  then  allows  that  Job  to  go  first,  which  is  selected 
by  the  machine  queue  discipline. 

Clernwnt  C!71  allowed  for  the  possibility  of  sevcr*sl 
machines  performing  a  given  job.  This  more  clearly  resembles  our  Job 
alternative  interdependency  case.  His  results  show  that  the  heuristic 
Russo  found  so  successful  only  improved  performance  of  the  simple  dis¬ 
patching  rules  wlten  the  total  machine  ioadings  were  highly  imbalanced. 

In  the  case  of  balanced  loads,  switching  actually  decreased  the  perform¬ 
ance  of  simple  decision.  For  the  balanced  case  only  the  simple  switch¬ 
ing  heuristic  “switch  if  the  alternate  queue  is  empty"  consistently 
improved  performance  of  the  simple  dispatching  rules.  His  results 
also  Indicate  that  a  dispatching  rule  “Covert'  derived  by  Carroll  Cl4] 
completely  dominates  all  conventional  rules  tested. 

Job  alternativi-  interdependency  Is  also  found  In  scheduling 
techniques  designed  for  the  1 srge  project  problem.  The  series  of  SPAR 
programs  by  Wiest  f80l  wi ' 1  now  be  discussed.  With  each  Job  he  essoci- 
ates  three  operating  levels,  maximum  crew  size,  norma)  crew  size  and 
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minimum  crew  size.  Of  course,  the  '  length  Is  inversely  a  function 
of  resource  level.  The  jobs  originally  at  normal  crew  size  are  origin¬ 
ally  ordered  by  early  start  time,  one  possible  technological  order.  Then, 
as  in  the  Hoder  and  Philip  routine  discussed  above,  a  sub-list  is  gen¬ 
erated  and  continually  updated  of  jobs  available  for  scheduling.  From 
this  list  jobs  are  selected  for  schedul’ng  with  a  probability  inversely 
related  to  the  available  slack  f  the  job.  if  a  job  is  selected  to  be 
scheduled  and  the  resources  are  unavailable,  it  is  left  to  be  schedule 
in  a  subsequent  period. 

Several  subsidiary  heuristic  routines  operate  within  the 
basic  framework  outlined  above.  If  the  slack  of  a  job  is  low,  an 
attempt  is  made  to  schedule  it  at  maximum  resource  usage.  If  the 
resources  are  not  available  for  this,  a  subroutine  attempts  to  borrow 
resources  from  jobs  operating  on  that  day  with  normal  or  maximum 
resource.  A  second  approach  is  to  find  jobs  using  the  tight  resource 
and  delay  their  start  for  one  or  more  periods.  This  frees  their 
resources  for  the  critical  jobs.  Finally,  if  ail  else  fails,  the 
critical  job  would  be  delayed  one  period.  After  the  application  of 
these  and  other  routines,  a  final  schedule  is  produced.  Since  there 
is  a  random  element  in  the  choice  of  jobs  to  schedule,  it  is  suggested 
that  the  process  be  repeated  several  times  and  the  best  schedule 
selected. 

An  important  addition  to  the  program  is  a  search  routine 
that  progressively  shifts  the  level  of  initial  resources  from  solution 
to  solution  in  attempt  to  find  those  limits  that  will  minimize  the  sum 


of  re  ource  cost  and  project  completion  time  cost. 


Discussion  of  Constraint  Categories 

A«  we  have  seen  the  categories  chosen  do  not  allow  a 
categorization  of  all  planning  models.  In  several  ca^es  a  good 
arniiTient  could  L'e  made  against  the  allocation  we  have  decided  on. 
Nevertheless,  for  several  reasons,  this  particular  categorization  is 
useful.  It  does  suggest  that  more  complex  models  may  have  relevance 
in  functional  areas  where  they  have  not  yet  appeared.  For  example, 
time  precedence  constraints  would  seem  relevant  to  the  capital  budget 
ing  problem,  further  use  of  alternative  interdependencies  the  job- 
shop  problem  and  the  assembly-li ne  problem. 

An  even  more  fruitful  use  of  such  a  framework  will  be  to 
suggest  that  researchets  examine  a  broad  range  of  literature  in  their 
search  for  suitable  solution  techniques.  Examples  of  this  certainly 
have  already  occirred.  For  example,  Wlest  LboJ  uses  the  Bowman  [Sj 
job-shop  L.F.  formulation  as  the  base  of  his  project  scheduling  formu 
lation.  A  thesis  oy  Kr’iht  L^b]  attempts  to  relate  heuristics  in  the 
job-shop  problem  to  ht:uri  sties  in  the  project  scheuulinq  problem. 
Finally,  Wilson  [82]  shows  how  several  line  balancing  algorithms  may 
be  applied  to  resource  leveling.  Tpir.  compa'-ison  also  allows  us  to 
make  some  general  statements  about  the  efficiency  of  various  sdiution 


techniques . 
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Solution  of  Planning  Models 

For  several  categories  of  problems,  algorithms  have  been 
developed  that  are  much  more  efficient  than  any  of  the  existing  program¬ 
ming  routines.  For  example,  the  longest  path  algorithm  very  quickly 
determines  the  critical  path  of  a  simple  time  constraint  problem. 
Similarly,  the  Ford-Fulkerson  technique  efficiently  finds  the  optimum 
Job  lengths  and  job  start  times  for  a  project  given  a  fixed  due  date 
and  a  bounded  inverse  linear'  relation  between  job  cost  and  job  time. 

This,  of  course,  is  a  special  case  of  tl>e  time-interdependency  constraint 
case. 

There  is  a  second  set  of  problems  for  which  special 
efficient  solution  techniques  have  been  developed,  but  because  of  the 
nature  of  the  problem,  the  techniques  may  be  considered  as  restricted 
integer  programming  routines.  In  this  category  I  would  include  the 
truncated  enumeration  methods  of  Weingartner  and  Ness  and  that  of 
Glover  applied  to  the  simple  resource  (or  knapsack)  problem.  Also 
Root's  algorithm  for  the  selection  problem,  which  is  a  combined 
resource- i nterdependency  problem,  is  in  this  group.  We  may  summarize 
the  above  cases,  then,  by  saying  that  we  may  obtain  optimum  solutions 
more  efficiently  with  existing  techniques  than  with  programming 
methods . 

We  now  shift  in  the  spectrum  to  those  problems  for  which 
opthnin  solutions  are  only  avaiioble  through  programming  methods.  It 
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should  be  emphasized  that  no  clear  line  can  be  drawn  between  problems 
in  this  group  and  those  in  the  previous  one.  An  example  of  this  would 
be  problems  with  pure  interdependency  constraints,  ^e  have  stated 
above  that  certain  simple  resource  or  resource- interdependency  con¬ 
straint  problems  ma  '  be  solved  with  special  computational  methods.  On 
the  other  hand,  if  many  resources  were  involved,  it  would  be  necessary  to 
turn  to  conventional  integer  programming  techniques.  Similarly,  for 
many  large  problems  with  pure  interdependency  constraints,  the  0-1  tree 
search  algorithms  wouid  give  the  most  efficient  solutions.  Finally, 
the  time-i nterdependency  problem, when  the  job  alternatives  are  discrete, 
is  best  solved  by  conventional  integer  programming  routines  [70]. 
Techniques  are  available,  however,  for  this  problem,  as  we  will  show 
later,  that  will  significantly  reduce  the  number  of  required  constraints 
and  the  number  of  variables.  In  this  regard,  it  is  interesting  to 
observe  that  for  integer  programming  problems  In  0-1  variables,  added 
interdependency  constraints,  by  eliminating  branches  in  the  feasible 
solution  tree,  actually  make  problems  easier  to  solve. 

Heuristic  routines  are  available  for  the  above  problems, 
and  we  have  covered  several  approaches  to  the  time- interdependency 
problem,  given  discrete  Jobs,  no  alternative  interdependency  and  common 
precedence-successor  relations  for  Job  alternatives.  It  can  be  shown 
that  existing  heuristic  methods  are  not  at  all  appropriate  for  problems 
with  any  reasonable  alternative  interdependency  complications. 

The  Final  category  finds  combinations  of  time  and  resource 
constraints  and,  as  we  have  discussed,  this  grouping  requires  large 
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numbers  of  constraints  and  variables  for  programming  solutions.  As  a 
result,  these  problems  are  solved  almost  exclusively  in  practice  by 
heuristic  methods.  As  noted  above  some  of  the  simpler  Job-shop  problem 
heuristics  have  been  adapted  to  the  project  scheduling  problem  with 
some  success.  This  suggests  that  more  complex  heuristics,  based  on 
variations  of  successful  Job-shop  rules,  should  now  be  tested.  When 
Job  alternative  interdependencies  are  added  to  the  problem,  as  Wiest 
[801  does,  it  is  possible  to  build  subsidiary  switching  routines  based 
on  local  (in  the  time  dimension)  resource  usage  information.  If  alter¬ 
native  interdependencies  were  added  to  the  problem,  such  local  informa¬ 
tion  might  no  longer  be  a  sufficient  base  for  a  switching  decision.  In 
fact,  as  we  have  seen,  any  purely  heuristic  method  might  have  difficulty 
approaching  an  optimum  solution  to  this  kind  of  problem. 


Chapter  ill 
DECISION  CPM  MODELS 


The  paper  hat  will  serve  as  the  basis*  for  this  chapter  is 
that  of  reference  TiSI  by  Crowston  and  Thompson.  This  decision  network 
formu'iation  contains  "time*'  constraints  and  both  "Job  alternative"  and 
"other"  types  of  Interdependency.  The  mathematical  basis  of  Decision 
CPM  will  be  discussed  as  well  as  several  alternate  integer  programming 
formulations  of  the  problem.  A  numerical  problem  is  introduced  in  this 
chapter  (FIGURE  lll>l)  that  will  also  serve  as  an  example  in  Chapters 
IV  and  VI. 

2.  The  Mathematical  Basis  of  Decision  CPM 
This  section  will  follow  in  part  the  article  t2J  by  Levy, 
Thompson,  and  Wiest.  Let  Jr  Ol.  52.  S3  ...)  be  a  set  of  igb  sets 
that  must  be  done  to  complete  a  project.  Some  job  sets  are  unit  sets 
$1  -  {^iljf  ^nd  other  job  sets  have  several  members,  Sj  *  5j2, 

Sjj,  ,..J.  In  order  to  complete  the  project,  one  of  the  jobs  from  each 
job  set  must  be  completed.  Associate  with  each  job  set 

^i  “  i?l  1 . ^ik(i)} 

k(i)  variables 

*Sectlons  2,  3,  ^  and  5  are  taken  essentially  verbatim  from 
fi9j  although  a  new  numerical  problem  is  introduced.  Furthermore,  this 
chapter  will  assume  that  exactly  one  job  alternative  is  selected  from 
each  Job  set ,  that  i s 
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Note  that  (a)  says  that  we  cannot  do  both  S{j  and  ;  (b) 

says  that  we  can  only  do  Sjj  if  we  also  do  ;  and  (c)  says 

that  we  either  do  S{j  and  or  v«  do  neither.  The  above 

discussion  illustrates  sotne  of  the  possible  complexity  of  problem 
formulation  that  is  possible  within  the  Decision  CPH  framework. 

In  addition  to  the  relations  described  above,  there  will  be 
precedence  relations  between  the  jobs  of  a  decision  project.  Let  ' 
denote  a  relation  between  pairs  of  jobs  in  J  such  that  S|j  ^ 
is  defined  for  some  pair  of  jobs  Sfj,  and  is  read  S]j  i s  an 

immediate  predecessor  of  The  Interpretation  of  this  statement  is 

that  ail  immediate  predecessors  of  a  job  must  be  completed  before  that 
job  can  be  started.  A  decision  project  is  the  set  J  together  with  the 
specified  interdependencies  and  the  relation  ^  defined  on  J. 

The  decision  project  graph  of  a  project,  G,  is  a  planar  graph 
with  nodes  representing  jobs  and  a  directed  line  segment,  connecting 
two  nodes  S]j,  it  and  only  if  Sjj  holds.  A  path  in  G 

is  a  set  of  nodes  connected  by  immediate  predecessor  relations,  a  cycle 
in  G  is  a  closed  path  of  tt«  form  Sjj  s  •!  4C.  a2  ^  •  •  •  ^^n”  ~ 

Sjj.  A  project  graph  is  acyclic  if  and  only  If  it  has  no  cycles. 
Definition:  Sjj  <  implies  Sjj  precedes  5^,^,  (or  alternatively 

succeeds  S|j)  if  and  only  if  there  is  a  set  of  jobs  aj,  a2  •••0^ 
n  2  such  that 

SiJ=  «2<^*3  •  •  •  «an= 


in  other  words,  Sjj  precedes 


if  and  only  if  theis  is  a  path  from 
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Sjj  to  In  the  decision  project  graph  G. 

Assumption  I  The  precedes  relation  Is  asymmetric,  that  Is  If  S|j 
then  It  is  false  that  for  all  Sjj  and  S^n  In  J.  Defini¬ 

tion:  A  relation  that  is  transitive  and  asymmetric  Is  said  to  be  a 
orefe rence  relation. 

Theoretm'f*  1 1 1-1:  If  assumption  1  holds,  then  the  predecessor  relation 
is  a  preference  relation,  and  the  graph  G  is  acyclic. 

Definition;  A  technologically  ordered  Job  list  J=  |a] ,  32-..  an) 
is  obtained  from  a  set  of  Jobs  J  » -^a,  b,  c  .  .  by  listing  them  so 
that  no  Job  appears  on  the  list  until  all  of  its  predecessors  have 
already  appeared. 


Theorem'*'?*  1 1 1-2:  Assumption  I  holds  If  and  only  if  it  is  possible  to 
list  the  job-  in  J  i n  a  technologically  ordered  Job  list  J. 

In  addition  to  these  definitions  and  theorems  from  reference  M. 
several  additional  conventions  are  necessary  because  of  the  fact  that 
some  Jobs  may  be  eliminated  from  the  decision  project  graph  as  the 
result  of  decisions  that  are  made,  if  we  decide  to  do  one  of  the  Jobs 
in  a  job  set,  then  all  immediate  predecessor  relations  that  the  Job 
satisfies  must  hold  in  the  flno!  graph.  If  vie  decide  not  to  do  that  job, 


*Note  that  this  assumption  differs  from  the  correspotidi ng 
assumption  in  (49j  in  that  the  requirement  of  K- i ntransi ti vi ty  is 
omitted.  For  this  reason  theorem  I  of  that  reference  does  not  hold 
in  the  prese.:t  context. 

**The  proofs  of  these  theorems  are  exactly  as  in  reference  [9^. 
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then  none  of  Its  Imnediate  predecessor  relations  hold.  In  the  decision 
project  graph,  if  we  decide  not  to  do  a  given  Job,  then  we  must  remove 
that  Job  togecher  with  ail  edges  that  impinge  on  it  from  the  decision 
project  graph  to  obtain  the  final  project  graph.  It  follows  from  this 
that  if  any  Job,.  Sjj  has  a  sole  immediate  predecessor  S^i^,  and  if 
that  predecessor  is  a  member  of  a  Job  set,  it  will  be  necessary  to 
create  a  dummy  immediate  predecessor  relation  between  S|j  and  a  Job 
which  is  a  predecessor  of  If  this  is  not  done,  then  it  would  be 

possible  for  the  path  containing  S;j  to  be  broken  and  Sjj  would 
lose  its  project  time  ordering.  Sim! lari 1y  a  dummy  immediate  successor 
relation  must  be  established  for  Jobs  having  only  one  immediate  succes¬ 
sor,  If  that  successor  is  a  member  of  a  Job  set.  in  addition  it  may  be 
necessary  to  create  a  dunny  relation  between  two  Jobs  even  if  both  have 
several  immediate  predecessors  and  successors.  If  on  any  path,  two 
Jobs  are  separated  by  a  Job  which  could  be  elimikiated,  and  if  it  is 
desired  to  maintain  a  technological  ordering  of  the  two  Jobs,  a  dummy 
immediate  predecessor  relation  must  be  established  between  them. 

For  a  given  project,  when  the  Jobs  are  technologically  ordered 
and  ail  planning  decisions  are  made  designating  the  Jobs  to  be  performed 
in  each  set,  the  nonnal  critical  path  analysis  may  then  be  carried  out. 
The  usual  concepts  of  early  start,  late  start,  critical  path,  etc., 
will  apply  to  this  reduced  graph.  These  terms  will  bo  used  throughout 
the  balance  of  the  thesis  without  further  definition. 
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3.  Decision  Project  Graphs 

A  graphical  representation  of  the  combined  planning  and 
scheduling  problem  is  shown  in  the  decision  project  graph  of  Figure 
lll-l.  In  this  graph  the  circular  nodes  represent  Jobs  and  the  triangu¬ 
lar  node-.'  introduce  the  mutually  exclusive  decision  Job  nodes  of  a  Job 
set.  In  addition  to  the  precedence  constraints  implied  by  the  directed 
line  segments,  we  impose  several  "alternative  Interdependencies". 

These  are 

(a)  dg  j  4  dg  2 

^6.2  —  ‘^12,1 

where  (a)  says  that  Sgj  may  only  be  performed  if  ^5  2  performed; 
and  (b)  says  that  2  ^j2  1  ^  performed  or  both 

not  be  performed.  The  problem,  then,  is  to  select  the  project  graph 
which  minimizes  total  project  cost. 

4.  Decision  Graph  Solution  by  Integer  Progr-vnminq 

Consider  job  set  (l)  and  its  associated  decision  vari¬ 
ables  with  constraints  given  by  (2)  (3),  and  (4).  Sesides  these, 

there  may  be  any  ’■he  other  constraints  descussed  in  the  second 
section  of  this  chapter  or  other  constraints  showing  various  types  of 
complicated  interdependencies  between  jobs  in  the  project. 

Witii  each  Job  S|,  we  associate  a  time,  t>  ,  and  a  cost. 

Cjj.  Also,  we  assume  a  reward  of  'r*  dollars  per  day  for  each  day  under 
D,  the  required  due  date  of  the  project,  and  a  penalty  payment  'p*  for 
each  day  beyond  D.  As  defined  earlier,  Wf  will  be  the  number  of 
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days  after  0  that  the  project  finishes  and  Wf"  will  be  the  number 
of  days  before  D  that  the  project  finishes.  We  can  now  formulate  the 
integer  prjgramming  problem  of  selecting  that  best  project  graph  and 
finding  its  critical  path 

H  k(i)  ^ 

S  ‘^ij^^ij  ■  ^  ^f~  +  P 
i=l  j=l 


The  first  term  calculates  the  cost  of  all  the  decision  jobs 
or  job  alternatives  that  are  to  be  performed.  It  is  governed  by  the 
constrai nts 


0  ^d.  .  <.1 

•J 


'  ®  J  »2 . m 


where  djj  i s  an  integer.  The  second  term  is  explained  by  the 
constrai nt 


+■ 

Wf  -  Wf  -f  Wf  -  0  =t  0 

where  is  the  early  start  time  of  Finish,  the  last  job  in  the  pro¬ 
ject.  If  D,  then  the  project  is  not  completed  until  after  the 

due  date  so  that  “  Uf  -  D  and  a  penalty  of  pW^  is  incurred. 

Other  constraints  must  hold  because  of  the  precedence 
relations.  For  instance,  if  S;  and  are  unit  set  jobs  and 

S|«S„.  we  have 


W.  -F  t;  4  W„ 

1  »  fn 


where  W;  is  the  early  start  time  cf  job  Sj  .  If  Is  a  unit  job 

set  j  and  Sjj  is  from  a  multi-job  set  and  Sjj  <<^5^^,  then 
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-M(l  -  djj)  -|"  W.j  +  t,j 

where  M  is  u  large  enough  number  so  that  the  inequality  does  not 
constrain  the  variable  d|j  unless  d|j  «  1.  Thus  all  paths  through 
Jobs  that  are  not  performed  will  be  broken. 

It  is  i>ow  possible  to  set  up  the  problem  of  .^igure  lll-i  as 
an  integer  programming  problem.  Note  that  it  is  not  necessary  to  include 
the  cost  of  unit  set  Jobs  in  the  functional  or  the  finish  node  (S^) 
in  the  precedence  constraints. 

Min.  400  dg  ,+’00  dg  2  3  dgj+IOO  dg  2 

■+0dc.  3+50d,2  ,+0d,2^2''‘'°°‘*15,f‘'°^l5,2‘*°^J7,i 

■+  _ 

+'l50d,^^2'^P'*^f  ” 

St. 

Precedence  Constraints  (link  formulation) 


<  W7 

*2 

^  % 

'3 

4^6,1 

<  V2 

<  '^6,3 

-M(|.d6j)+t6,,tW6,; 

t?  + 

-Md-d^ 

w,o 

-M(l-d^  j)+t^  j 

<  “10 

-Hd-d^  3)+t^  3 

t44-'^4  ^  *^8 

^8  ■^'^8  ^*^10 

ts+Wfi  <.  ^9,2 
t5+^5  <‘‘^9.3 

t8+'^8  ^  '^9,1 

^8'*'  '^8  ■^'^13 

-Md-dg  jHtg^l+^a.l  ^  ^3 
-M(l-d9^2^+^,Z^*^9.2  ^ 

-Hd-d,  <•  “l4 

-HCl-d,  ,Ht5_jt«9,3  <  ‘'15.1 

-Hd-d^  jj+tj  j+Wj  j  <  W,5_2 

‘io+''io  ^“n 

'ii+^n  ^ 

*10+*'l0  ^  “l2.' 

^  '^12,2 

tj4+Wj4  ^16 

-Hd-d,j_,Ht|j_|+W,j,,  ^  “it 

-Hd-d,,  2M-‘|2,2'"u,?  ^  ‘'l* 

-Hd-d,5,,Ht,5_,TW,5_,  4:  K|e 

-Hd-d,5.2Ht|5,2^*'l5.2  <  “W.l 
.Hd-d,5  2>-“|5,2  ‘'15.2  ^  "n-J 

.Hd-d,2_|)+t|2  |t<'|7.|  <  “18 
.Hd-d|7  j)+t|7  2+''|7,2  ^  “18 
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■+ 

4- 

-h(l-d,5  2)+ti5,2 

4  Wf 

4  Wf 

ti8  +W18 

4wf 

Due  Date  W^.  - 

i-w^-  D 

Interdependence 

Job  A| ternati ves 

+‘‘6,2  +‘‘6.3  = 
‘‘9.1  ^9.2  ^ 

‘*12,1+  ‘‘12.2  “ 

‘*15,1  ‘‘15,2 

‘‘l7.2  +  ‘‘l7.2  * 


Other  Interdependency 


'9.1 


6.2 


*'6,2  ^12,2 

0  <-d,j  4  > 

An  alte'’n«te  formulation  suggested  in  1191  and  the  formulation 

to  be  extended  in  this  chapter,  constrains  the  length  of  the  critical 

path  by  including  one  constraint  for  each  possible  path  from  to 

S^.  A  typical  path  might  pass  through  Sc,  S,,  ....  S,,.  S  and 

f  ■*  '  I J  mn 

S^.  The  constraint  for  this  path  would  be  written 
t,+  ...  -H(l-d,j)  +  ,,j  ...  ^  W, 
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The  precedence  constraints  on  the  length  of  the  critical  path 
for  Figure  1 1 1-1  will  now  be  shown  In  "path"  form. 

Precedence  (paths) . 


titt7+tio+tii-M(l-di2j)f-ti2j+t,g 

ti+t7+tio-M(l"di2  j)+t,2j+t,6 
tl  •ft7+tl0"^(>-<li2,2)+^12,2+^16 

ti-M(l-dg  i)+tg  i+tj,+  tio+tij-M(l-di2  ,)+ti2  itti6 
ti-M(l-d^  ,)i-t6j4-t7tt|o-M(l-di2j)+t,2j-Ki6 
tj-Md-dg  jKtg  j-rt7+t,o-M(l-d|2^2^+‘^12,2'*’M6 
t]-M(l-dg  2)‘»'t6, 2'‘’^10-‘-^ir^^* ■‘^12,1  J+^1 2,1  + *^16 
ti-M(l-dg  2)^t6, 2+*^10"”^' ”'^12,1  i+’^l 2, 1+^16 
t|-M(l  -d^^2)+l6,2+^!Q"^^^  “‘^1 2,2^'*'^i  2,2'+^l6 

tpMd-dg  jktg  3-tt;o»-tj,-M(l-dj2j)+ti2j’»t,6 
tpMd-dg  3)+t6^3+t)o-Md-di2j)+t|2j+t,6 

ti-Md-dg  3)^^  3^t,o-M(Udi2  2)+t,2,2+t,6 

tpMd-dg  3)-ht6  3+t0tt|Q-K,pMd-d,2j)tt,2,i+t,5 

t,-Md-d0  3)4.t5  3T't9i  t|o-Md-d, 2, 1^12, 1+1)5 
t,-Md-d0  3)+t0  3+t8»-t,o-Md-d,2^2^+^j2,2+'^16 

t)-M(l-d0  3)-»-t(,  3-ft8-Kd-d5  ,)+tg  ,+t,3i.ti4+t,0 

t ) -M  d -dg  3)+ 1^  3“tt8+t  ]  j+ 1 )  4+t )  0 

t2+t4+t8<-t  j  Q+t  1 ) -M  (1 -d  1 2  j )  H  i  2  j+t  1 6 

t2+t4+t8+^  1 0“'^  ^  1 2 , 1  1 2 ,  r+d  6 

t2ft44t8ftjo-H(l-cli  2^2)+^  12.2+^16 
1 2+t4+‘^8"^  ^  ’  ■‘^9 , 1  ^■+  *9 ,  1 3+^  1 4^  ^  1 6 

t2+^4-‘-^8+M3tM4 

-M(i-d,5  2)Tt5  5^2"”(’"‘^l7,l^'*‘‘l7,l*'‘l8 


<Wf 

<Wf 

<Wf 

<Wf 

<  Wf 

<  Wf 
^Wf 

^  Wf 

<wf 

'f 

< 

<'^f 

4  ''f 

<  Wf 

^  Wf 

<  Wf 

<  Wf 

<  Wf 

<  Wf 

^  Wf 
<.Wf 


<Wf 

-M(1-d,7  ,)ft,7  ,4-t,8 

-M(l-d|7 

<,Wf 

t3+ti^t8ft,(;ftipM((-d,2j)^-ti2j+ti8 

<Mf 

t3+t4+t8+t  1  o-M  (1 -d ,  2  j  )+t  j  2  j  ft  1 6 

t3ft4+.t84tjo-M(l“d  12,2^^12, 2+^6 

<«f 

f 3fttf«8"^  ( 1-d  5  j ) + 1 j  j  4.t ,  t  j  j  g 

<''f 

tjt  f4+t8^^13'^*l4 

1 3+^5"^  ^  •  “'■‘g ,  1 ^9 ,  1 3^*  1  1 6 

<*'f 

t3f tj-M  (  ) -dg^)-t  tg  ^  i  4+t  1  6 

1 3+t5-M  (I  -dg  ^  3)+  tg  ^  3^-t  1 4+t  j  6 

t3+t5-H(l-dg  2)+tg  j-Mvl-djj  2)+<^15,2 

t3ftr-H  ( 1  -dg  3)4  tg  3-M ( 1  -d ,  5  ^  15,2 

-fl(l-d|5  2^  +^15,2 

<“f 

ti  -i-  t^  ^-tjo  4-  tj| 

t]  -  Md-dg  |)  -f-  t6,l  *7  +*  ^10  -t-  ^J1 

t|  -  M(l-d6  2)  ^6,2  +  ^10  t  ^11 

<Wf 

t,  -M(l-d6^3)  +-t6^3  4-t,o4-t,, 

<“f 

t,  -  Md-dg  3)  4-  tg  3  -l-tg+  t|o-f.t|| 

h  '''■H  +  ^8  ^iO  ^11 

*3  ’‘■4  *8  ^10  +■  *1 1 

<  “f 

We  suggest  that  typically  in  the  second  formulation  of  the 
precedence  constraints  (and  it  is  also  true  of  the  first)  many  of  tfte 
constraints  -  ill  be  redundant.  This  may  be  observed  in  the  constrain 
representing  the  problem  of  this  chapter.  Th»  next  chapter  will 
develop  an  algorithm  which  will  select  and  eliminate  the  redundant 


constraints . 


Chapter  IV 


Decision  Network  Reduction 

The  previous  chapter  formulated  the  OCPH  problem  as  an  integer 
programming  problem  that  included  either  a  r-)nstraint  for  every  link  in 
a  precedence  graph  or  a  constraint  for  every  path  in  the  network.  For 
the  relatively  small  network  of  Figure  tll-l,  either  approach  gives 
approximately  forty  precedence  constraints.  In  this  chapter  will 
show  that  in  the  "path"  formulation  of  precedence  constraints,  many  ot 
the  constraints  are  dominated  and,  therefore,  'i>jt  many  of  them  may  be 
eliminated.  It  Is  also  shown  that  the  "other"  interdependency  relations 
may  make  some  precedence  constraints  infeasible  and  these,  too,  may  be 
eliminated.  Finally,  an  algorithm,  is  developed  to  reduce  decision 
networks  to  equivalent  networks  containing  only  decision  jobs  and 
maximal  distance  between  them.  This  reduced  network  could  be  used  to 
generate  a  set  of  undominated  paths  to  be  included  in  the  'nteger  pro¬ 
gramming  formulation  as  precedence  constraints. 

2.  Dominance  Tests  for  Constraint  Elimination 
Given  any  decision  network,  it  is  possible  to  order  the  jobs, 
including  the  decision  jobs,  from  to  Sf  so  that  no  job  appears 

on  the  list  until  all  of  its  predecessors  have  appeared.  We  may  now 
proceed  down  this  technologically  ordered  list  labelling  the  decision 
jobs  u(l),  u(2),  ...,  u(k),  ...,  u(h)  where  h  is  the  total 

number  of  decision  jobs  in  the  network.  Thus  we  may  uniquely  specify 
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any  decision  job  as  1=1,....,  h 

If  we  now  define  P  as  the  set  of  all  paths  from  S,  to 
Sf  In  the  network,  P  may  be  partitioned  Into  the  fuMowIng  subsets, 
P 


l^po 


r>l  p2 

»  t  ^  t  •  •  •  1 


such  that 

P®=  jo  order  paths,  those  containing  of  the  decision  jobs 
P^  =*  |ist  order  paths,  those  containing  one  of  the  decision  jobs^ 
('jth  order  paths,  those  containing  exactly  i  of  the  decision 

"  •  U 


(bth  order  paths,  those  containing  exactly 
pb  -  ) 

/jobs 


b  of  the  decision^ 


where  b  Is  the  maxlmun  number  of  decision  jobs  on  any  path  from 


to  Si 


It  is  now  possible  to  sub-partition  Into  h  subsets  in 


the  following  way. 

.1 


’■  !’■ 


p’  p’  p’ 

•  .  u(k)'  •••*'  u(h) 


(I)’  u(2) 


] 


where  that  subset  of  psths  containing  only  one  decision  job 


which  is  S 


u(k) 


Similany,  the  2nd  order  paths  may  be  partitioned  Into 

h! 


(h-2)!21 


■  j^i(l).u(j)| 


subsets  «s  follows, 

i  If  2f  1 

1*^  li  •••»  b 


and  pJ  may  be  partitioned  Into 
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h! 

(h-j;!j!  subset*. 

in  total,  the  rnaximum  number  of  subsets,  ec'h  containing  a  unique 

combination  of  decision  Jobs  will  be 

hi  •  .  h!  .  .u 

‘  ^  {h-2)!2!  ^  (h-j)!j! 

The  actual  number  of  subsets  required  would  be  much  less 

than  this  however.  The  set  of  interdependency  relations  from  Chapter 

ill  k(i) 

2  iij  “  I 

J=i 

guarantees  that  many  combinations  of  decision  nodes  are  not  feasible. 

Now  . .  u(k) ,  ...,  u(n) 

is  a  set  of  paths  containing  a  particular  combinatior  o>'  j  of  the  h 

decision  Jobs.  For  any  OCPH  problem,  any  solution,  that  is  any  selec> 

tion  of  decision  jobs  to  perform  that  meets  the  Interdependency 

constraints,  will  have  one  of  the  following  mutually  exclusive  prop* 

ert.es.  Either  all  of  the  j  decision  Jobs  in  P-i/  >  /,  * 

■'  ■*  U\m; . u(k)  ,  . . . 

u(n)  arc  perfonned  or  at  least  one  of  them  is  not  performed. 

If  one  of  them  is  not  performed,  then  the  precedence 
constraints  representing  all  of  the  paths  in  this  set  will  not  be  bind¬ 
ing  on  the  early  start  of  the  Finish  node. 

That  is  if 

‘*u(m)-^  •••’“^‘*u(k)-^'  •^‘^u(n)  <  ^ 

then  -.he  ccnstraint  representing  any  path  in  the  subset  will  meet  the 


condi tion 


5G 


+  “uW  -  ••••  -"<'-''u(k)>+  ‘u(k)  -  •••• 

+  tu(„)  +  T  <0  <W, 

where  T  represents  the  sum  of  the  times  of  non-decision  Job»  on  the 
path. 

If  ali  the  Jobs  in  the  subset  are  performed,  then  the  finish 
day  of  the  project,  wi I i  be  constrained  by  equations  representing 

all  paths  in  the  subset.  Since  the  constraints  are  identical  except 
that  the  value  for  T,  the  sun  of*  the  times  of  the  non-decision  jobs 
on  a  particular  path,  vary,  it  is  only  necessary  to  retain  the  con¬ 
straint  that  has  the  maximum  value  for  T.  All  constraints  representing 
other  paths  in  the  subset  will  be  dominated  by  the  one  chosen.  Thus  we 
require  only  one  path  for  each  possible  combination  of  v'rrision  nodes. 

We  wili  term  this  path  an  "undominated"  path. 


3.  Implementation  of  Dominance  Tests  for  Path  Elimination 

This  section  will  develop  an  algorithm  based  on  the  repetikive 
application  of  the  largest  path  calculation  to  generate  an  undominated 
set  of  paths  in  a  decision  network.  As  stated  above,  t.Se  undominated 
path  will  simply  be  the  longest  path  in  a  subset  of  paths  which  ali 
contain  a  particular  combination  of  decision  icbs.  For  example,  in 

the  subset  ^y(j)  .  j.  'fitis  not  empty,  we  must 

calculate  the  longest  path  through  jobs 


*u(i)'  u(2) 


•  (k)' 


'u(j)’ 


which  goes  through  no  other  decision  Job.  This  is  equivalent  to 
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findino  i.he  longest  partial  path  connecting 
^5  Sy(,).  and  Sy(2)‘ 

and  combining  them  to  form  the  complete  path.  Again  each  partial  path 
must  not  contain  other  decision  jobs  since  this  would  result  In  the 
calculation  of  the  longest  peth  through  a  different  combination  of 
dec! Sion  Jobs. 

A  routine  based  on  the  usual  longest  path  calculation  (CPM; 
was  written  to  generate  the  maximal  distances  between  all  decision 
wdes  and  S^,  Sf.  The  algoritnm  uses  two  time  values  for  each  job 
in  the  network.  The  "actual  time"  is  the  estimated  completion  time 
for  the  job  which  remains  constant  throughout  while  the  "current  time" 
is  oqual  cither  to  actual  time  or  -M,  a  large  iiegative  number.  As  we 
shall  see,  if  a  job  time  is  taken  as  -M  ,  then  no  path  through  the  job 
has  a  non-negative  length  and  ali  such  paths  will  be  ignored.  A  flow 
chart  for  the  algorithm  will  now  be  presented  and  reference  will  be 
made  to  Tables  IV-I,  2  and  3  which  show  some  steps  of  the  algorithm 
for  the  decision  network  of  Figure  1 1 1-1. 

1.  List  the  jobs  of  the  decision  network  in  technological  order 

2.  Form  a  matrix  with  a  row  for  each  job  in  the  network.  Colunun  < 
will  contain  the  job  number;  Column  2,  the  actual  time;  Column  3, 
the  job's  predecessors,  Colunns  4,  5  •nd  S,  the  current  time,  the 
early  start  time  and  the  early  finish  time,  respectively,  (ic. 
Table  IV.«) 

3.  Identify  J«cision  Jobs  in  Column  I.  (in  Tables  IV-1  brackets 
are  used) 
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4.  Set  the  current  time.  Column  4,  for  all  decision  jobs  at  -H,  for 
all  other  Jobs  set  the  current  time  equal  to  actual  time.  Set 
counter  SIsr  S^.  Go  to  Step  8. 

5.  Beginning  at  II,  search  down  Column  I  for  the  first  decision  job 

below  II,  say  this  is  Set  l|=.Sy^^j. 

6.  Set  the  early  start.  Column  5,  for  all  jobs  above  II  at  -H. 

7.  Set  the  current  time  of  II  and  the  early  start  of  II  at  0. 

8.  Calculate  the  early  start  time  of  all  jobs  beneath  II  on  the 
technologically  ordered  job  list.  If  the  early  start  value  cal¬ 
culated  Is  negative,  enter  -H  in  the  early  start  colunn.  To 
calculate  the  early  start  of  the  job  directiy  beneath  II,  list  all 
the  predecessors  of  the  job  and  for  each  predecessor  add  the 
predecessors  early  start  time  to  its  current  time.  The  largest 

of  these  values,  for  all  predecessors,  is  the  early  start  time  for 
the  job  being  considered.  Then  proceed  to  tlie  next  job  on  the 
technologically  ordered  list  for  the  next  early  start  calculation. 

9.  For  all  decision  jobs,  add  early  start  time  to  actual  time  to 
compute  early  finish  time.  Column  6. 

10.  For  all  decision  jobs  beneath  II  on  the  technological  list  and  for 
S^,  record  the  early  finish  time  if  It  is  non-negative.  These 
times  will  be  the  longest  sub-path  from  the  finish  of  decision  job 
II  to  the  finish  of  the  decision  Job  being  examined  (the  longest 
path  containing  no  other  decisio.i  Job). 

11.  if  II  Is  the  last  decision  job  in  the  technological  list,  go  to 
Step  12.  If  not,  to  Step  3. 
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12.  HALT. 

This  algorithm  has  been  applied  to  the  problem  of  Figure 
MI-1  with  results  shown  in  Table  IV-3<  Each  non-blank  entry  in  this 
matrix  indicates  the  length  of  the  maximal  path  between  pairs  of  deci¬ 
sion  jobs  dr  between  the  Start  and  Finish  nodes.  Note  that  the  entries 
for  line  one,  the  maximal  distances  from  to  all  decision  Jobs,  are 
taken  di recti/  from  Column  6,  Table  IV-1.  Similarily,  the  entries  for 
row  two,  the  distances  from  Sg  ^  to  all  decision  jobs,  are  taken  from 
Column  6  of  Table  IV-2.  Blank  entries  in  the  matrix  indicate  that  no 
path  not  containing  a  decision  job  connects  the  to  jobs.  The  output  of 
a  program  which  generates  data  for  the  matrix  of  IV-3,  and  estimates  of 
computation  times  for  van..js  decision  networks  are  given  in  Appendix  B. 

4.  Feasibility  Tests  for  Path  Elimination 
Consider  a  particular  subset  of  paths 

~  ^  •  •  •  >  j 

For  the  constraints  representing  these  paths  to  constrain  w^,  the 
following  condition  must  hold 

■‘.(0=  J 

This  condition  may  be  contradicted  by  the  "other"  interdependency 
constraints 

ie.  d/.,  \  -+-d,..—  1 

u(l)  u(2) 

If  the  relation  5'  *.  j  can  be  shown  to  be  infeasible,  then 

i  =  I 

all  paths  in  Fy(j)  i  —  I,  2,  ....  j,  that  is  all  paths  containing 
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1 

2 

3 

4 

5 

6 

Job 

Actual 

Time 

Predecessors 

Current 

Time 

Early 

Start 

Early 

Finish 

s 

0 

0 

0 

0 

0 

1 

12 

Ss 

12 

0 

12 

2 

10 

Ss 

10 

0 

10 

3 

8 

Ss 

8 

0 

8 

4 

4 

2.3 

4 

10 

14 

5 

4 

3 

4 

8 

12 

(6,1) 

4 

1 

~M 

12 

16 

(6.2) 

6 

1 

-M 

12 

18 

(6,3) 

8 

1 

-M 

12 

20 

7 

2 

l.(6.1) 

2 

12 

14 

8 

3 

4. (6.3) 

3 

14 

12 

(9.0 

5 

CO 

“M 

17 

22 

(9,2) 

10 

5 

-M 

12 

22 

(9,3) 

15 

5 

-H 

12 

27 

10 

10 

(6, 3). (6. 3) ,7.8 

1C 

17 

27 

11 

5 

10 

5 

27 

32 

(12,1) 

4 

10,11 

-H 

32 

36 

(12,2) 

5 

10 

-K 

27 

32 

13 

2 

8.(?,0 

2 

17 

19 

14 

8 

(9, 2). (9, 3), 13 

8 

19 

27 

(15.1) 

9 

93, S, 

-H 

0 

9 

(15,2) 

13 

9’.S^ 

-M 

0 

13 

16 

10 

(12. 0.(12. 2), ’i*. (15.1) 

10 

27 

37 

(17.1) 

1 1 

(15. 2), S, 

-M 

0 

II 

(17,2) 

3 

(15.2)55 

-M 

0 

3 

18 

5 

(17, 0.(17. 2). Sj 

5 

0 

5 

Sf 

0 

16, 18. (15, 2). 11 

0 

■*7 

il 

Table  iV-l 
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1 

2 

3 

4 

5 

6 

Job 

Actual 

Time 

Predecessors 

Current 

Time 

Early 

Starr 

Early 

Finish 

5$ 

0 

0 

0 

-M 

1 

1-2 

Ss 

12 

-M 

— 

2 

10 

Ss 

10 

■  rl 

— 

3 

8 

S 

s 

8 

-M 

— 

4 

4 

2.3 

4 

-M 

— 

5 

4 

3 

4 

-M 

— 

(6.1) 

4 

1 

0 

0 

— 

(6.2) 

6 

1 

-M 

-M 

-M 

(6.3) 

8 

1 

-M 

-M 

-M 

7 

2 

l.(6.l) 

2 

0 

8 

3 

4,63 

3 

-M 

-M 

(9,1) 

15 

5.8 

-M 

-M 

-M 

(9.2) 

10 

5 

-M 

-M 

-M 

(9,3) 

15 

5 

•'i 

-h 

-M 

10 

10 

(6, 2). (6, 3)  ,7, 8 

10 

2 

12 

n 

5 

10 

5 

12 

17 

(12.1) 

4 

10.11 

-M 

17 

21 

(12.2) 

5 

10 

-M 

12 

17 

13 

2 

8, (9,1) 

2 

-M 

-M 

14 

8 

(9. 2), (9, 3). 13 

8 

-M 

-M 

(15.1) 

9 

(9.3).Ss 

-M 

-M 

-M 

(15,2) 

13 

(9.3) .Sj 

-M 

-M 

-M 

16 

10 

(12.I).?12.2).14,{15.1) 

10 

-M 

-M 

(17,11 

11 

(•5.2).Sj 

-M 

-M 

-M 

(17.2) 

3 

(15.2)  .$5 

-M 

~'A 

-M 

18 

5 

(17,1).(I7,2).S^ 

5 

-M 

-M 

Sf 

0 

16. 13, (15. 2). 11 

0 

17 

17 

Table  IV-2 
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Reduced  Network  Hatrix 


T  0 


6,1  6,2  6,3  9,1  9,2  9.3 

12,1 

12,2 

15,1 

15,2 

17.1 

17.2  Sf 

16  18  20  22  22  27 

36 

32 

9 

13 

It 

3  37 

21 

17 

17 

19 

15 

15 

8 

22 

18 

23 

20 

18 

9  13  18 


10 

10 

10 

II  3  0 

5 
5 


Table  IV-3 
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the  combination  of  decision  jobs. 

^ii(i)  '  =  U  2,  j 

may  be  eliminated  from  the  integer  programming  formulation  of  the  OCPH 
problem.  Paths  not  eliminated  by  the  interdependency  constraints  will 
be  termed  '‘‘feasible"  paths. 

5.  Application  of  Pori  nance  and  Feasibility  Tests  to  a  DCPM  Network 

When  the  tests  of  the  previous  three  sections  are  applied 
to  precedence  equations  of  the  example  Figure  lil>l,  the  reduced  set  of 
equations  of  Table  Wl-k  result,  for  example,  the  set  |)  contains 

the  j  paths 

t]  +  *7  -l-tiQ  -M(|-d,2j)  f  t,2j  +  t|6  4  Wf 

^1  +  ^7  +  ^10  *11  *12,1  +  *16 
*2  t  *4  “*"*8  ■^*10  '*"*11  ■^^’■**12,1^  *12,1  +  *16  4'^f 

*2  +  *4  +  *8  •+•  *10  +■  *12,1  1*  *16  4^f 

*3  -h  *^4  +  *8  ^  *10  t  *11  "”^’"^12,1^  *12,1  1'*16  ^  S 

*3  *4  1-*8  +  *I0  "”^’’**12,I^  *"’'12,1  -t-‘l6  ^'^f 

These  may  be  rewritten 


-M(l-d 

12,1 

)  +  38 

< 

^f 

-M(l-d 

12.1 

)  +  43 

-M(l-d 

12,1 

)  +46 

< 

^f 

-M(l-d 

12.1 

)  41 

'^f 

-KKd. 

1 

2.1^ 

+  44 
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-H(l-dj  3) 

37 

4  43 

<.Mf 
<  Mf 

-H(|. 

4  42 

-Md-d, .,) 

4  40 

-H(I-d,_3) 

4  45 

<Wf 

.M(l-d|2_,) 

-4  46 

-H(l-d,^_^) 

4.  42 

-M(l-dj  ,)  -Md-dij  j) 

■4  43 

-Md-d^_2>  -"('-''la.i' 

+-37 

-Hd-dj  3)  -Md-d,j_j,) 

-4  kb 

-Hd-dj_3)  -Md-d|j ,) 

T  46 

-"d-d,.,;  -"('-<,5,3)  2) 

+  56 

-Hd-d.5  2)  -Kd-d|3 1) 

4  29 

<Wf 

.Hd.d,5_2)  -Md-d,3  2) 

21 

<Wf 

-M(:-d|2,|) 

4  16 

<«f 

'”*'■'‘17,2* 

4  8 

-«d-d3_3)  -Hd-d,5  2)  -Hd-iij.,) 

-+  52 

^  Wf 

Tabl«  iV.4 
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-M(l-d,2  ,)  t39 

With  -M(l-dj2  j)+  ^6^  y^  in  the  prcbietn,  the  other  5  constraints 
are  redundant. 

We  will  now  give  an  example  of  constraint  elimination  by 
feasibl 1 1 ty  testa. 

The  equation 

-M(l-d^^j)  -  M(I-dg  j) -4- 48  Vp  will  hold  only  If 

'*6,3  '^‘‘9,1  ~  ^ 

However,  the  Interdependency  constraints 

Sj  <  ‘^6.2 

and 

<*6,2'^  * 

violate  this  condition.  Therefore,  the  equation  may  be  dropped. 

6.  Lower  Bound  Calculation 

The  set  of  path  constraints  has  now  been  reduced  from  one  per 
path  in  the  original  decision  graph  to  one  par  combination  of  decision 
Jobs  lying  in  a  path  from  to  Sp.  The  number  may  be  even  further 

reduced  because  of  feasibility  tests  on  the  interdependency  relations. 

We  will  now  define  set 
M  =  ^all  feasible,  undominatcd  paths 
and  set 

M|  *  the  subset  of  paths  in  M  that  contain  any  of  the  decision 

jobs  !n  decision  set  S.  ,  a  multi-job  set'^  Each  path  in  Mj  contains 


exactly  one  of  S,j,  J  =  I ,  2,  . . .  k(l )  since  we  are  assuming  that 
only  one  job  from  each  multi -job  set  may  be  aelected..  but  they  may  also 
contain  other  Jecislon  Jobs. 

We  will  now  define  a  "feasible  solution"  to  a  DCPM  problem  as 
a  selection  of  a  set  of  decision  jobs  that  will  satisfy  all  "alternative" 
and  "other"  interdependency  constraints  on  the  problem.  This  set  of 
decision  jobs  will,  with  the  non-decision  jobs,  form  a  usual  project 
graph  of  the  critical  path  method.  Given  the  project  graph  resulting 
from  a  "feasible  solution",  we  state  that  one  of  the  paths  from  Mj 
will  be  included  In  that  graph.  If  no  such  path  exists,  then  no  path 
through  S|  exists  in  the  network.  This  would  violate  our  assumption 
that  Sj  is  a  predecessor  and  Is  a  successor  to  all  jobs  in  the 
project  graph. 

Each  constraint  representing  a  path  from  M|  wi 1 1  be  of  the 
general  form 

-H(l-dij)  .  -•'(1-0+  “f 

where  Tj  ^  i s  the  length  of  path  e  in  Mj. 

If  one  of  the  paths  from  Mj  must  be  performed,  then  a 
lower  bound  on  Wf  is  the  lowest  value  of  Tj for  all  paths  In  M, . 
Furthermore  if  we  calculate  a  lower  bound  for  every  decision  node,  we 
may  select  •^he  highest  ot  the  iower  bounds  as  a  lower  bound  on  W^. 

Any  undominated,  feasible  pati’  which  is  shorter  than  this  r<ound  may  be 
eliminated  since  it  is  redundant,  given  the  bound  on  W^. 
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An  Appiicatton  of  Dower  Bound  Calculation 

The  reduced  set  of  equations  from  Table  IV-4  will  be  used  to 
illustrate  the  dominance  tests  described  above.  The  only  feasible  paths 
through  Sg,  that  is  the  set  Mg,  are  the  following  paths, 


-HU-d,  J  •+•  43 

0|3 

4 

-Md-d^  ,)  .H(l.d,2,5)H-47 

4  w, 

-M(l-4^  2^  -  H(l-d,2  2)  +^3 

-M(l-d^^^)  -  -+-52 

The  minimum  of  the  T. 

1  ,C 

.,4,  is  43  and  therefore  4’  is 

a  lower  bound  on  the  length  of  the  critical  path.  Similarly,  the  lower 
bound  provided  by  is  40,  Mj2  Is  42,  i s  21  and  Mj^  is  8. 

These  calculations  are  shown  in  Table  IV-5.  The  highest  of  the  bounds 
is  then  43  and  all  paths  shorter  than  or  equal  to  this  may  be  removed. 
For  this  particular  problem  this  reduction  technique  does  not  give  any 
further  improvement.  The  final  set  of  predecessor  constraints  may  "ow 
be  wri tten 


(1) 

-H 

43 

4 

«f 

(2) 

-"<'-“9.3' 

-f- 

45 

4 

“f 

(V 

-"<'-“12. i> 

+ 

46 

4 

"f 

(• 

46 

4 

(5) 

k «  /  •  •  \ 

-nvi-u^  1/ 

-r 

47 

4 
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(6)  -M(r-d  ,)-M(l-d..  ,)  -H(l-d,_  J  4  W. 

9.3  i5.2  17,2  f 

(7)  '♦“52  4. 

(8)  3)-H(l-d|^  j)  -M(i-d|y  j)  +56  C 

By  selecting  appropriate  values  for  H  and  combining,  these 
may  be  rewri tten 


(1) 

(3) 

3<‘|2.1 

H-  43 

(2) 

(4) 

4  39 

4 

(5) 

(7) 

4d^  j  4  9d^  j4  4d,2  j 

4  39 

4  Wf 

(0  (8)  >5-5.3 +  '’^5,a+"V,-^’-, 7.2 + 

The  final  problem,  now  In  reduced  path  form,  may  be  solved  by 
any  standard  integer  programming  technique. 


Chapter  V 


A  Network  Algorithm  for  Restricted  Types  of  Integer  Piograms 


In  the  literature  review  of  Chapter  1 1 ,  we  referred  briefly 
to  several  integer  progranvning  routines  for  problems  containing  0,1 
variabUs.  In  our  terms,  these  techniques  could  optimize  planning 
problems  wi th 'job  alternative"  and  "other"  types  of  interdependency. 

We  will  now  develop  an  algorithm  based  on  the  calculations  of  the 
critical  path  method  which  will  solve  certain  restricted  types  of 
integer  programs,  A  numerical  example  of  the  algorithm  will  be  given 
and  then  the  algorithm  will  be  applied  to  a  problem  taken  from  the 
decision  jobs  and  interdependency  constraints  of  a  decision  network. 

2.  The  Network  A1 qori thm 
Consider  the  integer  programming  problem 
(I)  Maximize  wb +  bQ 

Subject  to  wA^  ^  c 

w  o  ,  w  integer 

where  b  i  s  ar;  m  x  I  column  ve«.tor,  b^  is  a  scalar,  c  is  a  1  x  n 
row  vector  and  o  and  w  are  I  x  m  row  vectors  with  values  assigned 
to  w  so  as  to  maximize  che  objective  function  subjec'.  to  the  constraints 
of  (l).  A°  i s  an  ir  x  r.  matrix  which  meets  the  following  assumption. 
Assumption  i:  All  rows  of  A®  contain  a  mcximiw  of  two  non-zero 
entries  and  these  are  found  in  adjacent  columns.  The  programming 
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problem  may  now  be  interpreted  as  a  network  problem  in  the  following 


way.  For  each  variable  Wj 


=  1,2,  ...,  m  define  a  set  o^ 


0-1  variables,  one  for  each  possible  integer  level  of  w- .  These 
would  be 

W:  =  7  w.'  ,  w.  , ,  . . . ,  w.  . ,  . . .  w  ,  .  I 

'  '.o’  >.r  i,j  i.k(i)J 

where  w,  .  =  X  '  variable  wj  =  j 
’  0  otherwise 

and  k(i)  is  the  maximum  possible  value  of  w: 


Tne re fore 


k(i) 

S  W.j  =  1 

j  0 


and  w.  =  2  >,j 


A  feasible  solution  to  the  programming  problem  will  then 
consist  of  a  choice  of  exactly  one  variable  from  each  of  the  following 
sets  such  that  the  constraints  of  (l)  are  met. 


{“i.o’  "i.r 
"  “2,1’ 


•  “l.k(l)/ 
'  “2,l<(2)_}- 


{w  ,  w 

m,o’  tr., 


w  ,  ^ 

'  m,k(m)  ) 


L i nk  each 


Now  let  the  Wjj  be  nodes  in  a  network  formed  as  follows, 
w,  I  ,  j=  1,2 . k(l)  to  each  w,  .  ,  j  =  1  ,2,  - k(2) 


by  a  directed  line  segment.  Then  link  each 

d  j=l.2 . k(t)  to  each 

1  >J 
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j,  j—  l(2(  k( t'f"  1 ) 

whe  re  t  =  2 ,  3 .  m- 1 

Since  every  variable  in  each  set  is  linked  to  all  variables 
in  adjacent  sets,  the  resulting  network  contains  a  path  for  all  combine-’ 
tions  of  Integer  values  of  Wj ,  . .  w^.  Since  there  are  no 

links  within  a  variable  set,  the  restriction  k(i) 

.2  • 

J=0 

is  maintained.  Thus  every  feasible  and  non-feasible  solution  is  repre¬ 
sented  by  a  path  in  the  network  and  every  path  is  a  feasible  or  non- 
feasibie  solution.  Here  a  feasible  solution  is  defined  as  the  selection 
of  a  level  for  each  of  the  w.  variables  that  meets  the  constraints 
of  (1).  The  network  may  now  be  modified  to  eliminate  paths  represent¬ 
ing  non-feasible  solutions. 

First  vje  will  show  how  to  calculate  k(!),  the  maximum  integer 
value  of  variable  wj .  Assume  ail  variables  Wj ,  i  =•  2,  3,  .  m 

have  the  value  o,  that  is  Wj  ^  =  i . 

Then  from  the  constraints  in  wA°  c  of  (1)  that  have 
only  positive  entries,  it  is  possible  to  calculate  a  maximijm  value  for 
W|  ,  that  is  k(l).  Similarly,  we  may  calculate  maximin  values  for 
al I  other  var i ables . 

We  now  have  the  graph 
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Where  at  each  level  wt;  have  nodes  representing  all  the 
feasible  integer  values  of  one  variable,  and  each  of  these  nodes  is 
connected  to  all  nodes  of  the  suceedirg  level  (or  variable). 

The  constraints  of  (l),  wA°  c  may  now  be  introduced 
exactly  into  the  network.  From  our  earlier  discussion,  it  is  known 
that  each  path  in  the  network  represents  a  particular  combination  of 
integer  values  of  the  variables,  Wjj  and  each  path  is  therefore  a 
potential  solution  to  the  problem.  For  example,  one  path  might  be 


If  such  a  solution  (path)  violates  one  of  the  constraints  of  (1),  we 
remove  the  path  by  eliminating  particular  d'rected  line  segments  that 
connect  the  integer  levels  of  the  two  variables  contained  in  the  con¬ 
straint.  Assumption  i  guarantees  that  only  two  variables  will 
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involved  in  the  constraint  and  that  they  v^i  1 1  be  adjacent  in  the  yraph. 
Graphical  interpretations  of  certain  constraints  follow. 

(i)  '^i  4.  1  “  ® 

Mil  links  between  variables  W|  and  will  be  removed 

with  the  exception  of  those  now  listed. 


Where  the  sum  of  the  values  represented  by  the  connected 
variables  equals  "a" 

(iij  wj  =r  w^, 

all  links  between  variables  Wj  and  will  be  removed  with  the 

exception  of  those  now  listed. 
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where  K(i )  ^  k(i  +  l) 

When  all  the  constraints  of  (l)  have  been  introduced  into  the  network 
(by  link  elimination),  we  will  term  the  resulting  graph  a  "feasible” 
network. 

Theorem  1:  Every  path  in  the  feasible  network  is  a  feasible  solution 
to  the  problem  and  every  feasible  solution  may  be  represented  by  u  path 
in  that  network. 

PROOF:  Between  each  level  of  the  graph,  non-feasible  links 

are  removed.  If  the  constraining  equations  have  been  interpreted 

correctly,  remaining  links  connect  feasible  pairs  of  variables. 

Now  assume  a  path  exists  in  the  reduced  netv/ork  from  the  set 

of  variables  W]  to  the  set  w^^^  and  variables  W|  . 

c  path.  Variable  set  Wj  may  only  be  combined  in 

constraint  equations  with  variable  set  if  consideri  ng  al  1  the 

constraints  on  w,  ,  w.  the  combination  of  w,  ,  w.,  .  is  feasible 

12  I  , a’  2,0 

then  Wj  ^  may  appear  in  any  solution  that  contains  W2  {j.  The 
argument  may  be  extended  to  w, ,  wk  and  finally  to  w  . 

Therefore  any  path  in  the  reduced  network  is  a  feasible 

solution. 

Now  assume  a  feasible  solution  exists  for  which  a  path  does 
not  exist.  This  implies  that  bctvx:en  two  adjacent  variable  sets  Wj  , 
w.  I  there  is  no  connecting  link  in  the  network  berv-een  the  variables 
from  these  sets  which  arc  in  the  fcas.bi.!  solution,  that  is  between 


76 


”i  a’  '^itl  b’  links  arc  only  removed  if  in  combination  they  arc 

infeasible  under  the  constraints,  then  a  solution  containing  w.  and 
Wi  1  b  is  not  feasible.  Therefore  every  feasible  solution  must  be 
represented  by  a  path. 

Theorem  2:  If  each  node,  Wjj ,  of  the  network  graph  is  evaluated  at 
jbj ,  the  solution  to  maximization  problem  is  the  longest  path  in  the 
network. 


PROOF:  By  Theorem  1  every  path  is  a  feasible  solution  to  the 
problem.  Since  Wjj  represent  the  original  variable  w;  at  i  iteger 
level  j,  the  value  of  w.j  in  the  functional  will  be  jbj  where  b| 
is  the  contribution  of  a  unit  of  w.  to  the  functional.  The  path 
length  for  any  solution  is  simply  the  value  of  the  functional.  The 
variables  on  the  longest  path  s.c  those  that  give  an  optimum  solution 
to  the  maximization  problem.  Similarly  the  solution  to  a  minimization 
problem  r  iy  be  found  by  determining  the  shortest  path. 

The  familiar  rules  of  Critical  Path  Scheduling  will  calculate 
the  longest  path  in  network.  In  this  context  the  early  start  (E.S.) 


for  a  variable  wi ;  is  interpreted  to  be  the  maximum  value  of  the 


criterion  function  for  variables  Wj  to  w._j  that  may  be  feasibly 
combined  with  w.j  (maximization),  late  start  (L.S.)  is  the 
maximum  feasible  value  of  the  criterion  function  less  the  maximum  value 


of  the  criterion  function  for  variables  Wjj  to  that  may  be 

feasibly  combined  with  '^jj- 

In  addition  the  “sla<.k"  values  calculated  by  that  technique 
have  meaning  in  this  context.  The  slock  values  of  a  job  indicates  the 
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differuice  in  length  between  the  longest  path  on  v/hich  the  job  is  found 
and  the  longest  path  of  the  network.  Since  in  the  integer  programming 
problem,  the  path  length  is  equivalent  to  the  value  of  the  functional, 
and  a  solution  is  any  path  in  the  reduced  network,  the  slack  value 
indicates  the  minimum  change  in  the  functional  which  will  be  realized 
if  Wjj  is  brought  into  the  solution.  Thus  the  "slack*'  values  are 
exact  evaluations  for  every  Wjj,  that  is,  for  every  variable  w.  at 
every  feasible  level.  These  evaluations  may  be  calculated  for  either 
maximization  or  minimization  problems.  As  we  would  expect,  evaluators 
for  variables  in  the  solution  vn  1 1  be  0. 

Theorem  3:  If  the  variables  in  a  problem  which  meets  Assumption  1  may 

be  listed  so  that  no  constraint  exists  betv/een  Wj  ,  w.^j  i  si . . 

m-1  then  the  optimum  solution  to  a  problem  made  up  of  the  variables 
Wj  .  .  .  Wj  and  the  co;istraint5  relating  to  them  and  the  optimum  of 
problem  containing  the  variables  Wj^j  .  .  .  w^  wi*l  together  give  an 
optimum  solution  to  the  combined  problem  and  the  problems  are  independent. 

PROOF:  Let  the  Of)timum  solution  to  the  first  problem  be  a 
path  wj  ^2  b  •  •  •  '^i  r  ^ solution  to  the  second  be 

w.  ,  ^  .  w  ,  .  Si  nee  no  consti  ai  nt  exi  sts  bet\'«cn  w.  .  w.  ,  ,  .  all 

I  I  .3  m,b  I  '  I-*- 1  • 

nodes  arc  connected.  Thercrorc  the  longest  path  in  the  complete  net¬ 
work  will  consist  of  the  longest  patii  ab'we  the  level  w.^^l  plus  the 
longest  path  bciow  Wj  .  This  will  bf  the  sum  of  the  individual  solu¬ 
tions  to  the  two  problems. 

For  purposes  of  the  algorithm  it  is  possible  to  have  any  set 
of  mutually  exclusive  variables,  rather  than  a  scries  of  mutu-illv 
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exclusive  levels  of  the  same  variable  in  each  variable  set  of  the  graph. 
Assumption  1  must  hold  for  the  set  of  variables  in  their  relation  to 
other  sets  of  variables  however.  We  also  note  that  the  contribution 
of  the  variable  to  tha  functional  need  not  be  a  linear  function. 


3.  A  Numerical  Example 
Max.  6.tj  7x2  8xj  ~f*2x2^ 

St.  Xj  4 

3Xj  4-^X2  >0 

><2  +  ^4  ~  ' 

xi  .j.  X3  =  6 

X.  ^  0,  'nteger 

The  maximum  values  for  the  variables,  determined  from  the  constraints, 
are 

X,  =  3.  1.  x^  ~  6,  x^=  1 


Figure  V-I  illustrates  the  beginning  graph  and  Figure  V-2  the  graph 
with  links  removed  to  satisfy  the  constraints.  In  these  networks,  the 
variables  when  listed  In  the  order  xj.  Xj ,  x^,  and  x^,  satisfy 

Assumption!.  The  optimum  solution  consists  of  those  variables  for 
which  E.S.=  L.S.  It  is  possible  that  several  such  solutions  would 


exist.  In  Figure  V-2  these  variables  are  seen  to  be  w,  , 

I  kO 


w,  ,  and  w,  that  is  w  —  o, 
3,6  4,0  I 


w.,  =  I  ,  Wj  =»  6  and  w^^  =0. 


The 


value  of  the  functional  for  this  solution  is  (6) (o)  (7) ( I )  *4  (8) (6) 
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•+  (2){o)—  55*  In  addition,  we  see  that  variable  at  levels  o,  I 
and  2  is  infeasible  since  variables  wj  j ,  2  I '*  on  no 

complete  path.  Finally,  we  can  calculate  that  the  best  solution  con¬ 
taining  variable  X2  at  level  o  (v/2  q)  would  have  a  functional 
value,  (^-53)’=  “5»  five  units  lower  than  the  optimal  solution  of 
55. 

4.  Violations  of  Assumption  I 


)f  Assumption  1  is  violated,  it  may  be  possible  to  adapt  the 
algorithm  to  the  resulting  problem.  If  to  the  problem  of  Section  3» 
we  add  the  constraint  x/^  <  then  Assumption  1  is  violated.  Since 

X|  is  involved  in  constraints  with  the  three  other  variables.  It  is 
not  possible  to  list  the  variables  so  that  all  constraints  involve 
adjacent  variables.  The  problem  may  still  be  solved  in  several  ways. 

(i)  Define  a  new  variable  xj^  representing  all  possible 

combinations  of  variables  x,  and  x...  These  are  x^»°,  x^»\  x?'?, 

^  4,0  4,0  4,1 

Constraint  X2  x^  1  rules  out  Tbe  resulting  network  and 

its  solution  is  shown  in  Figure  The  functional  is  now  valued  at 

53  wi tb 

X,  =  I,  ^2  *  '•  "^3  =  5.  x^  =  0 

(ii)  The  technigue  of  truncated  enumeration  may  be  used  on 
thet  part  of  the  problem  which  cannot  be  introduced  into  the  graphical 
network.  A  straightforward  branching  strategy  would  be  t.  '  select 

in  turn  the  variables  that  (a)  appear  In  constraints  and  (b) 


Figure  V-3 
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cannot  be  listed  in  adjacent  levels  In  our  graph.  At  each  level  we 
would  develop  the  node  with  the  maximum  evaluation  (Maximization).  At 
each  node  the  bounds  would  be  calculated  using  our  pfoblem  network  with 
the  network  modified  to  represent  conditions  on  the  path  rf  the  enumera¬ 
tion  tree,  above'  the  node  in  question.  The  problem  of  this  section 
will  now  be  solved  with  this  procedure  (that  i s  x^  X|  is  added  to 
the  problem  of  Figure  IV-2). 


Therefore  the  selection  of  <•’»  ^  gives  an  optimal 

sol uticn  of  53. 

5.  Application  to  a  OCPM  Problem 

The  method  of  this  chapter  has  been  developed  to  handln 
problems  arising  in  Decision  CPM  networks.  A  project  having  six 
decision  sets  might  give  the  following  problem. 

ox|  -h  lOx^  •+■  20Xj  f  -h  7x5  9x^ 

0x7  +  xg  f  OX|g  4-  ]  r  I5X|2 


Minimi /e 
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ox, 3 +6x14  +  ox,5  +  I3x,5  +  k 
Subject  to  "job  alternative"  interdependency 


Job 

1 

X,  X2  X3  =• 

Job 

2 

X4  4-  X5 

Job 

3 

4-  X8  +  “ 

Job 

4 

’‘10  T  ^11  -f  ’‘U 

Job 

5 

x,3  4  x,4  =  I 

Job 

6 

^15  +  ^16=  ’ 

and  "other"  Interdependency. 

Xj*  Xg 

’‘lO 

XjSSXg 

X2'b  xj  4  I 

Xg 


If  the  job  sets  are  ordered  6,  I,  2,  3t  5,  then  ail  constraints  may 
be  Included  In  the  first  network.  Figure  V-4  illustrates  the  reduced 
graph  and  the  evaluations  of  the  variables.  The  minimum  cost  solution 
to  thi s  problem  is  Xj ^  =s  1  *  X2  I ,  x^  =  1  ,  x^  —  I  ,  Xj^  =r  1  and 
X|2  »  I  for  a  total  cost  of  10. 


FIGURE  V-4 


Chapter  VI 


Branch  and  Bound  Algorl thms  for  the  DCPH  Problem 

In  Chapter  IM  the  OCPM  problem  was  formulated  as  an  Integer 
programming  problem  but  for  large  problems  this  approach  would  require 
substantial  computation  time.  Forthis  reason,  it  was  decided  to 
examine  the  applicability  of  restricted  enumeration  solution  techniques. 
These  methods,  under  the  geneial  name  "branch  and  bound"  have  been 
used  for  solving  a  series  of  lather  difficult  cominatoria!  problems. 

The  applications  include  the  travelling  salesman  problem  L27,  Sij  , 
a  version  of  the  plant  location  problem  [78l,  truck  routing  [sz] 
and  production  sequencing  problems  [42,  47,  6l].  This  work  has 
demonstrated  that  for  a  wide  range  of  problems  these  techniques  are 
efficient. 

The  work  "branch"  in  this  context  refers  to  a  specific 
decision  rule  for  enumerating  a  tree  of  all  possible  solutions  to  the 
problem.  The  "bound"  term  implies  that  at  each  node  of  the  tree  a 
maximum  value  (for  maximization  problems)  oi  the  criterion  function  is 
calculated.  This  bound  on  the  value  of  the  criterion  function  for  a 
partial  solution  may  indicate  that  no  optimal  solution  will  be  found 
that  contains  the  particular  partial  solution  and  therefore  search  on 
the  path  may  be  terminated.  Alternately,  the  calculation  may  show  that 
the  partial  solution  could  be  contained  in  an  optimal  solution  to  the 
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problem  and  that  the  path  (solution)  should  be  developed  further. 

Two  truncated  enumeration  schemes  for  the  OCPM  problem  dis¬ 
cussed  in  Chapter  III  will  be  presented  here  and  experimental  results 
for  one  of  the  schemes  will  be  given.  The  first  approach,  the  "reduced 
constraint"  algorithm,  assumes  that  the  set  of  precedence  constraints 
are  reduced  by  the  dominance,  feasibility  and  lower  bound  tests  of 
Chapter  tV.  The  second  approach,  the  "fixed  order"  algorithm  assumes 
only  that  project  information  such  as  job  times  and  precedence  rela¬ 
tions  are  avai I  able. 

2.  Reduced  Constraint  Algorithm 

This  algorithm  will  be  based  on  the  assumption  that  a  reduced 
set  of  path  constraints  has  been  determined  from  the  problem  to  be 
solved.  For  the  numerical  example  of  Chapter  IV  the  resulting  con¬ 
straints  are  shown  by  Table  Vl-I  .  If,  for  example,  iob  S  is 
performed,  then  the  minimum  length  of  the  critical  path  would  be  45 
days.  Since  the  lower  bound  on  project  length  calculated  In  Chapter 
IV  was  43  days,  the  reduced  constraint  matrix  assumes  that  will 

be  at  least  43  days. 

The  problem  we  must  solve  then  is 
(I)  Minimize  m  k(i)  _  ^ 

iri  j^l  'J  'J  f  ^ 

Subject  to 

Precedence:  As  given  in  the  reduced  constraint  matrix 
■4-  - 

Oue  Date:  f- Wf  •- 0 
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kU) 

Interdependency:  2  *  l=l,...,h 

J  I  ^ 

plus  "other"  Interdependency  constraints. 


Reduced  Constraint  Matrix 
Table  VI-1 


The  bounding  process  involves  a  separation  of  the  full 
problem  into  two  serially  related  subproblems.  We  first  solve  to  find 
the  minimum  total  job  cost  at  a  particular  point  in  the  tree,  consider¬ 
ing  only  decision  jobs  accepted,  d,j  ac  |,  and  decision  jobs  excluded, 
‘^mn~  o''  jobs 
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accepted  and  the  set  included 

Jl, 

(.)  C,=  2  s  c, 


Subject  to 

ACCEPTANCE  djj  =  I 

EXCLUSION  d„„r;  o 

mn 

Interdependency  k(i) 


S  =  |Sp  > 
mn  ]  E  \ 


plus  "other"  interdependency  constraints. 

An  optimal  solution  to  this  problem  is  the  selection  of  a  set 
of  the  decision  jobs  which  includes  all  of  the  jobs  in  S^,  none  of 
the  jobs  in  S^  meets  all  interdependency  constraints  and  minimizes 
job  cost  Cj.  If  no  such  solution  exists,  then  that  point  of  the  tree 
may  be  labelled  as  infeasible  and  no  further  search  is  required.  If  a 
solution  is  found,  then  we  can  state  that  no  set  of  decision  jobs  can 
be  selected,  given  existing  constraints,  that  will  have  a  lower  total 
job  cost  than  Cj . 

We  have  now  calculated  minimum  job  cost  at  a  particular  point 
iti  the  tree  and  we  now  show  that  it  is  possible  to  ado  to  it  Cp,  which 
is  a  lower  bound  on  penalty  cost  tor  the  date  of  project  completion. 
Given  the  constraints  implied  by  the  jobs  in  and  Sf  and  given 

the  information  from  Table  VI-1,  it  is  possible  to  calculate  a  minimum 
length  for  the  critical  path  of  the  project.  For  example,  if 
contains  job  j  and  S£  contains  job  Sj^  2* 


then  we  must  do  jobs 
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Sq  ,  and  S  (  o  die  i  ”  I  )•  From  Table  IV-I  we  see  that  t-  se 

15,1  j_) 

two  jobs  Introduce  a  path  46  days  long  Into  the  network.  Thus  we  can 

calculate  a  lower  bound,  C^,  of  due  date  penalty  or  premium.  Thus 

the  lower  bound  on  total  cost,  Cj,  for  the  partial  solution  is 

Ci-HC^.  In  addition,  if  It  is  found  that  the  complete  solution  to  the 
j  P 

job  cost  minimization  problem  also  gives  a  path  length  resulting  in 
cost  Cp,  we  may  terminate  further  search  down  the  p  th.  The  optimum 
solution  containing  the  partial  path  has  been  found  since  we  have  a 
complete  solution  with  minimum  job  cost  and  path  length,  given  our 
previous  selections. 

The  general  branching  strategy  to  be  used  here  is  as  follows. 
From  the  set  of  all  decision  paths,  the  set  of  paths  with  the  maximum 
number  of  accepted  and  excluded  jobs,  Pp,  will  be  taken.  For  each  of 
these  path  a  minimum  bound  will  be  calculated  and  the  path  with 

the  minimum  lower  bound  will  be  selected  for  further  elaboration.  Once 
the  path  has  been  chosen,  a  decision  rule  is  required  co  select  deci¬ 
sion  jobs  that  will  be  examined  at  the  next  stage.  The  general 
approach  to  be  used  here  is  similar  to  that  used  by  Eastman  [27l  and 
later  by  Shapi ro  [6/]  in  the  travelling  salesman  problem.  At  each 
node  of  the  problem  they  solve  an  assigi*n«nt  problem  on  the  cost  matrix. 
The  solution  that  results  may  imply  several  sub-tours,  rather  than  a 
complete  all-city  tour.  For  example  In  a  five  city  problem  the  solution 
1-4,  4-1 .  2-3,  3-5,  5-2 

contains  two  subtours,  1-4-1  and  2-3-5-2.  Shapiro  takes  the  smallest 
of  the  tours,  i.e.  1-4,  4-1  and  branches  on  these  variables 
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i nfeasi ble 

In  this  noiat'on  implies  that  the  route  selected  will  include 

/Os 

a  trip  from  city  1  to  city  4  and  \J“V  implies  that  such  a  trip 
will  not  be  made.  Since  the  combination  of  1-4,  4-1  constitutes  a 
subtour,  by  problem  definition  they  both  cannot  appear  in  a  feasible 
solution.  Thus  either  the  condition  "not  1-4"  or  "not  4-1"  (or  both) 
must  hold.  If  it  is  possible  to  close  all  nodes  leading  from  1-4  and 
4-1  ,  then  all  solutions  have  been  discovered  and  bounded. 

This  suggests  the  following  branching  strategy  for  out 
problem.  At  a  particular  rode  of  the  tree,  solve  the  following  integer 
programming  problem  as  in  (2) 


Mi nimize 


S  ■  I,  I' 


Subject  to 


Acceptance 

d,j  =  i 

Exc I  us i on 

d  —  0 

S  ^ 

mn 

mn  / 

I  nterdependency 


k(.) 

>  d. 


-  1 


plus  "other"  interdependency  constraints. 
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A  solution  to  this  problem  will  consist  of  the  cheapest  feasible 
selection  of  the  decision  jobs.  It  is  then  possible  to  test  In  the 
original  decision  graph  for  the  length  of  the  critical  path  associated 
with  this  selection  of  decision  jobs.  Alternately,  If  the  reduced  set 
of  path  constraints  exist,  for  example  see  Table  Vl-l,  the  length  of 
the  critical  path  may  be  determined  directly  from  that  set  of  equations. 
In  addition,  we  may  identify  those  variables  on  the  critical  path  of 
the  network.  These  critical  path  variables  then  form  one  branch  from 
the  new  node,  and  each  variable,  individually  excluded,  forms  the 
alternate  paths  from  the  node.  Any  solution  to  the  problem  must  con¬ 
tain  either  the  set  of  all  decision  nodes  on  the  original  critical 
path,  or  have  at  least  one  of  them  excluded.  If,  for  example,  at  a 
particular  node  on  the  tree,  the  cheapest  feasible  set  of  decisions 
contained  j  and  S|2  |  and  both  of  these  jobs  were  on  a  critical 
path  of  h7  days,  we  would  branch  as  shown  below 


In  the  application  of  the  idea  to  the  DCPH  problem,  the  combination  of 
jobs  on  the  critical  path,  say  j  is  feasible  rather  than 
infeasible  as  in  Shapiro's  problem.  However,  although  it  is  feasible. 
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it  is  never  necessary  to  develop  that  particular  node  further.  Given 
the  exclusions  and  acceptances  to  that  point  on  the  path,  we  solved  a 
progranvni ng  problem  to  minimize  the  total  job  cost,  Cj.  The  minimum 
cost  feasible  set  of  jobs  is  then  tested  in  the  matrix  of  fable  VI-1  to 
determine  project  length,  project  length  cost  C^,  and  the  set  of  jobs 
on  the  critical  path.  Given  that  we  do  the  set  of  jobs  on  the  critical 
path  and  that  prior  acceptances  and  exclusions  be  enforced,  ;.''e  total 
cost  can  not  be  less  than  Cp  +  Cj  and  the  solution  that  gives 
Cj  is  a  feasible  solution  to  the  problem.  Thus  the  node  is  completely 
developed. 

If  at  any  stage  of  the  development  of  the  tree,  all  new  paths 
generated  are  bounded  by  a  previous  complete  solution  or  are  infeasible 
or  if  a  complete  solution  has  been  ootained,  it  is  necessary  to  i.ove 
backwards  up  the  tree  and  pick  new  paths  to  develop.  Again  we  find 
the  unbounded  paths  that  have  the  maximum  number  of  job  acceptances 
and  exclusions  and  cf  these  select  the  one  with  the  minimum  bound.  This 
implies  that  we  move  up  the  tree  to  the  first  open  node.  The  algorithm 
will  now  be  presented  in  flov.'  chart  form,  then  applied  to  the  sample 
problem  of  Figure  lli-1. 

Step  I.  Reduce  the  network  by  the  dominance  feasibility  and  lo\>er 
bound  tests  of  Chapter  IV  to  obtain  the  reduced  constraint  matrix. 

Step  2a,  Solve  an  integer  programming  problem  to  find  the  minimum  cost 
set  of  decision  jobs  that  will  meet  all  interdependency  constraints. 

Set  the  value  of  Cp  to  t^ual  the  cost  of  this  solution.  Go  to  Step  3. 
Step  2b.  Select  the  complete  solution  and  the  job  cost  Cj  calculated 
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at  Step  7  for  the  path  to  be  elaborated. 

Step  3.  Calculate  the  length  of  the  critical  path,  associated 

wi  til  the  decision  Jobs  selected  at  Step  2.  This  may  be  determined 
either  from  the  original  project  graph  or  from  the  reduced  set  of 
constraints.  Let  -rWf"  -I-  pWj:  given  f  —  D.  Thus 

Cp  is  either  the  early  premium  or  late  pen.lty  associated  with  finish 
day  and  due  date  D.  Now  establish  one  branch  from  the  first 

node  which  includes  the  complete  solution  of  Step  2,  that  is  a  particu¬ 
lar  set  of  decision  Jobs,  with  a  total  project  cost  of  Cj -f*  Cj  “  ^p* 
Step  4.  Establish  alternate  branches  from  the  current  node  by  adding 
nodes  which  specifically  exclude  each  of  the  decision  Jobs  on  the 
critical  path  determined  at  Step 

Step  5.  for  each  excluded  decision  Job,  ^ij»  solve  an  integer 
programming  problem,  (2),  which  minimizes  the  total  su-n  of  Job  cost 
C j ,  considering  all  jobs  excluded  to  that  point  on  the  path.  The 
problem  may  have  no  solution  and,  if  so,  the  path  may  be  immediately 
termi nated . 

Step  6.  For  each  excluded  decision  Job,  ^(j.  use  the  Job  alternative 
interdependency  constraints 

y*,  d.  .  =  I  to  see  if  the  series  of 
.'•^1 

exclusions  have  forced  the  “acceptance"  of  some  decision  jobs.  Given 
the  “acceptances",  the  original  project  graph  on  the  reduced  set  of 
path  constraints  will  enable  us  to  calculate  a  minimum  bou’d  on  the 
length  of  the  critical  path,  r^f.  Use  Wf  to  calculate  a  lower  bound 
on  completion  dale  cost,  Cp.  For  each  decision  job  excluded,  add  the 
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job  cost  of  Step  5.  to  the  minimum  completion  date  cost  calculated 

above,  Cp,  to  determine  a  minimum  bound  on  total  cost  Cj  —  Cj  + 

Step  7  Record  the  full  solution  to  the  progrenming  problem  of  Step  5 
and  calculate  the  critical  path  length  associated  with  each  solution. 
Step  8.  Test  all  paths  from  the  last  node  for  feasibility,  complete 
solution  or  bounded  solut'on. 

(a)  If  the  program  of  Step  5  has  no  solution,  then  the  path 
being  tested  is  infeasible  and  no  further  search  on  that  patr  is 
requi red, 

(b)  If  for  any  path  the  bound  of  Step  6  ic  equal  to  the  com¬ 
plete  solution  cost  of  Step  7,  then  we  have  a  complete  solution  to  the 
problem  and  no  lower  cost  /ill  be  found  down  this  particular  path.  If 
the  total  cost  of  the  full  solution  is  less  than  the  cost  of  the 
e/;isting  best  solut 'on,  update  the  existing  best  with  the  new  value. 

(c)  If  the  lower  bound  calculated  in  Step  6  is  higher  than 
an  existing  complete  scluticn  to  the  problem,  the  path  is  bounded  and 
need  not  be  considered  further. 

Step  9.  Choose  from  the  feasible,  unbounded  partial  paths  at  the 
current  node  the  path  with  the  minimum  cost  bound.  If  a  tie  exists, 
chx)Se  the  path  with  tt,  maximum  nu.mber  of  forced  acceptances.  Go  to 
Step  2b. 

If  no  leaslble,  unhoundeu  partial  path  exists,  then  the  node 
is  closed  and  we  backtrack  one  level  in  the  graph.  Go  to  Step  8. 

If  all  nodes  in  the  network  are  closed,  then  the  problem  is 
solved  and  the  current  best  solution  is  the  cptimal  solutior,. 
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3.  Application  of  the  Reduced  Constraint  Algorithm 

The  steps  of  the  algorithm  will  now  be  illustrated  with  the 
DCPM  problem  of  Fi gure  1 1 1-1  with  D  =  45  days ,  r  $20  and  p  =  $40. 
The  complete  solution  tree  for  the  problem  is  shown  in  Figure  VI-2. 

Step  1.  Table  VI-1  illustrates  the  reduced  set  of  constraints  for  the 
sample  problem. 

Step  2.  The  integer  programming  algorithm  of  Chapter  V  will  be  used 
to  solve  the  problem  of  selecting  the  minimum  cost,  feasible  set  of 
decision  jobs.  Figure  Vl-1  shows  the  initial  programming  network 
modified  to  include  the  "other"  interdependency  constraints.  The 
minimum  cost  solution  is  2»  3»  ^12  1*  ^15  2  ^17  1 

with  a  cost  of  Cj  =  $50. 

Step  3.  With  the  Jobs  given  in  Step  2,  we  find  from  Table  Vl-1  that 
the  critical  path  is  56  days  long  and  it  contains  decision  jobs  S^ 

Sjj  2*  ^17  I*  6p  =  (56  -  45)  (40)  =  $440.  The  total  cost  for  the 

solution  Cj.  =  50  -r 440  ~  $490. 


complete  solution 


Step  4.  We  now  establish  alternate  branches  excluding  all  jobs  on  the 
critical  path  of  Step  4. 


FIGURE  VI- i 


exclude  S 


complete  solut'on 


Step  5.  For  each  excluded  job,  determine  the  remaining 

minimum  cost,  C  ,  feasible  solution.  The  values  would  be  obtained 

J 

by  removing  the  node  Sjj  from  the  graph  of  Figure  Vl-1  and  recalcu¬ 
lating  the  shortest  path. 


Step  6,  For  each  exclusion  calculate  the  minimum  bound  for  the  length 
of  the  critical  path.  If  we  exclude  S^  then  we  force  the  accept¬ 
ance  of  no  Job  and  therefore  the  minimum  length  of  the  critical  path 
would  be  the  absolute  niinimum,  43  days.  This  gives  a  reward  of 
(45  "  43)  (20)  =  -$4v1.  If  we  reject  S  ,  then  we  are  forced  to 

I  p  ,  4 

accept  j.  However,  no  path  that  includes  only  decision  Job 

S  is  longer  than  4}  days.  Therefore  tne  minimum  completion  cost 

■  5 1  ^ 

i s  aoai n  -$40 , 


zi  e 
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figure  VI- 
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Step  7>  Compute  the  complete  cost  for  each  node  developed  in  Step  5> 


i  j 

Solutlcn 

'j 

c 

p 

9,3 

S  S 

6,3'  ^.2' 

^12, r 

^15,2'  ^7,1 

150 

52 

280 

430 

15,2 

^^'5,1^ 

^6,3'  ^9,3' 

^12,1* 

^15,7'  ^17,1 

150 

52 

280 

430 

‘l7,l 

^6,3'  S,3’ 

^12,1  ' 

^5,2'  ^17,2 

200 

52 

280 

480 

Step  8.  All  three  paths  S  ,  S  ,  S  are  feasible  and  all 

but  Sj^  j  are  unbounded.  The  minimum  cost  solution  available  costs 


Step  9.  Select  node  Sq  ,  for  elaboration  since  it  is  tie  with  S 

j  1 5 1 2 

for  minimum  bound. 

The  final  solution  as  shown  in  Figure  Vl-2  contains  jobs 
^6,2’  ^12.2’  ^5,2  ^7.1  + '-P  = '^T. 


300  -  hQ  -  $260. 
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k.  Fixed  Order  Algorithm 

The  second  scheme  proposed  works  directly  from  the  matrix  of 
Table  IV>3*  the  reduced  network  matrix,  which  shows  maximal  distances 
between  decision  Jobs  in  a  decision  gr<jph.  For  convenience  this  table 
will  be  reproduced  as  VI-2. 


Reduced  Network  Matrix 


6J  S/  63  9,1  3/  12J  12;z  15,1  15^2  17,1  17,2 

El 

Ss 

16  18  20  22  22  27  36  32  9  13  11  3 

37 

21  17 

17 

V 

19  15 

15 

8  22  18 

23 

9,' 

20 

9/ 

18 

V 

9  13 

18 

'V 

10 

12/ 

10 

10 

15j2 

il  3 

0 

17,) 

5 

1J2 

5 

Table  VI-2 


It  contains  implicitly  a  reduced  set  of  project  paths  that  exist  In  the 
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origindl  network  Figure  lll-l  .  In  this  algorithm  the  order  of 
selection  of  decision  sets  is  strictly  determined,  that  Is  for  our 
problem  an  ordering  of  decision  set  S^,  S^,  Sj2,  Sj^  and  Sj^ 

might  be  chosen.  In  elaborating  a  particular  decision  set,  the  paths 
to  be  considered  are  simply  all  possible  decision  jobs  in  the  set  under 
consideration.  At  Sj  we  consider  branching  on  all  Sjj  ^  j  ^1,2, 
...,k(i}.  This  involves  a  feasibility  test  to  see  i the  new  job  S|j 
is  consistent  with  jobs  previously  chosen  on  the  path  given  the  "other" 
interdependency  constraints.  A  feasible  path  is  then  selected  for 
elaboration  based  on  a  calculation  of  a  total  cost  bound  for  each 
feasible  path.  As  in  the  first  algorithm,  the  path  selected  will 
always  be  from  the  most  fully  developed  paths  so  that  always  push 
directly  down  to  a  complete  solution  or  halt  because  a  path  is  bounded, 
then  backtrack  up  the  tree  closing  open  nodes  as  we  find  them. 

At  any  node  in  the  tree,  it  is  possible  to  calculate 
directly  the  total  cost  of  decision  jobs  assigned  on  the  path  to  that 
point.  If  we  add  to  this  the  total  cost  of  minimum  cost  jobs  in 
decision  sets  from  which  no  assignment  has  been  made,  we  have  a  lower 
bound  on  total  job  cost  C^.  In  our  problems  the  costs  for  jobs  in 
each  decision  set  have  been  normalized,  that  is  the  smallest  cost  is 
subtracted  from  all  costs,  so  that  the  total  minimum  cost  for  unassigned 
sets  will  always  be  zero. 

At  the  node  under  consideration,  it  is  also  possible  to 
calculate  a  minimum  length  for  the  decision  path  under  consideration. 

A  project  graph  can  be  constructed  from  the  information  of  Table  V!'2 
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as  we  have  done  in  Figure  VI -3.  The  reduced  graph  is  broken  at  all 
unassigned  decision  nodes,  and  the  longest  path  algorithm  is  applied 
to  find  the  shortest  possible  length  of  the  network  given  previous 
assignments.  This  value  is  a  lower  bound  on  project  length  and  may  be 
used  to  calculate  — >  minimum  completion  cost.  Cj  »  Cj -h  Cp  gives 

a  lower  bound  for  the  particular  path  under  consideration. 

The  steps  of  the  aigcrithm  are  as  follows. 

Step  I.  Normalize  the  cost  of  all  jobs  in  decision  sets  by  subtracting 
the  lowest  cost  for  any  mob  in  a  set  from  all  jobs  in  that  set. 

Step  2,  Apply  the  routine  of  Chapter  iV,  Section  3.  to  find  a  reduced 
network  containing  only  decision  jobs  and  maximal  length  paths  between 
dec! Sion  jobs. 

Step  3.  Sequence  the  decision  sets  in  a  fixed  order. 

Step  4.  Elaborate  the  first  decision  set  in  the  fixed  order  list,  for 
which  no  assignment  has  been  made  and  calculate  total  job  cost,  Cj, 
for  each  resulting  path.  The  integer  programming  algorithm  of 
Chapter  V  wi H  find  the  minimum  cost  solution,  given  interdependency 
constraints  and  decision  job  acceptances,  if  no  feasible  solution 
exists,  then  search  on  this  path  may  be  terminated. 

Step  5.  For  each  newly  developed  oath,  find  a  minimum  bound  for  the 
length  of  the  critical  path,  and  evaluate  the  cost  of  this  path 

length  as  Cp.  The  reduced  network  from  Step  2  may  be  used  to  compute 
minimum  project  length.  As  in  the  algorithm  of  Chapter  IV,  Section  3, 
we  break  ihe  reduced  network  at  all  decision  jo‘'5,  then  introduce  only 
those  decision  jobs  accepted  on  the  partial  path  which  is  under 
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consideration.  A  critical  path  calculation  which  ignores  all  unaccepted 
decision  jobs  then  gives  a  low«r  bound  on  project  length.  A  lower 
bound  on  the  cost  of  the  partial  path  will  then  be 

Step  6.  Record  the  full  solution  to  the  programming  problems  of  Step  4 
and  calculate  the  critical  path  length  associated  with  each  full  solu¬ 
tion  using  the  methods  of  Step  5.  Compute  the  total  cost  of  each 
solution. 

Step  7.  Test  each  newly  developed  path  for  feasibility,  complete 
solution  and  bounded  partial  solution. 

(aj  A  path  may  be  er^.cluded  based  on  a  feasibility  test  of 

Step  4. 

(b)  If  for  any  path  'he  cost  of  a  complete  solution  from 
Step  6  is  equal  to  the  lower  bouno  calculated  in  Step  5,  no  further 
development  is  necessary.  No  other  so’ut>on  on  this  path  can  have  a 
lower  total  cost.  If  the  total  cost  of  the  full  solution  is  less  than 
the  cost  of  the  existing  best  solution,  update  the  existing  "best" 

wi th  the  new  value. 

(c)  If  the  lower  bound  of  Step  5  is  higher  than  an  existing 
complete  solution  to  the  problem,  the  path  is  bounded  and  need  not  be 
considered  further. 

Step  8.  (a)  Choose  from  the  feasible  unbounded  partial  paths  at  the 

current  node,  the  path  with  the  minimum  cost  bound.  Break  ties  with 
random  selection.  Go  to  Step  3. 

(b)  If  no  feasible,  unbounded,  partial  path  exists,  then  we 
backtrack  one  level  in  the  graph.  Go  to  Step  8(o)  . 


lOA 

(c)  If  no  open  nodes  are  found,  HALT.  The  optimal  solution 
Is  the  current  best  solution. 

5.  Application  of  the  Fixed  Order  Algorithm 

The  steps  of  the  algorithm  will  now  be  illustrated  with  the 
DCPM  problem  of  FI gu re  lll-l  with  D  *  45  days ^  r  ^ $20  and  p  ~  $40. 
The  complete  solution  tree  for  the  problem  is  shown  in  Figures  VI-4,  5- 
Step  1.  The  Job  costs  in  this  problem  are  normalized. 

Step  2.  Table  VI-2  and  Figure  VI-3  Illustrate  the  reduced  network. 

Step  3.  Choose  sequence  Sg,  S^,  Sj2,  Sj^,  Sj^  although  any 
sequence  is  permissible. 

Step  4.  Elaborate  the  three  al  ternati  >res  for  decision  job  S^  and 
for  each  calculate  Cj  the  minimum  Job  cost  for  a  feasible  solution 
containing  Sq  .. 


Step  5.  Calculate  the  minimum  path  length  associated  with  each  partial 
path.  From  Figure  VI-3,  we  can  determine  the  longest  path  through 
Sq  ,  and  only  Sq  ,  is  42  days  long.  Similarly,  bounds  on  _  and 
Sq  ,  may  be  calculated  as  40  and  45  respectively.  The  associated 
values  for  Cp  would  be  -?60,  -$I00  and  0.  Therefore,  the  total 


lower  L)Ounds,  Cj,  are  $3^0,  $50  and  $50. 


Step  6.  The  comolete  solutions  of  Step  h  are 


Acceptance 

Sol  ut 

•on 

Wf 

C-r 

S 

A  1 
^  >  * 

Vr 

S 

12,2 

.  s 

15.2' 

^7.1 

43 

-4Q 

360 

S,2 

'6,3* 

S.2’ 

^12,1 

•  ^15.2' 

^17,1 

150 

52 

280 

430 

^9,3 

^9.3 

^12,1 

‘  ^15.2* 

^17,1 

50 

56 

0 

490 

Step  7.  Choose  either  2  ^9  3  fuitha.-  elaboration  since 

they  have  identical  lower  bounds  of  $>0. 

The  optimal  solution  is  shown  in  Figure  Vl~4  to  be  3^ 
n  -»  1C  -5  '  With  a  total  cost  of  $260. 

7i‘  lift.  17  ,i  '/|l 

6.  Computational  Results  with  the  '^ixed  Order  Algorithcn* 

The  cfficienv.y  of  this  algorithm  for  any  given  problem  will 
depend  on  the  fixed  order  established  for  the  decision  sets.  In  our 


*  The  computational  results  reported  in  this  section  were  obtainad 
with  the  col  labo'’ition  of  M,  Wagner.  The  details  of  four  competitive 
algorithms,  along  with  computional  results,  arc  giver  'n  ^23;.  The 
programs  used  are  listed  in  Appendix  E. 


FIGURE  VI- 
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example,  the  fact  that  S  ,  S,,  and  S  appeared  at  the  head  of 

y  ^  I 

the  list  coupled  with  the  interdependency  constraints  on  these  job 
sets,  allowed  substantial  truncation  based  on  feasibility  tests.  We 
may  also  observe  from  this  example  that  a  combination  of  S,  -  and 
Sj2  I  always  gave  a  good  project  completion  cost  bound  since  they 
allowed  a  long  path  to  form  in  the  network  (52  days).  This,  of  course, 
could  be  observed  directly  in  Table  Vl-l.  This  suggests  that  we  might 
cut  search  time  if  decision  sets  that  appeared  on  long  paths  were 
placed  at  the  front  of  the  fixed  order  list.  This  would  certainly  be 
in  the  spirit  of  the  first  algorithm  presented.  In  order  to  measure 
criticality  of  the  decision  sets,  the  following  procedure  was  used. 

Step  1.  At  each  decision  set,  the  cheapest  decision  job  was  chosen. 

Step  2.  The  resulting  CPM  problem  is  solved  and  the  slack  for  each 
job  calculated. 

Step  3.  Sequence  the  decision  sets  in  increasing  order  of  slack 
measured  for  the  member  of  that  set  in  the  CPh  problem. 

The  algorithm  was  tested  on  10  problems  of  approximately  210 
jobs  and  15  three-job  decision  sets  each.  The  total  number  of  combina¬ 
tions  of  decision  Jobs  would  therefore  be  3*^~  1.5  x  10^.  Character¬ 
istics  of  the  10  jobs  are  given  in  Appendix  A. 

Each  problem  was  srived  twice,  the  first  time  with  a  simple 
technological  ordering  of  decision  sets,  and  the  second  tiine  with  the 
slack  ordering  described  above.  Computational  results  are  given  in 
Table  VI-3.  The  algorithm  was  programmed  in  Fortran  IV  for  an  IBM  7094. 
The  program  was  run  under  a  time-sharing  system  —  the  time  reported 
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derived  from  the  system  Interval  timer  and  do  not  include  swap  times 
imposed  by  tl. ;  time-sharing  system.  Times  required  to  obtain  the 
reduced  network  matrix  are  not  included  here,  but  reported  separately 
in  Appendix  6.  in  addition  to  the  times  reported,  these  problems 
required  2-6  seconds  for  reading  data  and  performing  certain  initiali¬ 
zing  functions  in  preparation  for  the  application  of  th..  branch  and  bound 
procedure.  The  computation  times  reported  are  exactly  repeatable. 

Computation  results  show  an  impressive  superiority  for 
choosing  decisions  based  on  initial  slack.  Ratios  of  computation  time 
as  high  as  50:1  were  found.  In  addition,  the  absolute  amount  of  time 
required  by  the  algorithm  using  a  slack  ordered  list  suggests  that  the 
algorithm  is  erficient  for  these  problems. 


Ill 


Problem 

Number 
(Append! x  A) 

Order!  ng 

Heur i Stic 

Optimal 

Found 

Computation 

Soi ution 

Proven 
time  (sec.) 

78 

Technological 

330-T 

330-T 

S  1  ack 

6.4 

9.8 

79 

Technological 

73.7 

213-T 

S  lack 

16,4 

36.8 

80 

Technological 

51.0 

157-T 

S  1  ack 

13.8 

23.5 

81 

Technologi cal 

139.4 

261-T 

S  1  ack 

101.3 

130.5 

82 

Technologi cal 

7.7 

141-T 

S  1  ack 

12.8 

26.6 

83 

Technologi cal 

14.3 

18.5 

S  i  ack 

1 .8 

2.2 

84 

Technological 

40.8 

73.2 

S  1  ack 

5.8 

9.4 

85 

Technologi cal 

2.2 

5.6 

S  1  ack 

2.0 

2.4 

86 

Technological 

98-T 

98-T 

Slack 

2.2 

2.4 

87 

Technologi cal 

81  .0 

95-T 

S  1  tick 

i  .6 

1  .6 

T  -  computation  terminated  before  completion 


Table  VI-3 
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RESOURCE  CONSTRAINED  DECISION  NETWORKS 

The  decision  net.vorks  we  have  considered  in  previous  chapters 
have  included  precedence  and  interdependency  constraints.  To  evaluate 
a  particular  design,  that  is  a  particular  solution  to  the  problem,  it 
was  simply  necessary  to  calculate  the  length  of  the  critical  path  for 
the  resulting  network.  Then  the  project  cost  could  be  determined 
exactly.  If  resource  constraints  are  added  to  this  problem,  i t  i s  no 
longer  a  simple  matter  to  find  the  mirimum  length  of  tha  project. 

Branch  and  bound  techniques  have  been  proposed  to  solve  the  problem 
exactly,  but  for  problems  of  reasonable  size,  the  computation  time  is 
excessive.  Johnson  (_^2j  reports  that  a  100  task,  single  resource 
problem  ran  79  minutes  without  proving  an  optimal  solution.  In  a  five 
decision  set  problem,  we  would  have  3^  =243  possible  solutions  to 
evaluate  in  this  manner. 

For  this  reason  it  was  decided  to  evaluate  proposed  solutions 
to  the  decision  problem  by  loading  them  under  the  limited  resource  with 
a  "good"  heuristic  rule.  The  heuristic  to  be  used  should  I?  efficient 
so  as  to  keep  computation  time  within  reasonable  limits,  but  it  should 
provide  tight  schisdules.  In  th's  ^Isapter  we  will  examine  nine  possible 
heuristic  leading  ■'ulcs  and  chocsr  the  best  of  them  for  use  in  solving 
resource  constrained  decision  me  dc  I  s  .  «/e  then  examine  three  techniques 

for  choosing  an  optimal  set  of  decisions.  These  include  complete 

1  12 
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enumeration,  pairwise  switching  and  multiple  pairs  switching. 

^ •  Project  Scheduling  Heuri s ti cs 

A  graphical  technique  for  project  scheduling  was  first 
proposed  by  Gantt  [3^.  In  this  technique  each  resource  is  shown  as  a 
bar  on  a  bar  chart  where  tne  horizontal  dimension  is  time.  Each  task 
can  then  be  identified  as  a  rectangle  of  resource  use  continuing  for  a 
specific  length  of  ti.ne.  Precedence  relations  between  tasks  and  a 
loading  heuristics  determine  the  relative  position  of  the  tasks  in  time 
and  the  completion  date  of  the  project.  Given  this  visual  display,  it 
is  possiblejfor  small  projects^to  experiment  with  various  sequences  that 
satisfy  the  precedence  constraints,  so  as  to  determine  the  optimal, 
that  is  minimal  length,  schedule. 

For  large  projects  this  graphical  technique  is  inefficient. 

Not  only  is  it  difficult  to  attempt  many  sequences,  but  it  would  require 
much  time  to  keep  job  i  nfori.ation  up  to  date.  Therefore,  in  most 
current  applications,  a  computer  model  of  the  Gantt  chart  is  maintained. 
With  a  computer  system,  it  is  relatively  easy  to  up-date  project  infor¬ 
mation  as  jobs  are  completed,  calculate  job  slacks  and  test  various 
loading  heuristics.  As  we  have  discussed  in  Chapter  II,  various 
heuristics  have  been  proposed  for  the  locoing  ^.^ob!em  [^3,  25,  43,  44^, 
but  there  has  been  little  comparison  ■'ate  of  tnc  effectiveness  of  these 
rules.  The  rules  to  be  tested  here  are  those  that  can  be  impletnented 
quickly  and  that  have  been  well  regarded  in  previous  work.  They  are 
all  basically  serial  loading  tcchniqt;es.  The  Jobs  in  the  project  are 
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ordered  first  in  technological  order,  then  within  technological  order 
by  a  secondary  measure.  The  loading  routines  take  the  first  job  on  the 
list,  schedule  it  at  Its  early  start  date,  then  proceed  down  the  ordered 
list  scheduling  each  job  in  turn.  The  particular  sequence  of  jobs  in 
the  list  will,  therefore,  strongly  influence  the  completion  date  of 
the  project.  Two  routines  which  will  be  discussed  below  modify  the 
serial  loading  slightly  to  allow  a  previously  scheduled  job  to  be 
shifted  forward. 

Nine  rules  are  to  be  tested  here.  These  consist  of  three 
basic  job  sequences  within  the  technological  order,  each  sequence 
loaded  by  three  heuristic  programs.  The  three  orders  to  be  tested  are 

1 .  random 

2.  i ncreas i n;5  early  start 

3.  increasing  late  start 

The  three  serial  murines,  LOADC  ,  LOADN  and  LOAD  will  now  be  presented. 

LOAOC  takes  each  job  as  it  appears  on  the  technologically 
ordered  list  and  places  it  in  the  schedule  at  the  earliest  possible 
time.  The  start  time  of  any  Job  will  be  constrained  by  the  finish 
time  of  its  predecessors  and  availability  of  resources.  The  flow  chart 
for  the  routine  is  as  follov/s. 

1.  Technologically  order  the  jobs. 

2.  Solve  for  the  critical  path  of  the  network.  Find 
Early  Start,  Late  Start  and  slack  for  each  job. 

3.  Select  the  job  sequence  to  be  tested  (random,  E.S.,  L.S. 
withm  the  technological  order)  and  reorder  the  jobs. 
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4.  Set  i  ■=-  I  where  i  denotes  a  position  on  the  ordered 

job  list. 

5.  Calculate  ES. ,  the  Early  Start  day  of  the  job  in  the 

i position  on  the  list.  Since  the  list  is  technologically  ordered, 
al 1  predecessors  of  the  job  will  be  scheduled  and  Early  Start  can  be 
calculated  in  the  usual  way. 

6.  Attempt  to  schedule  the  job  in  the  I  position  on  doy 
ESj,  If  the  job  cannot  be  scheduled  because  of  insufficient  resources, 
go  to  Step  8.  If  it  can  be  scheduled,  do  so. and  set  i=  i+  1. 

7.  If  all  iobs  have  been  scheduled,  go  to  Step  9.  If  not, 
go  to  Step  5. 

8.  ESj  —  ES.  -t-  1  .  Go  to  Step  6. 

9.  Halt. 

Routin«.s  LOADN  and  LOAD  are  similar  to  LOAnC  except  that  jobs 
previously  scheduled  may  be  shifted  forward  to  reduce  the  den.  ind  for 
resources  o.t  a  given  day.  A  flow  chart  for  LOADN  will  now  be  presented 
>.  Technologically  order  the  jobs. 

2.  Solve  for  the  critical  path  of  the  network.  Find  Early 

Start,  Start  and  Slack,  (SL)  ,  for  each  Job. 

3.  Select  the  job  sequence  to  be  tested  and  reorder  the  jobs 

4.  Set  i  =»  I  where  i  denotes  a  position  on  Che  ordered 

job  list. 

5.  Test  to  sec  if  the  job  in  the  i position  is  scheduled 
If  it  is  scheduled,  set  i  3=.  i  -f- 1  and  go  to  Step  5.  If  :t  is  not 
scheduled,  go  to  Step  6. 
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b.  Calculate  ESj ,  the  Early  Start  day  of  the  i i Job. 

7.  Schedule  the  job  in  the  i^^  position  of  the  ordered 
list  to  begin  on  day  ES|.  If  more  than  the  available  resource  are 
required,  measure  the  excess  resource  required,  then  go  to  Step  9.  If 
sufficient  resources  were  avialable,  set  i  I  -)- I  . 

8.  If  all  Jobs  have  been  scheduled,  go  to  Step  15*  If  not, 
go  to  Step  5. 

9.  List  all  Jobs  scheduled  to  operate  on  day  ESj  which  use 
fit  least  as  much  resource  as  the  excess  resource  .i«asured  in  Step  7- 

Of  these  Jobs,  select  the  i  tt.  maximum  Job  slack,  SL.  Assume  this 

Job  holds  portion  j  on  the  ordered  list. 

10.  Remove  from  the  schedule  any  successor  of  Job  J  which 
has  been  scheduled. 

11.  Set  ESj=  ES,  -hi  . 

12.  Atteoipt  to  schedule  Job  j  or  day  ESj .  If  the  Job 
cannot  be  scheduled  because  of  insufficient  resources,  go  to  Step  13. 

If  It  can  be  scheduled,  do  so  and  set  i  =  J  +  |.  Go  to  Stup  8. 

'3.  S‘-J=  '• 

Id.  ES ,  ES  -hi.  Go  *o  Step  12. 

J  J 

15.  Halt. 

In  the  LOADC  routine  discussed  earlier,  all  resource  conflicts 
were  resoNad  by  shifting  the  job  currently  being  scheduled  forward. 
LOAON  TOdlfles  this  so  that  given  excessive  demands  for  resource  on  e 
given  day.  the  job  with  the  greatest  amount  of  slack  was  shifted 
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forward.  The  routine  LOAD  which  was  tested  is  identical  r.o  LOADN  in 
the  flow  chart  above,  except  that  Step  13  is  omitted.  This  means  :hat 
the  slack  of  a  job  is  not  reduced  as  it  is  shifted  forward. 

The  nine  heuristic  loading  technique'’-  tested  are  illustrated 
by  the  following  matrix. 


JOB 

ORDEll 


LOADING  HEURISTIC 


LOADC 

LOADN 

LOAD 

Random 

L.S.  ! 

.  1 

E.S. 

1 

J 

3.  Experimentation  with  Project  Schedul i ng  Heuristics 


Sixty-five  projects  were  generated  for  this  scries  of  tests. 
These  projects  varied  in  size  from  40  to  230  jobs,  in  length  of  the 
critical  path  from  54  to  over  200  days,  in  scheduled  length,  given 
resource  constraints,  from  6(j  to  over  600  days.  Specific  characteristics 
of  individual  projects  arc  -given  in  Appendix  A.  The  operating  results 
of  the  heuristics  programs  are  given  by  problem  number  in  Appendi;  C. 

The  results  may  be  summori.'cd  briefly  as  fellows.  Of  the 
nine  rules  tested,  the  LOAOC  code  oper.3ting  on  j  LS  ordered  list  was 
clearly  superior  to  all  other  heuristics.  Tables  Vll-I  and  Vi  1-2  show 
that  In  56  of  the  65  problems  tested  ^his  routine  gave  the  shortest 
project  completion  time  and  that  in  47  of  the  56  cases  no  othei  routine 
found  the  minimum  length  schedule.  Furthermore.  Tabic  VI 1-3  shows  that 
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on  no  occasion  was  the  worst  schedule  generated  by  this  particular 
combination  of  Job  ordering  and  loading  routine.  An  examination  of 
particular  project  results  in  Appendix  C  is  interesting.  The  differ¬ 
ence  between  the  best  and  worst  solution  in  many  cases  is  as  high  as 
30  percent  to  SO  percent  of  the  best  .chedule  achieved. 

In  general,  the  late  start  ranking  is  superior  to  either 
random  with  technologically  ordered  list,  or  an  t.S.  ordered  list. 

The  LOAOC  loading  method  is  superior  to  LOAON  and  LOAD  for  L.S.  order¬ 
ings,  but  the  evidence  is  not  so  clear  for  either  a  random  or  an  E.S. 
orde".  The  test  results  clearly  show  that  the  use  of  LQADC  loading 
routine  on  a  list  of  jobs  ordi^red  by  L.S.  wl  I !  provide  good  solutions, 
relative  to  the  other  heuristics  tested.  This  method  wHI  now  be  used 
to  evaluate  alternate  solutions  to  the  resource  constrained  decision 
network  problem. 


LOAOC 

LOAON 

LOAO 

RANDOM 

if 

4 

4 

L.S. 

b6 

10 

2 

E.S. 

4 

3 

2 

Results  of  65  projects 
Number  of  ''best'*  schedules 


Table  Vll-I 
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LOAOC 

LOADN 

LOAD 

RANDOM 

0 

1 

0 

L.S. 

47 

3 

0 

E.S. 

2 

0 

1 

Resul ts  of  65  projects 

Number  of  unique  best  scheduling 

Table  VI 1-2 

LOAOC 

LOADN 

LOAD 

RANDOM 

7 

2 

16 

L.S. 

0 

1 

20 

E.S. 

1 

21 

Results  of  65  projects 
Number  of  worst  schedules 


Table  VI 1-3 

k.  Total  En’jmeratlon 

The  technique  iray  be  explained  as  the  evaluation  of  all 
feasible  combinations  of  tiie  decision  variables  in  a  combinatorial 
problem,  in  general,  the  method  is  nor  useful  because  of  the  large 
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amount  of  computer  tFme  required  to  solve  problems.  We  use  the  method 
here  on  a  series  of  iO  small  problems  in  order  to  h^ve  some  means  of 
evaluating  the  heuristics  we  propose  to  use  on  the  resource  constrained 
decision  network  problem. 

The  problems  to  be  tested  have  three  active  resource  categories, 
40-60  Jobs,  5  decision  sets  of  3  jobs  each.  Thus,  there  are  3^=»  243 
possible  combinations  cf  decision  jobs  to  be  tested  in  each  problem. 

The  routine  iteratively  activates  a  combination  of  five  decision  jobs, 
one  f'om  each  job  set,  then  calculates  the  critical  path  in  the  result¬ 
ing  project  network,  orders  the  jobs  by  late  start  and  toads  then,  under 
specified  resource  iimits.  The  cost  of  the  decision  jobs  in  a  particular 
combination  and  the  resulting  length  of  the  project  are  recorded. 

Table  VI  1-4  reports  the  computation  times  for  the  complete 
enumeration  routines  programmed  in  Fortran  IV  and  run  under  a  time¬ 
sharing  system  on  an  13^  7094.  The  times  reported  are  derived  from 
the  system  intervai  times  and  do  not  include  swap  times  imposed  by  the 
time-sharing  system.  Tite  best  solution  found  is  also  reported. 

As  explained  above,  the  schedule  lengths  are  heuri stical ly 
determined  so  that  for  any  combination  Of  decision  jobs,  we  almost 
certainly  are  not  reporting  the  optimal  schedule.  It  is  our  point, 
however,  that  large  problems  of  the  type  we  are  attempting  to  solve 
cannot  practically  be  solved  by  existing  algorithms.  To  illustrate 
this  point,  we  note  that  since  the  running  time  for  our  3^  problem 
Is  approximately  129  seconds,  a  problem  with  13  decision  sets  would 
requi re  at  least 
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1 

-  .  129  X - 5=2J2  hours 

35  60x60 

if  we  were  to  enumerate  all  possible  solutions  and  solve  them  heuristi- 
cally.  This  number  can  be  multiplied  by  30  to  iOO  if  the  project  has 
60-100  tasks  and  we  wish  an  optimal  loading  for  each  design.  Appendix 
0  reports  all  undominated  solutions  and  a  distribution  of  critical 
path  lengths  and  resource  constrained  schedule  ler;jths  for  the  projects. 

We  now  report  a  technique  that  will  substantially  reduce  the 
search  required  —  but,  as  wi  1 1  be  explained,  will  r<ot  guarantee  that 
the  best  set  of  decision  jobs  are  selected. 

5.  Pairwise  Interchange 

The  pairwise  interchange  technique  is  a  techod  of  parttaTiy 
enumerating  combinations  of  variables  in  0  problem.  We  will  illustrate 
the  method  in  terms  of  a  simple  OCPK  problem.  If  the  decision  network 
contains  three  decision  nodes  Sj,  $2,  Sj  with  respectively  three, 
two  and  two  Job  al  ternati  ves,  then  Figure  VN-1  shows  an  enumerated 
tree  of  all  possible  combinations  of  the  variables. 

A  pairwise  enumeration  scheme  would  begin  with  one  decision 
job  from  decision  set  Sj,  one  from  and  one  from  S^,  In  our 
routine  the  Job  with  the  lowest  cost  in  each  set  is  selected  as  e 
starting  solution.  Assume  the  initial  solution  Is  ^2  I 

|.  The  oairwise  interchange  routine  begins  with  one  of  the  decision 
Jobs,  say  S^,  and  iterates  It  through  ail  possible  alternatives. 
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S  ,  S  and  S  while  holding  al'  other  jobs  as  they  were  In 

1,1  1,2  1,3 

the  initial  solution  ($2  p  t*’®  l^-TSt  decision  job  is 

returned  to  its  original  setting  Sj j  and  then  Job  S2,  then 
are  iterated  through  their  possible  alternatives.  The  combinations  of 
jobs  examined  are  shown  in  Figure  VI 1-2.  Note  that  each  solution  con¬ 
tains  only  one  job  that  is  different  than  the  set  of  jobs  in  the  initial 
solution  and  in  total  we  examine  only  four  solutions  in  addition  to  the 
initial  solution. 

For  each  combination  of  decision  jobs,  the  initial  solution 
plus  the  set  of  four  combinations  generated  by  the  pairwise  interchange 
routine,  the  project  is  scheduled  using  the  heuristic  discussed  in 
Section  2  of  this  chapter.  For  each  of  the  fcur  schedules,  the  sum  of 
job  cost  and  completion  date  cost  is  calculated.  The  combination  of 
jobs  with  lowest  total  cost  is  compared  to  the  cost  of  the  .nitial 
solution  if  it  is  lower,  this  combination  by,  •^mes  the  new  "initial" 
solution  and  jobs  in  that  solution  are  exchanged.  The  process  continues 
until  no  improvement  is  found. 

The  ten  problems  tested  by  complete  enumeration  were  tested 
again  with  the  pairwise  interchange  routine.  The  cost  of  the  best 
solutions  generated,  along  with  computation  times,  are  given  in  Table 
Vli-4.  Complete  details  of  the  solutions  are  given  in  Appendix  D,  A 
review  of  results  shows  that  for  each  of  the  ten  projects,  the  optimal 
solutions  were  found  and  that,  in  every  case,  the  solutions  stepped 
from  one  undominated  solution  to  another. 

For  comparison  with  a  later  heuristic  two,  fifteen  decision 
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set  projects  were  attempted.  The  solutions  obtained  to  these  problems 
are  shown  in  Table  VI 1-4.  It  will  be  noted  that  these  problems  ran  268 
and  465  seconds.  These  times  suggest  that  a  more  efficient  technique 
Is  required  for  large  problems. 

6.  Multiple  Pairs  Switching 

The  pairwise  interchange  routine  when  applied  to  a  five 
decision  set  problem,  three  Jobs  per  set,  examines  ten  alternative 
solutions,  then  accepts  the  best  of  these  for  further  exploration.  In 
some  instances,  several  good  exchanges  are  discovered  on  the  first  pass, 
but  only  the  best  exchange  is  accepted.  Then,  in  the  second  set  of 
pairwise  exchanges,  a  previously  discovered  change  is  found  to  be  good 
and  at  the  end  of  the  second  pass.  It  is  accepted.  This  snggests  that 
a  more  efficient  routine  should,  at  each  stage,  accept  all  exchanges 
that  appear  to  be  beneficial. 

The  routine  developed  proceeds  as  follows,  it  begins  by 
applying  a  simple  pairwise  switching  routine  to  all  decision  variables. 
The  base  solution  used  is  simply  that  which  contains  the  clieapest  Job 
from  each  decision  set.  For  each  new  Job  brought  into  the  solution,  we 
calculate  a  project  length  based  on  the  loading  heuristic  described 
earlier  in  this  chapter.  The  total  cost  of  this  solution  Is  compared 
to  the  previous  lowest  cost  solution  and  the  difference  Is  defined  as 
a  price  for  the  decision  job. 

When  all  decision  Jobs  not  in  the  starting  solution  have  been 
evaluated  in  this  way,  we  have  "prices"  for  all  decision  Jobs.  The 
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implicit  price  of  Jot*s  in  the  originel  solution  is,  of  course,  zero. 

Now,  using  the  algorithm  of  Chapter  IV,  with  joh  prices  as  definsd  here, 
we  can  solve  fer  the  best  set  of  jobs  to  perform.  A  solution  might 
imply  that  orie  or  more  variables  are  tl>  be  tched.  This  design  can 
then  be  evaluated  by  our  loaoing  heuristic  and  a  total  cost  can  be 
calculated.  If  this  cost  is  lower  than  our  previous  optimum,  we 
recalculate  job  "prices"  with  a  pairwise  interchange  routine,  otlierwise, 
we  stop.  Results  of  this  algorithm  for  our  twelve  problems  are  shown 
In  Table  VI I -A  and  in  more  detail  in  Appendix  0. 

7.  Discussion  of  Results 

Table  VI 1-4  compares  the  computation  times  and  quality  of 
solutions  for  the  three  routines  applied  to  the  full  problem.  For  small 
problems  (five  decision  sets),  the  complete  enumeration  routine  takes 
approximately  6-9  times  as  long  as  pairwise  interchange  methods.  Pair¬ 
wise  interchange  techniques  take  approximately  the  same  time  as  the 
multiple  pairs  approach.  !n  eight  of  the  ten  problems  examined,  both 
the  pairwise  interchange  and  the  multiple  pairs  exchange  method  found 
the  optimal  solution  as  proved  by  complete  enumeration.  In  one  of  the 
remaining  two  problems,  the  pairwise  interchange  routine  found  a  better 
solution  than  multiple  pairs. 

For  large  problems  (fifteen  decision  sets),  the  pairwise 
interchange  routine  took  longer  but  found  superior  solutions  to  the 
multiple  pairs  method.  For  these  problems,  it  was  not  practical  to 
enunerate  all  possible  designs  and  so  the  optimal  solution  is  not 
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known  as  cotnputatlon  tltr-^s  for  these  probleTiS  of  6  to  II  minutes  suggests 
that  even  the  simple  routines  ould  not  be  practical  for  large  problems. 


TOTAt  ENUMERATION  PAIRWISE  MULTIPLE  PAIRS 
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Table  Vi  I 


Chapter  VIII 


APPLICATIONS  OF  THE  DCPH  MODEL 

The  decision  nodes  of  the  OCPM  model  have  been  primarily  used 
throughout  this  thesis  to  represent  job  alternatives  in  the  discrete 
time-cost  tradeoff  problem.  This  chapter  will  show  that  the  decision 
nodes  may  also  represent  any  given  Job  performed  at  different  points  in 
time,  or  at  different  physical  locations  and  with  these  interpretations, 
the  model  may  ba  used  to  formulate  the  resource  constrained  project 
scheduling  problem  and  the  single  product  assembly  line  balancing 
problems  as  integer  programming  problems.  In  addition,  the  application 
to  project  time  cost  trade-off  problems  will  be  extended  to  projects 
under  incentive  contracts  with  non-linear  criterion  functions. 

2.  The  mxn  Jch-Shop  Scheduling  Problem 

This  formulation  assumes  that  the  job-shop  problem  has  been 
solved  heurf sticai ly  [20,  2^  and  that  3  feasible  finish  date  the 

early  start  of  artificial  finish  job  has  been  determined.  For  each 

job,  '«j,  it  is  then  possible  to  calculate  an  early  start  f.Sj  and  an 
early  and  a  late  finish  time,  EF|  »LF|,  using  the  usual  rules  of  the 
critical  path  method. 

If  the  criterion  function  is  to  mi(  imize  the  make  span  for 
the  fixed  job  file  under  consideration,  then  job  S;  must  begin  on  day 
ESj  or  later  and  mu..(  finish  on  or  before  LSj.  A  stort  before  ESj 
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is  not  possible  given  precedence  constraints  and  simiiariy  a  finish 
iater  than  LF|  wouid  deiay  the  completion  of  the  project  beyond  W^. 
Since  Wf  is  a  feasible  solution,  a  schedule  that  delays  cannot 
be  optimal . 


Subscripts: 

f  machines 
t  day 
;  job 


f  I  ,2,**a, m 

t  =  1,2 . D 

i  *  i  ,2 . f 


of  ’.Tnediate  predecessors  of 
of  all  Jobs  performed  on  day 


job  ij- 
i  on  machine 


r 


Variables: 


U 


Job  i  beginning  on  day  J 

I  if  the  job  is  performed 
0  otherwise 


j  ■”  ^^i . . 


Constrai nts : 

A 


fd 


ES 


ij 


number  of  machines  of  type  f  available  on  day  d 

time  length  of  Job  i.  It  is  assurwd  that  each  job 
requires  only  one  machine. 

the  start  time  associated  with  each  Job  alternative 
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This  criterion  function  -ifempts  to  minlnulze  the  d.iy  on  which  the 
artificial  finish  Job,  Sf,  begins.  This  effectively  minimizes  the 
day  on  which  all  jobs  are  ftnisi<ed  vn  ai  •  machines. 
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This  formulation  is  related  to  that  of  Bownan  [jbJ  and  may  be 
extended  to  the  resource  constrained  project  in  the  same  way  that  Wiest 
Iso]  has  extended  Bowman's  model.  The  project  formulation  may  also  be 
extended  to  the  resource  constrained  discrete  time-cort  trade-off 
problem  by -expanding  the  decision  set  'or  Cc-ch  job  to  include  several 
(k(i))  job  al ternati ves. 

Wiest  (,80]  has  estimated  that  a  project  with  S5  jobs  in  4 
shops  with  a  time  span  of  3^  days  would  require  S2.2S  equations  and  1650 
variables.  If  job  splits  were  not  allowed,  the  number  of  equations 
would  rise  to  6870.  The  formulation  suggested  here  requites  an  inter¬ 
dependency  constraint  for  each  job  (55),  a  resource  constraint  for  each 
resource  for  each  day  (4  x  30  120)  and  a  constraint  for  each  precedence 

relation  (lOO  -  500).  We  car.  then  estimate  that  this  formulation  would 
require  300  *  700  constraints  —  approximately  onr -tenth  of  those  with 
Bovxnan-Wiest  formulation.  !i  this  formulation  there  would  be  at  least 
one  variable  for  each  d.y  of  slack  in  the  original  heuristic  schedule. 
This  number  (X>uld  be  substantially  below  the  figure  1650,  estimated  by 
Wiest.  An  estimate  in  the  range  50-200  would  not  seem  unreasonable  here. 
The  exact  number  of  constraints  and  variables  will,  of  course,  depend 
on  the  specific  problem. 

It  is  difficult  to  make  an  exact  comparison  with  Manne's 
formulation  [54]  since  his  approach  is  not  suitable  for  the  resource 
constrained  projec.  problem.  He  implicitly  assures  a  resource  level 
of  one  by  his  use  of  non-interference  constraints.  He  does  estimate 
that  -t  job  file  of  ten  tasks  to  be  nci*o>  ed  on  five  machines  would 
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require  250  variables.  There  would  be  a  conftraint  for  each  precedence 
relation  and  two  constraints  for  each  pal r  of  Jobs  which  must  use  every 
machine.  This  would  Involve  approximately  500  constraints  for  his 
problem. 


3.  The  S ingle'^Hroduct  Assembly>Line  Balancing  Problem 

This  formulation  will  again  assume  that  the  combinatorial 

problem  has  been  solved  heurl stical ly  [72,  7?]  and  an  unper  limit 

M  has  been  set  on  the  number  of  stations  to  be  used, 
max 


Subscripts: 

j  -  stations 
I  -  Job 


i  =  1  H 
j  ....  »'^max 

i  =  l....,Sf 


Sjj  -  Job  i  performed  at  station  J 

P.=  ^set  of  immediate  predecesslons  of  job  i  J 


Lj-  -|^set  of  all  Job  alternatives  (Sjj)  that  may  be 
performed  at  station  j .  J" 


Variables: 


d[j  Job  i  performed  at  station  J  j  =  1 . M 

djj  _  ^1  if  job  is  performed  at  station  J 

ip 


max 


otherwi  se 
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Constrai nts : 

C  cycle  time  —  maximum  amount  of  work  a  station  may  perform 
tj  time  length  of  job  1 


Constraints : 

1)  Interdependence 

djj~l  1=1, 

j=l 

each  job  must  be  performed  once  and  only  once. 

2)  Resource  limits 


S;jELj 


tid.j  4C 


1 . M, 


max 


only  C  units  of  time  can  be  performed  at  any 
station 


3)  Precedence 
^max 

^  ■’^pj  ^  j'^ij 

J=1  j=l 


p£Pi 

I  =  1.2 . 


A  job  cannot  be  allocated  to  a  station  unless  all 
its  ,)redecessors  eie  assigned  to  that  station  or 
an  earl ier  one. 


CrI  terion: 


Minimize 


max 


j^I 


i^f.j 
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Here  we  attc.ipt  to  minimize  the  number  of  the  station  in 
v^ich  the  final  job,  S^,  appears.  This  effectively  minimizes  the 
total  number  of  stations  used. 

For  the  II  Joj  ,  rol  lem  of  Jackson  with  C=  10  and  an 
initial  heuristic  solution  of  six  stations,  this  formulation  would 
require  eleven  interdependency  constraints,  six  resource  constraints 
and  fourteen  precedence  constraints  for  a  total  of  31*  The  maximum 
number  of  variables  would  be  66,  but  this  could  be  reduced  by  an  ES-LS 
argument  to  approximately  39.  This  is  a  substantially  smaller  problem 
th^n  competitive  formulations. 

4.  An  Application  of  Decision  CPH  to  Incentive  Contracts 

Recently  several  government  agencies  have  changed  their 
contracting  procedures  from  predominantly  cost-plus- fixed  fee  to 
incentive  fee  contracts.  The  incentive  contracts  are  written  so  that 
the  maximum  fee  obtainable  decreases  as  cost  of  the  project  increases 
[59]  and  the  per  cent  of  the  maximum  fee  actually  paid  decreases  with 
decreasing  performance.  Performance  polr.ts  may  be  awarded  for  success- 
full  performance  or  quality  tests,  and  for  meeting  a  series  of 
specified  due  dates  (or  mile  stones)  within  the  project  network.  A 
sample  contract  fee  vructure  is  shown  in  Figure  Vlll-I. 

A  manufacturer  faced  with  a  time-cost  trade-off  problem 
within  an  incentive  contract  has  an  especially  difficult  problem.  If 
a  Job  is  "crashed",  it  is  possible  that  extra  points  will  be  earned  as 
a  result  of  meeting  a  particular  due  date.  At  the  same  time,  the 
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increase  in  project  cost  will  cause  all  "points"  to  be  slightly 
devalued.  The  following  integer  programming  problem  is  solved 
iteratively  to  determine  the  optimal  selection  of  jobs  to  be  per¬ 
formed.  The  exact  procedure  will  be  described  after  the  model  is 
presented. 

Subscripts: 

i  job  sets  i  s  I , . . .  ,m 

j  decision  jobs  in  a  job  set  j  =  1 ,2 . k(i) 
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RELATION  BETWEEN  FEE  PAID  AND 
PERFORMANCE  ON  POINTS  AND  COST 


\ 

\ 


I’lGUKi'; 


VIII-2 
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Constar.ls: 


■ij 


u 


'm 


P 


m 


'm 


the  time  required  to  perform  S 


tj 


the  cost  of  S 


due  date  given  to  Job  i,. 


B  maximum  project  cost 


point  loss  for  each  day  after  D 


m 


daily  point  reward  for  early  comply 


Constraints: 

interdepsndence 

k(i) 


di .  ==  I 


Pre'edance 

As  In  [19I  If  Sj  precedes  and 
job,  the  preceJ«i)ce  relation  is  shown 
W;  tj  4  W. 


m 


and  if  Sjj  precedes  and  it  is  i 

-Md-djO  -h  t..+  W, ,  4  W 
•j  IJ  IJ  ^  m 


Resource  (Budoet)  Constraint: 

n  k 


n  ^ 

^  S  ^'1  l*^i  •  4-  ® 
i«l  J«I  ^ 


Job  Is  completed 
itlon  of  Job  S„ 


S|  is  a  uni t  set 

I  multi-set  Job 
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Due  Date  Constraints: 

“m  -  “I  t  “»,  =  “m  ^  £  ffl. 

Cri  terlon: 

Maximize  ^  "  Pm^^m 

This  formulation  maximizes  the  number  of  performance  points 
obtained  subject  to  cost  limits  on  the  complete  project,  initially, 
the  problem  Is  solved  without  a  budget  constraint  to  determine  tlte 
maximum  available  points  and  the  cost  associated  with  this  solution. 

B  is  then  set  one  unit  below  the  cost  found  above  and  the  problem  Is 
resolved. 

The  routine  is  applied  iteratively  until  the  minimum  cost 
point  is  reached.  At  each  solution  the  contbination  of  performance 
points  obtained  and  budget  cost  will  allow  total  fee  to  be  calculated. 
Figure  VI 1-3  shows  a  series  of  such  points  with  an  optimal  solution 
marked.  This  approach  has  been  applied  to  an  actual  oroject  by  Smylie 
and  is  reported  by  Crows  ton  and  Smylie  [22j. 


Chapter  IX 


Summary  and  Conclusions 

Our  original  goal  was  to  determine  efficient  procedures  for 
the  solution  of  design  problems  in  situations  where  the  design  could 
only  be  evaluated  in  term'  of  optimal  operating  decisions  for  the 
system.  The  study  was  confined  to  the  area  of  project  scheduling 
although  the  models  we  used  were  related  to  many  other  planning  problems. 
The  relationship  between  these  problems  was  shown  in  the  literature 
review  of  Chapter  II. 

The  articles  reviewed  there  were  categorized  by  the  type  of 
constraint  found  in  the  model.  These  Included  time  or  precedence 
constraints,  showing  the  sequence  In  which  the  jobs  were  to  be  per¬ 
formed,  or  putting  time  constraints  on  the  start  or  finish  time  of  a 
particular  job.  The  use  of  resource  constraints  implied  that  the  jobs 
required  resource  inputs  to  be  performed  and  that  the  availability  of 
the  resource  was  limited  in  each  time  period.  Finally,  interdependency 
constraints  between  Individual  tasks  were  Introduced.  One  particular 
type  of  interdependency,  that  is,  the  mutually  exclusive  relation 
between  sets  of  jobs,  was  developed  In  some  detail  throughout  the 
thesis  since  the  job  alternatives  could  be  used  to  represent  design 
alternatives  in  a  planning  problem. 

Models  with  simple  time  constraints,  the  usual  Critical  Path 
problem,  could  be  solved  optimally  with  longest  path  algorithm.  Problems 
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tliat  involve  simple  interdependency  relations  form  a  special  class  of 
0  -  1  integer  programming  problem.  Tree  search  schemes  have  been 
developed  for  this  class  of  problem.  With  time  and  interdependency 
constraints  in  a  model,  as  we  have  in  DCPM,  the  literature  suggest  a 
sta'^dard  integer  programming  routine  is  required. 

When  resource  demands  and  constraints  are  added  to  a  model,  it 
becomes  a  difficult  problem  to  solve  it  optimally.  For  example,  given 
time  and  resource  constraints  on  a  set  of  tasks,  it  is  a  difficult 
combinational  problem  to  find  the  sequence  of  jobs  that  minimize 
project  length.  Branch  and  bound  techniques  have  been  developed  for 
small  problems,  but  large  problems  must  be  solved  heuri stical ly .  Only 
a  few  models  include  all  three  types  of  constraints  that  we  have  dis¬ 
cussed,  and  these  are  solved  heuri stical ly. 

The  design  problem  that  is  the  central  concern  of  the  thesis 
is  the  problem  of  selecting  which  jobs  to  perform,  from  a  decision  set 
of  mutually  exclusive  alternatives.  The  optimal  set  minimizes  the 
total  of  job  cost  and  project  completion  date  cost.  This  problem  is 
formulated  in  Chapter  Mi  as  an  integer  programming  problem.  It  is 
shown  in  this  chapter  that  the  number  of  precedence  constraints  required 
for  any  network  would  be  one  for  each  precedence  link,  or  one  for  each 
path  in  the  network.  This  number  could  be  so  large  that  the  problem 
could  not  be  solved  by  current  integer  programming  routines. 

In  Chapter  IV  dominance,  feasibility  and  lower  bound  tests 
are  developed  to  eliminate  many  non-binding  precedence  constraints  from 
the  integer  programming  model.  This  is  equivalent  to  the  elimination 
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of  many  paths,  which  can  under  no  condition  become  critical,  from  the 
decision  network.  The  remaining  set  of  paths  which  may  become  the 
critical  path  are  termed  the  "reduced  constraint"  set.  An  algorithm 
is  developed  to  implement  the  dominance  tests  referred  to  above.  In  a 
decision  graph,  the  algorithm  determines  the  longest  path  between  all 
pairs  of  decision  Jobs,  If  any  path  exists,  and  eliminates  all  paths 
but  the  longest  one.  The  result  Is  a  "reduced  network"  containing  only 
decision  Jobs  and  maximal  length  paths  between  decision  Jobs. 

If  we  consider  the  set  of  tasks  In  any  DCPH  problem  and  the 
set  of  interdependence  constraints  defined  on  that  set,  we  observe 
large  numbers  of  mutually  exclusive  relations,  "Job  alternative"  inter¬ 
dependencies,  and  some  number  of  constraints  between  Individual  members 
of  Job  sets,  "other"  interdependencies.  In  Chapter  V  an  integer  program¬ 
ming  network  algorithm  Is  developed  to  efficiently  solve  the  problem  of 
selecting  the  minimum  cost  set  of  decision  Jobs  given  these  types  of 
interdependency  constraints.  The  optimization  technique  used  is  the 
"longest  path"  calculation  of  the  critical  path  method  applied  to  a 
network  in  which  each  patn  through  the  network  is  a  feasible  solution 
to  the  integer  programming  problem.  The  length  of  each  path  is  exactly 
equivalent  to  the  value  of  the  criterion  function  for  the  solution  the 
path  represents.  In  this  algorithm,  the  usual  "slack"  measure  of  the 
critical  path  method  may  be  interpreted  as  the  dual  evaluation  of 
variables  which  are  not  in  the  solution.  When  the  structure  of  inter¬ 
dependence  In  a  problem  does  not  allow  the  algorithm  to  be  directly 
applied,  it  is  Shown  that  It  may  be  coupled  with  a  branch  and  brand 
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algorithm  to  solve  the  problem. 

In  Chapter  VI  two  branch  and  brand  routines  are  developed  to 
solve  1*^6  DCkM  problem.  Both  of  these  use  the  algorithm  of  the  previ¬ 
ous  chapter  to  handle  the  interdependency  constraints.  In  the  first 
routine  at  each  node,  we  calculate  the  minimum  cost  Job  selections, 
given  Job  selections  and  rejections  to  that  point  on  the  tree  and 
determine  the  critical  path  given  these  job  selections.  The  branching 
decisions  are  made  so  as  to  progressively  break  all  critical  paths,  by 
prohibiting  the  use  of  decision  Jobs  on  the  path.  The  second  algorithm 
sets  a  fixed  order  in  which  the  decision  sets  will  be  considered,  and 
at  each  node,  evaluates  all  alternatives  within  that  set.  Given  any 
selection  of  alternatives,  a  lower  bound  on  total  project  cost  is  given 
by  the  sum  of  job  cost,  for  those  Jobs  previously  committed,  and  a 
project  completion  cost,  based  on  a  finish  date  determined  from  the 
reduced  network,  with  only  the  Jobs  that  were  previously  committed 
active  in  the  network. 

The  second  routine  was  coded  and  tested  with  two  different 
orderings  of  the  Job  sets.  The  first  fixed  order  was  simply  a  tech¬ 
nological  order  taken  from  the  position  of  the  decision  sets  in  the 
original  network.  The  second  fixed  order  was  arranged  so  that  the 
decision  sets  were  listed  in  order  of  increasing  slack.  This  ilack 
was  determined  by  choosing  from  each  decision  set,  the  cheapest  decision 
Job  and  solving  the  resultim  critical  path  problem.  The  slack  on  these 
Jobs  was  then  determined  and  its  value  used  to  rank  the  decision  sets. 

In  every  case,  the  slack  order  was  superior  to  a  technological  order 
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and  the  superiority  ranged  from  5^1  to  over  Using  the  slack  order 

problems  of  fifteen  decision  nodes  or  1.5  x  I0^»  possible  designs  were 
solved  in  1.6  to  130  seconds.  This  suggests  that  large  decision  networks 
may  be  solved  efficiently  with  this  routine. 

In  Chapter  VII,  decision  networks  with  resource  limits  were 
introduced.  Before  the  combined  "design",  "operating"  problem  was 
attempted,  it  was  necessary  to  determine  a  good  "operating"  rule,  or,  in 
this  case,  a  heuristic  loading  rule  for  sequencing  project  jobs  under 
limited  resources.  Several  rules  from  the  literature  were  tested  and 
of  these  a  serial  loading  rule,  operating  on  a  job  list  ordered  by  late 
start,  with  no  job  bumping,  proved  superior. 

This  rule  was  then  used  in  combination  with  three  tree  search 
technique^.  These  were  complete  enumeration,  to  provide  proof  of  the 
optimal  solution,  pairwise  Interchange  and  multiple  pairs  exchange. 

The  pairwise  interchange  was  superior  to  mi'>tipie  pairs  routine  in 
that  the  computation  times  were  approximately  the  same  but  the  solutions 
were  superior.  In  ten  problems  the  pairwise  routine  solved  eight 
optimally.  It  was  concluded  that  none  of  the  techniques  were  efficient 
for  large  problems. 

In  the  last  chapter,  the  decision  network  integer  programrlng 
formulation  was  aopi led  to  the  m  x  n  job  shop  scheduling  problem,  the 
single  product  assembly-line  balancing  problem  and  the  OCPH  incentive 
contrast  problem.  For  the  first  applications,  the  formulation  was 
markedly  more  efficient  than  the  one  existing  Integer  programming  formu¬ 
lations  used  t,s  a  basis  of  comparison.  The  lost  formulation  was 
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interesting  in  that  the  inodel  was  applied  to  a  problem  with  a  non-linear 
criterion  function. 


Suggestions  for  Future  Research 

The  research  questions  that  have  arisen  in  this  study  fall 
into  three  broad  categories.  First  we  will  suggest  specific  areas  for 
development  of  algorithmic  and  heuristic  reoutines  for  the  solution  of 
related  problems.  Then  we  will  discuss  the  need  for  additional  tests 
of  the  models  developed  here  on  actual  problems  and,  finally,  we  wi 1 1 
discuss  the  extension  of  these  models  to  other  planning  problems. 

It  is  clear  that  much  more  work  remains  to  be  done  on  the 
development  of  truncated  enumeration  schemes  for  the  DCPM  problem; 

We  have  suggested  two  here  and  tested  one  of  these.  It  would  be 
desirable  to  develop  other  apprO'''-:hes  to  this  problem,  specifically 
one  that  could  operate  on  the  original  network  and  save  the  time 
required  by  our  network  reduction  scheme.  All  models  could  then  be 
tested  against  a  range  of  DCPM  problems.  In  Chapter  V  we  developed 
an  integer  programming  routine  for  restricted  types  of  problems,  and 
suggested  that  it  could  be  -''uplcd  with  a  tree  starch  scheme  for  more 
complex  problems.  It  would  be  useful  to  develop  a  program  for  such  e 
combined  trodel  and  test  i t  on  a  scries  of  problems. 

The  area  that  requires  the  most  work  would  be  the  solution 
of  resource  constrained  decision  networks.  We  have  attempted  several 
methods  that  give  reasonable  solutions,  but  use  excessive  amounts  of 
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time.  A  routine  patterned  after  the  first  branch  and  bound  scheme  of 
Chapter  VI  would  have  some  hope  of  success.  In  this  application,  the 
critical  sequence  as  defined  by  Wiest  would  appear  to  be  useful  in  the 
branching  rule,  rather  than  the  critical  path  used  here. 

One  application  of  the  DCPM  model  Cbs]  has  been  attempted 
and,  in  that  instance,  the  netvK>rk  reduction  rules  suggested  here  were 
extremely  powerful  in  reducing  the  size  of  the  network  that  had  to  be 
considered.  This  experience  suggests  that  many  real  networks  have 
large  areas  with  substantial  slack,  so  that  large  reductions  is  possible. 
If  this  is  true  in  general,  then  the  network  reduction  notions  of  this 
thesis  may  be  extremely  useful  to  managers  even  though  they  do  not  use 
full  range  of  optimization  techniques  suggested  here.  If  a  large 
problem  is  reduced  to  a  small  one,  then  the  manager's  heuristics  may 
perform  better. 

Finally,  we  believe  that  these  models  will  be  usefui  at  the 
detail  level  of  job-planning,  that  is,  to  the  optimization  of  man- 
machine  and  related  process  charts  and  to  the  most  aggregate  level  of 
policy-making  within  a  firm.  Wherever  there  is  a  variety  of  things  to 
be  done,  connected  with  the  types  of  constraints  we  have  discussed 
throughout  this  thesis,  our  models,  or  some  version  of  them,  should  be 
relevant.  We  hope,  then,  that  this  research  will  influence  the 
development  of  a  wide  range  of  planning  models. 


APPENDIX  A 


149 


APPENDIX  A 


Problem  Generation 

The  problems  used  throughout  this  thesis  were  generated  from  six 
basic  precederce  networks.  We  shall  label  these  Type  A,  B,  C,  D,  E  and 
F.  All  problems  of  type  A  contained  an  Identical  tree  structure  (precedence 
ordering)  and  Identical  decision  Jobs.  The  Job  lengths,  Job  costs  and 
resource  usages  for  all  decision  Jobs  In  all  projects  were  predetermined. 

All  these  nodes  contained  exactly  three  alternatives.  All  other  Job 
times  and  resource  usages  were  randomly  determined  for  all  non-declslon 
Jobs  In  all  problems.  Finally,  the  Cw5t  of  all  non-declslon  Jobs  was  set 
at  zero  since  these  costs  were  constant. 

The  subroutine,  HYPO,  used  to  generate  the  final  project  from  the 
basic  networks  is  shown  at  the  end  of  this  Appendix.  The  following  table 
gives  some  data  on  the  precedence  networks.  This  Is  number  of  unit  Jobs, 
number  of  decision  sets,  number  of  precedence  links. 


Network 

Jobs 

Decision 

Sets 

Precedence 

Links 

A 

34 

5 

72 

B 

46 

5 

87 

C 

49 

5 

92 

D 

49 

5 

92 

E 

171 

15 

364 

F 

202 

15 

343 

Table  A-2  will  present  details  on  projects  used  throughout  the  thesis. 
The  Information  Includes  network  type,  ctlClcsl  path  (given  cheapest  deci¬ 
sion  Jobs),  resource  limit,  total  resource  ussge  (cheapest  decislcn  Job 
solution)  and  project  deadline  Information  where  relevant. 
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PROJECTS 


Problem  Type 
No. 

Critical 

Path 

Resource 

Limit 

Total 
Res.  1 

Resource  Usage 
Res.  2  Res. 

1 

A 

61 

12 

899 

514 

579 

2 

A 

50 

12 

715 

403 

463 

3 

A 

57 

18 

747 

332 

266 

4 

A 

66 

18 

919 

85 

81 

5 

A 

70 

18 

815 

217 

68 

6 

A 

63 

12 

891 

580 

654 

7 

A 

53 

15 

526 

286 

265 

8 

A 

64 

12 

1,031 

225 

263 

9 

A 

62 

15 

875 

93 

56 

10 

A 

62 

12 

765 

523 

580 

11 

A 

58 

12 

769 

67 

18 

12 

A 

54 

12 

713 

549 

410 

13 

A 

67 

12 

736 

206 

102 

14 

A 

62 

12 

519 

199 

142 

15 

A 

62 

12 

750 

58 

48 

16 

A 

60 

12 

942 

674 

678 

17 

A 

58 

12 

723 

208 

2';6 

18 

A 

58 

12 

709 

83 

36 

19 

A 

63 

12 

928 

749 

638 

20 

A 

64 

12 

538 

162 

403 

21 

A 

56 

12 

785 

0 

94 

22 

A 

55 

15 

766 

4b  8 

717 

23 

B 

75 

12 

1,254 

713 

952 

?4 

B 

74 

15 

1,027 

775 

811 

25 

B 

63 

18 

999 

7^0 

700 

26 

B 

94 

12 

861 

296 

0 

27 

B 

70 

15 

732 

346 

199 

28 

B 

88 

18 

934 

307 

364 

29 

B 

67 

12 

851 

106 

143 

30 

B 

79 

15 

973 

164 

117 

151 


Problem  TVpe 
No. 

Critical 

Path 

Re80>trce 

Limit 

Total 
Rea.  1 

Reaource  Uaage 
Rea.  2  Rea. 

31 

B 

r,3 

18 

1,279 

78 

114 

32 

C 

67 

12 

1,355 

988 

974 

33 

C 

61 

12 

696 

302 

351 

34 

C 

63 

12 

1,051 

130 

121 

35 

C 

67 

15 

914 

601 

738 

36 

C 

69 

15 

1,044 

0 

232 

37 

C 

61 

15 

999 

136 

121 

38 

C 

58 

18 

1,075 

690 

782 

39 

C 

63 

18 

600 

427 

412 

40 

C 

63 

18 

985 

185 

122 

41 

D 

76 

12 

1,392 

785 

1,034 

42 

0 

63 

12 

1,082 

432 

510 

43 

E 

213 

12 

4,723 

4,230 

4,063 

44 

3 

194 

10 

4,700 

578 

579 

45 

E 

195 

18 

4,561 

298 

357 

46 

E 

192 

12 

4,164 

3,954 

3,265 

47 

E 

195 

18 

2,456 

1,666 

1.229 

48 

E 

202 

12 

2,333 

1,695 

1,726 

49 

E 

207 

12 

4,268 

341 

342 

50 

E 

201 

15 

4,  m 

4,048 

4,054 

51 

E 

221 

15 

2,185 

1,397 

1,996 

52 

E 

207 

15 

4,350 

560 

467 

53 

E 

200 

18 

5,101 

4,445 

4,043 

54 

E 

218 

12 

5,849 

3,022 

2,465 

55 

E 

218 

14 

3,936 

4,142 

3,846 

56 

S 

210 

12 

3,645 

3,258 

3,158 

57 

P 

102 

12 

4,397 

3,610 

3,536 

58 

F 

95 

18 

5,227 

661 

684 

59 

f 

109 

12 

2,322 

1,627 

1,808 

60 

T 

118 

12 

4,467 

456 

491 

152 


Problem  Type 
No. 

Critical 

Path 

Resource 

Limit 

Total 
Res .  1 

Resource  Usage 
Res.  2  Res.  3 

Reward 

Due-Date  Penalty 

6: 

F 

100 

15 

4,383 

3,616 

3,626 

62 

F 

95 

15 

2,320 

1,429 

’  ,594 

63 

F 

100 

15 

4,401 

485 

421 

64 

F 

113 

18 

4,638 

4,086 

3,953 

65 

F 

103 

12 

2,226 

1,223 

1,565 

66 

D 

69 

18 

1,155 

466 

732 

40 

70 

50 

67 

B 

87 

10 

1,079 

522 

543 

50 

135 

50 

68 

D 

67 

10 

1,212 

29S 

269 

30 

145 

100 

69 

D 

76 

12 

1,392 

785 

1,034 

30 

142 

60 

70 

C 

76 

10 

1,111 

494 

621 

50 

134 

35 

71 

3 

78 

1  e 

1,191 

1,054 

895 

50 

105 

80 

72 

B 

63 

18 

99« 

770 

700 

10 

50 

100 

73 

C 

67 

12 

1,354 

988 

974 

10 

160 

20 

74 

A 

53 

10 

523 

118 

0 

10 

66 

200 

75 

B 

63 

18 

999 

770 

700 

30 

70 

oO 

76 

F 

102 

12 

4,397 

3,610 

3,536 

30 

470 

70 

77 

E 

202 

12 

2,333 

1,695 

1,726 

30 

280 

50 

78 

E 

- 

- 

- 

- 

- 

150 

136 

50 

79 

E 

- 

- 

- 

- 

- 

150 

186 

50 

80 

E 

- 

- 

- 

• 

- 

150 

187 

50 

81 

E 

- 

- 

- 

- 

- 

150 

195 

50 

S2 

E 

- 

- 

- 

- 

- 

150 

185 

50 

F 

- 

- 

- 

- 

- 

150 

150 

50 

84 

F 

- 

- 

- 

- 

- 

150 

90 

50 

35 

F 

- 

- 

- 

- 

- 

150 

97 

50 

86 

p 

- 

. 

- 

- 

150 

92 

50 

87 


150 


100 


50 


HYPO  MADTRN 

04/21  1356.5 

00010 

SUBROUTINE  HYPO  (CP, 1 P, K,TIM) 

00020 

DIMENSION  CP(260,25), IP(50,15),K(30),TIM(15) 

00030 

PRINT  4 

00040 

4 

FORMAT  (35H  PUNCH  (SEED  1 5 , VRES , VMl, T 1 M5- 7F5 . 

00050 

READ  1, ISEED,VRFS,VM1,(TIM(I ), 1=11,15) 

OOOGO 

1 

FORMATC I5,7F5.3) 

00070 

print  37 

00080 

37 

FORMAT(13H  PUNCH  K14-I5) 

00090 

READ  1S,(K( 1 ), 1=1,14) 

00100 

13 

F0RMAT(14I 5) 

00110 

PRINT  3 

00120 

3 

FORMATdSH  PUNCH  IFILE  15) 

00130 

READ  5, 1 F 1 LE 

00140 

5 

FORMAT  (15) 

00150 

DO  7  11=1,300 

00100 

REAn(IFILE,G)(TIM( 12), 12=1,10) 

00170 

6 

FORMAT  (10F5.0) 

00180 

1  F(TI.M(1)  .EQ.O. )  no  TO  9 

00190 

C?( 11,12)=! 

00200 

IF  (TIM(2) .EQ.O. )  GO  TO  16 

00210 

CP( ll,12)=3 

00220 

16 

K(1)=K(1)+1 

00230 

DO  S  13=1,10 

00240 

CP( 1 1, 1 3)=TIM( 15) 

00250 

8 

CONTINUE 

00200 

7 

CONTINUE 

00270 

9 

CONTINUE 

00280 

REWIND  IFILE 

00290 

CALL  SETUF( ISEED,USEED) 

00300 

C 

PICK  UP  DJ  INFO 

00310 

1  2  =  1 

00320 

K1=K(1) 

00330 

1 p: I 2,1)=1 

003G() 

DO  10  11=1, K1 

00350 

CP( I 1,1S)=1 

00300 

CP( ll,22)=CP( 11,1) 

00370 

IF(CP( ll,12),E0.1.)  no  TO  10 

00580 

1  2=12  +  1 

00400 

! P( 1 2,2)=CP. 1 1,1) 

00410 

IP(I2,15)=CP(I1,2) 

UU420 

1 P( 1 2, 5)=CP( 11,2) 

00430 

1 P{ 1 2, l)=i 2 

00't40 

1 P( I 2, 14)=CP( 11,12) 

00  45  0 

CP( 1 1, 25)=I 2 

OOliGO 

10 

CONTINUE 

0  0  4  7  0 

Kf 5  )  =  l 2-1 

0  0  480 

C 

NOW  PREP  RFLATION-ASSUME  ALL  3 

00  490 

1  P(2,4)  =  l 

00500 

1 P(3,4)=l 

00510 

1 P(4,4)=l 

005  20 

IS=I 2-5 

00530 

DO  15  1 5=2, IS, 3 

C0540 

16=15+3 

0  (i  5  5  0 

18=15+5 

00560 

DO  12  14=16,18 

005  70 

DC  11  15=1,5 

0  05  8  0 

17=15+5 

00590 

IP( 14, I7)=I5+I5-1 

00600 

11 

CONTINUE 

153 


oo&io 
0OG2O 
00C5U 
OOG'jO 
00  65  0 
OGGGO 
00G70 


12  CONTINUE 

13  CONTI  NUT: 

ir=i 2+1 

DO  14  11=1/5 

12=11+5 

IP( IF/ I 2)=IF-4+ll 

14  CONTINUE 
IP( IF/1)=IF 

TIM(1)=.07 
TIM(2)=.21 
TIM(5)=.54 
TIM(4)=. 67 
T  1 1 1 ' 5  )  =  .  7  0 
T 1 11  {  G )  =  .  8  6 
T1M(7)=.02 
TIM(8)=.07 
TIM(D)=.00 
Tni(10)=1.0 
KlLl  =  Kl-'i 
DO  330  1=2, Kill 
1F(CP( 1 ,12) .MF.l .)  00  TO  530 
CALL  RANNOFCUGFED/VRAN) 

DO  320  ITfl  =  l,15 

T  =  Tir';(lTM) 

IF(VRAr'-T)  321/  321/320 

320  CONTINUr 

321  CONTINUE 
CP( 1 /3)=1 IN 

330  CONTINUE 

T  !  M  ( 1 )  =  .  0  4 
T  1  N  (  2 )  =  .  1 2 
TIf1(3)  =  .22 
Tlf:(4)  =  .34 
T  I  f  K  5  )  =  .  4  8 
T!M(C)=.G4 
T  1 1 K  7 )  =  .  7  8 
T I  ^'  (  8  )  =  .  8  8 
T  1  M  (  0 )  =  .  0  6 
TI!I(1U)=1.0 
T I  (1 5  )  =  1  .  0 
DO  3  8  0  1=2/KU1 
IFCCI'C  ,12).NE.l.)  CO  TO  330 
CALL  r.niNOFfUSFEP/VRAN) 

I F ( VRAM-VRES )  370,450,450 
570  rONTIf^Ur 

CALL  ilAUr.'OFd'SFFn,  VRAM) 

no  3';o  102  =  11,15 

T=Tir( 102) 

IF(VRAN-T)  4r,C,  400,3''0 
300  CONTlNUr 
|,0G  CONTINUE 
102  =  10.2-4 
MO  CONTI  N'LT 

TAi  L  i.;M.'rOF(U5FPP,vrAN) 
pn  Iv?()  !Tfl=l,10 
T=Tir( ITU) 

ir(VF'AN-T)  4  3  0,430,4  20 
It?')  CONTINUE 
1;3'0  CONTINUE 


03200 

CP(  Mr:2)  =  ITf! 

01210 

CALL  PANNOF(ur, ECO, VRAM) 

01220 

1  F(VRAN-Vfll)  370,380,380 

01250 

450 

CONTI MUE 

01240 

380 

CONTINUE 

01250 

PRINT  17,K1 

01260 

17 

FORMAT  (15) 

01270 

C 

DO  15  1=1, K1 

01280 

C 

PRINT  ri00,(CP(  l,,)),.J  =  l, 

01290 

COO 

FORMAT(i5F5.0) 

01500 

15 

CONTINUE 

01310 

■RETURN 

01320 

END 

n  5.283+1.56C 
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COMPUTATIONAL  RESULTS  -  REDUCED  M':TM0RK  ALGORITHM 


The  iterative  application  of  the  longest  path  algorithm  as  des¬ 
cribed  in  Chapter  IV-3  gives  maximal  distances  between  all  decision  jobs. 


The  time  required  for  this  routine 
NETWORK 

A 
B 
C 
D 
E 
F 


on  six  baxlc  networks  is  shown  below. 

COMPUTATION  TIME 
Sec. 

7.7 

12.5 
15.3 
14.7 

38.5 
35.2 


A  sample  output  from  the  program  for  network  A  follows: 


Decision  Jobs 


From  To 

KSg)  3 

1  4 

1  5 

1  9 

1  10 

1  11 

1  12 

1  13 

I  14 

1  26 

1  27 

1  30 


Maximum  Path  Length 
(days) 

1 

1 

1 

11 

11 

11 

8 

8 

8 

14 

14 

33 
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Decision  Jobs 

Maximum  Panh  Length 

From 

To 

(days) 

3 

26 

11 

3 

27 

11 

3 

30 

11 

3 

30 

4 

26 

16 

4 

27 

16 

4 

30 

16 

4 

35 

5 

26 

26 

5 

27 

26 

5 

30 

26 

5 

38 

45 

9 

^F 

26 

10 

^F 

29 

11 

"f 

31 

12 

22 

16 

12 

23 

16 

12 

24 

16 

12 

Sp 

35 

13 

22 

17 

13 

23 

17 

13 

24 

17 

13 

^F 

36 
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Decision 

Jobs 

Maximum  Path 

From 

To 

(Hays) 

14 

22 

18 

14 

23 

18 

14 

24 

18 

14 

37 

22 

^F 

11 

23 

13 

24 

21 

F 

26 

15 

27 

17 

30 

^F 

20 
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APPENDIX  C 

PROJECT  LOADING  HEURISTIC  RESULTS 


problem 

No. 

Order 

Load”C 

Load-N 

problem 
Load  No, 

Order 

Load-C 

Load»N  Load 

1 

Random 

L.S. 

E.S. 

lOA 

97 

109 

104 

113 

112 

124 

113 

110 

12 

Random 

L.S. 

E.S. 

59 

59 

62 

62 

61 

61 

61 

61 

62 

87 

95 

105 

13 

67 

67 

67 

4. 

79 

82 

97 

67 

67 

67 

78 

92 

92 

68 

67 

67 

66 

64 

64 

14 

77 

74 

69 

3 

60 

76 

79 

63 

63 

71 

60 

66 

62 

68 

63 

80 

A 

94 

83 

83 

15 

78 

90 

94 

78 

83 

87 

62 

79 

67 

93 

95 

99 

74 

74 

78 

5 

84 

70 

70 

81 

70 

16 

94 

76 

89 

79 

89 

86 

73 

71 

70 

73 

93 

84 

6 

111 

102 

107 

111 

112 

119 

115 

128 

130 

17 

66 

60 

63 

62 

60 

66 

62 

62 

66 

7 

63 

57 

63 

61 

58 

58 

61 

61 

61 

18 

62 

58 

59 

58 

63 

58 

58 

60 

58 

8 

92 

84 

86 

84 

93 

116 

19 

117 

116 

138 

120 

144 

129 

86 

102 

146 

118 

120 

150 

9 

111 

96 

92 

104 

99 

99 

110 

lOl 

83 

20 

75 

73 

74 

78 

73 

83 

80 

92 

88 

10 

87 

71 

86 

84 

93 

97 

2i 

96 

87 

99 

90 

99 

133 

78 

80 

105 

92 

89 

89 

11 

76 

72 

69 

88 

73 

73 

83 

73 

103 

22 

86 

76 

78 

86 

32 

95 

91 

112 

107 

162 


Problem  Problem 


No. 

Order 

Load-C 

Load-N 

Load 

No. 

_  Order 

Load-C 

Load-N 

Load 

23 

Random 

181 

178 

184 

35 

Random 

95 

90 

92 

L.S. 

148 

167 

167 

L.S. 

84 

89 

89 

E.S. 

157 

162 

172 

E.S. 

89 

98 

98 

24 

128 

113 

131 

36 

89 

91 

91 

98 

112 

140 

87 

90 

89 

115 

109 

145 

96 

98 

122 

25 

101 

96 

119 

37 

82 

88 

87 

77 

83 

97 

81 

'  t 

O'. 

85 

80 

95 

110 

88 

94 

127 

26 

138 

122 

122 

38 

79 

90 

94 

102 

98 

123 

75 

83 

34 

116 

123 

177 

83 

96 

94 

27 

87 

87 

97 

39 

64 

63 

63 

72 

74 

74 

63 

63 

03 

87 

84 

85 

64 

63 

63 

28 

99 

92 

92 

40 

68 

68 

68 

92 

92 

92 

68 

71 

73 

99 

99 

99 

77 

81 

77 

29 

113 

122 

141 

41 

167 

160 

174 

96 

103 

122 

153 

160 

189 

97 

97 

97 

170 

171 

232 

30 

107 

108 

108 

42 

131 

131 

133 

CO 

84 

102 

111 

119 

156 

100 

106 

112 

132 

143 

138 

31 

123 

112 

Hi 

^3 

603 

600 

632 

83 

97 

99 

568 

595 

911 

100 

104 

109 

614 

613 

817 

32 

170 

206 

209 

44 

627 

623 

729 

155 

162 

159 

601 

619 

673 

167 

169 

170 

611 

640 

764 

33 

89 

94 

94 

45 

337 

370 

351 

72 

81 

83 

298 

329 

431 

81 

“6 

IOC 

327 

344 

355 

34 

120 

117 

137 

46 

545 

513 

827 

109 

115 

175 

500 

564 

723 

112 

107 

130 

534 

549 

639 
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problem  Problem 


?Jo. 

Order 

l.ond-C 

Lon'»-N 

I.ond 

No. 

Order 

LOrid-C 

Load-N 

Load 

47 

Random 

233 

223 

255 

59 

R.indom 

261 

267 

321 

L.S. 

211 

203 

214 

L.S. 

237 

265 

324 

E.-S. 

230 

241 

242 

L.S. 

250 

298 

321 

48 

307 

308 

341 

60 

466 

445 

651 

294 

309 

382 

432 

428 

489 

313 

330 

368 

430 

456 

559 

49 

451 

476 

552 

61 

379 

410 

493 

417 

454 

607 

363 

358 

454 

433 

490 

644 

369 

379 

486 

50 

435 

447 

522 

62 

197 

245 

230 

416 

448 

520 

184 

199 

235 

440 

497 

540 

190 

193 

271 

51 

262 

256 

258 

63 

333 

371 

378 

245 

247 

253 

321 

359 

441 

271 

276 

288 

325 

348 

388 

52 

353 

385 

403 

64 

323 

306 

410 

338 

347 

562 

318 

343 

333 

353 

379 

493 

333 

357 

382 

53 

398 

399 

482 

65 

153 

183 

195 

349 

409 

490 

147 

164 

176 

368 

415 

47S 

149 

188 

188 

54 

543 

563 

625 

531 

583 

762 

5  56 

592 

713 

55 

453 

488 

663 

410 

4  7 1 

573 

464 

509 

574 

56 

4HH 

470 

56  5 

44  7 

4  73 

488 

5  78 

.84 

■■i7 

4'M) 

551 

■.07 

4. M2 

.82 

t.38 

50  7 

529 

1. 1  8 

:i;  n-. 

.'*,1  3111 

i  :’/H 
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APPEH^DIX  D 

COMPUTATIONAL  RESULTS  -  RESOURCE  CONSTRAINED 
DECISION  NETITORK  PROBLEMS 

This  Appendix  reports  experimental  results  for  three  routines, 
total  enumeration,  pairwise  exchange  and  multiple  pairs  exchange  applied 
to  resource  constrained  decision  network  problems. 


Undomlnatcd  Solutions 


Critical  Path 

Project  Length 

Cost 

Decision 

Jobs 

69 

79 

250 

^1.3 

^2,3 

^3,3 

^.3 

^5,3 

62 

76 

300 

^1,2 

^2,3 

^3,3 

"4,3 

^5,3 

62 

73 

350 

^1,2 

^'2,3 

^3,3 

^4,2 

S.3 

50 

71 

500 

^1.1 

h,2 

^3,3 

^4,3 

^^5,2 

50 

70 

550 

^1.1 

h,2 

^3,3 

^4,2 

®5,2 

49 

69 

850 

^1.1 

^3,1 

^4,3 

^5,1 

Critical 

Path 

Proj  ect 

Length 

Lciigch 

Number 

Length 

Numbi 

49 

24 

69 

4 

50 

24 

70 

9 

52 

60 

71 

21 

58 

18 

72 

44 

60 

18 

73 

25 

62 

18 

74 

30 

69 

54 

75 

26 

72 

27 

76 

13 

77 

24 

78 

39 

79 

8 

Time  259  sec. 

Total  Enumeration 

Problem 

66 

16^ 


Undomlnaced  Solutions 

Critical  Path  Project  Length  Cost  Decision  Jobs 


87 

153 

250 

^1.3 

^2,3 

^3,3 

=1..3 

^5,3 

87 

149 

300 

®1,3 

^2,3 

h,2 

®4,3 

^.3 

87 

145 

350 

®1.2 

^2,3 

^3,2 

®4,3 

^5.3 

87 

143 

450 

^.3 

^2,3 

^3,2 

^4.3 

^5.3 

87 

137 

550 

'>1.1 

^2,3 

■‘^3,1 

®4,3 

^5,3 

Critical  Path  Project  Length 


Length 

Number 

Length 

Number 

68 

162 

137 

1 

153 

15 

87 

81 

138 

1 

154 

18 

139 

4 

155 

13 

140 

2 

156 

14 

141 

4 

I*;? 

6 

142 

2 

158 

11 

143 

5 

159 

5 

144 

4 

160 

10 

145 

10 

161 

6 

146 

10 

162 

8 

147 

15 

163 

3 

148 

12 

164 

3 

149 

13 

165 

1 

150 

12 

166 

2 

151 

13 

167 

1 

152 

15 

168 

2 

169 

1 

170 

1 

Time  '*10  sec. 


Total  Enumeration  Problem  67 
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Undomlnated  Solutions 


Critical  Path  Projec 

t  Length 

Cost 

Decision  Jobs 

67 

169 

250 

‘^.3 

‘^2,3  ^3,3 

^4.3 

^5.3 

57 

169 

300 

^1,2 

^2,3  ^3,3 

^4.3 

^5,3 

57 

152 

350 

^1,2 

^2,3  ^3,3 

^4,2 

S.3 

55 

150 

400 

^1,2 

^2,3  ^3,2 

^4,2 

^5,3 

57 

147 

450 

^1.1 

®2,3  ^3,3 

^4.2 

^.3 

55 

145 

500 

^1.1 

^2,3  ^3,2 

^4.2 

^.3 

53 

143 

600 

^1.1 

^2,3  ^3,1 

^4,2 

S,3 

53 

142 

750 

^1,1 

^2,1  ^3,1 

^4.2 

S.3 

Critical  Path 

Project  Length 

Length 

Number 

Length 

Number 

47 

72 

142 

2 

:6 

143 

2 

156 

6 

53 

18 

144 

5 

157 

14 

55 

18 

145 

5 

158 

5 

57 

18 

146 

6 

159 

15 

67 

54 

147 

8 

160 

8 

68 

27 

148 

8 

161 

9 

149 

11 

162 

9 

150 

13 

163 

6 

151 

10 

164 

9 

152 

15 

165 

8 

153 

9 

166 

5 

154 

13 

167 

9 

155 

11 

168 

3 

169 

7 

170 

3 

171 

2 

Time  345 

172 

3 

sec  . 

174 

3 

176 

1 

Total  F.nvimcrntton  Problem  68 
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Undomlnated  Solutions 


Critical 

Path  Project  Length 

Cost 

Decision 

Jobs 

76 

153 

250 

®1,3 

®2,3 

"3,3 

"4,3 

^5.3 

76 

149 

300 

^1.3 

=2.3 

"3,3 

"4,3 

^5,2 

76 

147 

350 

^.3 

"3,2 

"4.3 

h,2 

64 

145 

400 

"3,3 

"4.3 

■‘’5.3 

62 

143 

450 

‘^.1 

"3.2 

"4.3 

®5,3 

60 

141 

550 

^1,1 

h,3 

"3.1 

"4,3 

^5.3 

Critical 

Path 

Project  Length 

Length 

Number 

Length 

Number 

56 

72 

141 

18 

58 

36 

142 

9 

60 

18 

143 

18 

62 

18 

144 

12 

64 

18 

145 

20 

76 

54 

146 

45 

78 

27 

147 

14 

148  36 

149  18 

150  33 

151  16 

153  4 

Time  345  sec. 


Total  Enumeration  Problem  69 


Undomlnatcd  Solutions 


critical  Path  Project  Length  Cost 


necislon 


76 

141 

76 

139 

76 

137 

56 

136 

56 

135 

56 

133 

56 

131 

Critical  Path 
Length  Number 

53 

48 

54 

24 

55 

36 

56 

30 

61 

18 

76 

81 

Time  354  sec. 

F.mimeratlon  Problem  JO 


250 

^1.3 

CM 

300 

^1,3 

^2.3 

350 

^1.3 

^2,3 

400 

S-t  ^ 

^2,3 

450 

^2,3 

500 

^1.1 

^2,3 

600 

-1,1 

^2,3 

project 

Length 

Length 

Number 

131 

1 

133 

3 

134 

3 

135 

5 

136 

8 

137 

10 

138 

11 

139 

17 

140 

13 

14' 

18 

19 

143 

15 

144 

19 

145 

18 

146 

14 

157 

15 

143 

13 

149 

8 

150 

12 

151 

5 

"3,3 

^3.2 

^3.2 

^3,2 

^3,2 

^3,2 

^3.1 


Jobs 
S 


4,3 

’4,3 

®4,2 

^4,2 

^4,3 

^4,2 

^4,2 
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Undomlnated  Solutions 

Critical  Path  Project  Length  Cost  Decision  Jobs 


78 

115 

250 

^.3 

=2,3 

^3.3 

=4.3 

^5,3 

78 

109 

300 

^,2 

=2,3 

^3,3 

=4,3 

®5,3 

78 

108 

350 

^1.2 

^2,3 

S,2 

=4,3 

^5,3 

78 

107 

450 

^,2 

^2.3 

®3,1 

=4,3 

^5,3 

78 

106 

500 

^1.2 

^2,3 

^3.1 

=4,2 

^5,3 

Critical  Path  Project  Length 


Length 

Number 

Length 

Number 

59 

108 

106 

3 

70 

36 

107 

3 

74 

18 

108 

21 

78 

81 

109 

12 

110 

11 

111 

9 

112 

25 

113 

20 

114 

48 

115 

18 

116 

22 

117 

11 

118 

14 

119 

6 

120 

7 

121 

6 

l?2 

4 

123 

2 

124 

1 

Total  Enumeration  Problem  71 


Time  304  sec. 


Undomlnated  Solutions 


Decision  Jobs 

^2,3  ^3,3  ^4 
®2,3  ^3,3  h 
^2,3  ®3,2  ^4 
^2,3  ^3,2  ^4 


Critical 

Length 

Path 

Number 

Project 

Length 

Length 

Number 

47 

108 

70 

8 

52 

36 

71 

17 

53 

18 

72 

34 

63 

81 

73 

42 

74 

37 

75 

27 

76 

25 

77 

24 

78 

18 

79 

9 

80 

2 

Time  228  sec. 

Total  Enumeration  Problem  72 


Critical  Path 

Project  Length 

Cost 

63 

77 

250 

^1.3 

63 

72 

300 

^1,2 

63 

71 

350 

63 

70 

AOO 

®1,2 

Undomlnated  Solutions 


Critical  ?ath  Project  Length 

Cost 

Decision 

Jobs 

67 

155 

250 

*1.3 

^2,3  ®3,3 

67 

149 

300 

*1.3 

^2,3  ^3,2 

“4. 

67 

147 

400 

^.3 

^2,3  ®3,1 

'4. 

Central 

Path 

Project  Length 

Length 

Number 

Length 

Number 

46 

4 

147 

2 

47 

4 

143 

2 

48 

16 

149 

13 

51 

24 

150 

13 

53 

24 

151 

29 

54 

36 

152 

34 

56 

18 

153 

34 

58 

18 

154 

38 

64 

18 

155 

25 

67 

54 

156 

18 

68 

27 

157 

17 

158 

7 

159 

7 

16C 

3 

161 

1 

Time  364  sec. 


Total  Enumeration  Problem  73 
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Undomlnated  Solutions 
Critical  Path  Project  Length  Cost 


53 

69 

250 

^1.3 

43 

68 

350 

^1.2 

41 

63 

450 

^l.l 

40 

62 

550 

^1.1 

Decision  Jobs 

^2.3  ^3,3  ^4,3 

^2.3  ^3,2  ^4,3 

^2,3  ®3.2  ^.3 

^2,3  ^3,1  ^4,3 


Critical  Path 

Project 

Length 

Leng-.h 

Number 

Length 

Number 

38 

54 

62 

3 

40 

9 

63 

3 

41 

9 

65 

3 

42 

9 

66 

9 

43 

81 

67 

12 

53 

81 

68 

6 

69 

12 

70 

15 

71 

21 

72 

15 

73 

21 

74  18 

75  21 

76  15 

77  18 

78  18 

79  15 

80  8 

81  3 

82 


Time  129  sec. 


b 
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Undonvlnated  Solutions 


Critical  Path  Pro.lect  Length 

Cos'i. 

Decision  Jobs 

67 

77 

250 

^1.3 

®2.3 

^3,3 

*4,3 

^5,3 

63 

72 

300 

^1.2 

^2.3 

®3,3 

*4.3 

^5,3 

63 

71 

350 

^..2 

^2.3 

S.2 

*4,3 

^5,3 

63 

70 

400 

^1.2 

^3,2 

*4,2 

®5,3 

critical  path 
Length  Number 

project  Length 
Length  Number 

47  108 

70 

8 

52  36 

71 

17 

53 

72 

34 

63  81 

73 

42 

,A 

37 

75 

27 

76 

25 

77 

24 

78 

18 

79 

9 

80 

2 

Time  27S  sec. 


Total  F.nutneraclcn  Problem  73^ 
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Palrwlsf  ’rchange  Results 


Problem  66 

Due  Date  70 

Decision  Jobs 


Premium  40 
Schedule  Length 


Penalty  50 

Job  Cost  Total  Cost 


^1,3  ^2,3  ^3,3  ^4.3 
^1,1  ^2,3  ^3,3  ^4,3 


s, 

79 

250 

5,3 

S,  - 

73 

400 

5,3 

700 

550 


OPTIMUM  SOLUTION  -  NOT  FOUND 
^1,1  ^2.1  ^3,1  ^  3 


f 


TIME  36,8  sec. 


5.1 


63 


850 


450 


Problem  67 

Due  Date  135 
Decision  Jobs 


1.3 

^2.3 

^3,3 

1.1 

^2,3 

^3.3 

1.1 

^2,3 

^4.1 

Premium  50 

Penalty 

50 

Schedule  Length 

Job  Cost 

Total  Cost 

S,  , 

152 

250 

1,150 

4.3 

5,3 

S  , 

145 

400 

90D 

4.3 

5,3 

S  „ 

147 

550 

650 

4.3 

5,3 

OPTIMUM  SOLUTION  -  FOUND  TIME  56  sec. 


Problem  68 


Due  Date  145 

Premium  30 

Penalty 

100 

Decision  Jobs 

Schedule  Length 

Job  Cost 

Total  Cost 

^1.3 

^2,3 

^3,3 

^4,3 

^5.3 

169 

250 

2,650 

C 

s 

s,  , 

s,  . 

Sc  •, 

154 

400 

1,300 

1,1 

^2,3 

j,3 

4,3 

5.3 

c 

s 

S, 

Sc  o 

147 

450 

650 

1.1 

‘2,3 

'3,3 

4,2 

5,3 

^1.1 

h,3 

S^  2 

1  ^ 

^4,2 

145 

500 

500 

OPTIMUM  SOLUTION  -  FOUND 

TIME  68.5 

sec . 
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Pairwise  Interchange  Results 


Problem  69 


Due  Date  142 

P7.esalum  30 

Penalty  60 

Decision  Jobs 

Schedule  Length 

Job  Cost 

Total  Cost 

^1,3  ^2,3 

^3,3 

^5.3 

153 

250 

910 

^1,1  ^2,3 

^3,3 

S.3 

145 

400 

580 

^1,1  ^2,3 

^3,2 

=4,3 

^.3 

143 

450 

510 

OPTIMUM  SOLUTION  -  FOUND  TIME  59*2  sec. 


Problem  70 

Due  Date  13^  Premium  50  Penalty  35 

Decision  Jobs  Schedule  Length  Job  Cost  Total  Cost 


1.3 

^2,3 

^3,3 

=4,3 

=3,3 

1.3 

^2,3 

^4.2 

=4,3 

^5*3 

1.3 

“2.3 

^3,2 

=4,2 

=5,3 

l.i 

^2.3 

"J.2 

=4,2 

=5..3 

OPTIMUM  SOLUTION  -  FOUND  TIME  70  sec. 


141 

250 

495 

139 

.  "^0 

475 

137 

350 

455 

133 

5or 

450 

Problem  71 


Due  Date  105 

premium  50 

Penalty 

80 

Decision  Jobs 

SchtJule  Length 

Job  Cost 

Total  Cost 

c 

r*  c 

s. 

Sr  T 

115 

250 

1,050 

"2,3  3,3 

4,3 

5,3 

c 

s 

S,  „ 

Sr  , 

109 

300 

5?0 

1*2 

2.3  ' 

4,3 

5,3 

c 

S 

S,  , 

Sr  , 

IOC 

350 

590 

1,2 

2,3  3,2 

4.3 

5,3 

OPTIMUM  SOLUTION  -  POLAND  TIME  48 

sec . 

Prol  1 

cm  72 

Due  Date  30 

Fromium  10 

Penalty 

100 

Decision  Jobs 

Schedule  Length 

Job  Cost 

Total  Cost 

c 

s 

s 

77 

250 

2,930 

"1.3 

"2,3  3,3 

4.3 

5.3 

=  2,2 

s 

4-  ,  J  J  ,  J 

^4.3 

^5.3 

72 

300 

2,500 

Q 

s  s  ^ 

S.  , 

71 

350 

2,450 

1  .2 

2»3  3,2 

4,3 

3.3 

c 

S 

S,  „ 

Sr  , 

70 

400 

2.4C0 

^1.2 

^2.3  ^3.2 

4.2 

5 , 3 

OPTL. 

i  M  :;oLuriON  -  n 

it  1  > 

riME  42 

sec . 
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Pairwise  Interchange  Results 


Problem  73 

Due  Date  160 
Decision  Jobs 

®1,3  ^2,3  ^3,3  ^4,3 


’5.3 


^1,3  ^2.3  ‘’3.2  "4.3  "5.3 

OPTIMUM  SOLUTION  -  FOUND  TIME  34.4  sec. 


Premium  10 

Schedule  Length  Job  Cost 
155  250 

149  300 


Penalty  20 

Total  Cost 
200 
190 


Problem  74 

Due  Date  66 
Decision  Jobs 


^1.3  ^2.3 

^3.3 

=5.3 

69 

OPTIinJM  SOLUTION  - 

NOT  FOUNE 

1  TIME 

18  sec. 

^1,1  ^2,3 

^3,2 

=4,3 

S.3 

63 

Problem  75 

Due  Date  70 

Prei 

Decision  Jobs 

Schedule  ' 

^1.3  ^2.3 

^3.3 

=4,3 

^5.3 

77 

^1.3  ^2,3 

^3,3 

=4.3 

^5.3 

72 

^1.2  ^2,3 

^3.2 

=4,3 

^5.3 

71 

^1,2  ^7,3 

^3.2 

=4,2 

^5,3 

70 

Premium  10 

Schedule  Length  Job  Cost 

250 


450 


250 

300 

350 

400 


Penalty  200 

Total  Cost 
850 


420 


Penalty  60 

Total  Cost 
670 
420 
410 
400 


OPTIMUM  SOLUTION  -  FOUi,:  TIME  54.5  sec. 
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Problem  76 

Due  Datd 


OPTIMUM  SOLUTION 

Problem  77 

Due  Date 


Pairwise  Interchange  Results 


premium  30 

Penalty 

70 

Schedule  Length 

Job  Cost 

Total  Cost 

482 

836 

1,676 

473 

932 

1,142 

465 

936 

836 

-  JNKNCTW.  TIME  ?*?  sec. 


Premium  30 

penalty 

50 

Schedule  Length 

Job  Cost 

Total  Cost 

294 

836 

1,536 

284 

878 

1,078 

268 

1,070 

710 

OPTIMUM  SOLUTION 


-  UNKNOWN,  TIME  441.0  sec. 
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Multiple  Pairs  Interchange  Results 


Problem  56 


Due  Date  70 

Premium  40 

Penalty  50 

Decision  Jobs 

Scnedule  Length 

Job  Cost 

Total  Cost 

^1.2  ^2.3  ^3,3  ^4,3 

^5,3 

79 

250 

700 

^1,1  h,2  ^3,2  ^4,2  ^5,3 

OPTIMUM  SOLUTION  -  NOT  FOUND  TIME  40.9  sec. 

550 

600 

^1,1  ^2,1  ^3.1  ^4.3 

^5,1 

63 

850 

450 

Problem  67 

Due  Date  135 

Premium  50 

Penalty  50 

Decision  Jobs 

Schedule  Length 

Job  Cost 

Total  Cost 

^1.3  ^2.3  ^3.3  ^4,3 

^5,3 

152 

250 

1,150 

^1,1  ^2,3  ^3,1  ^4.3  ^5,3 

OPTIMUM  SOLUTION  -  FOUND  TIME  38. 

137 

,6  sec. 

550 

650 

Problem  68 

Due  Date  145 

Premium  30 

Penalty 

100 

Decision  Jobs 

Schedule  Length 

Job  Cost 

Total  Cost 

^1,3  ^2,3  ^3.3  ^4,3 

169 

250 

2,650 

S,  ,  S,  ,  S,  ,  S.  , 

1,1  2,3  3,1  4,-: 

S.3 

143 

600 

540 

^1,1  "2,3  ^3,2  ^4,2 

^5,3 

145 

500 

500 

OPTIMUM  SOLUTION  -  FOUND  TIME  52.8  sec. 
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Multiple  Pairs  Interchange  Results 

Problem  69 

Due  Dat^  142  Premium  30  Penalty  60 

Decision  Jobs  Schedule  Length  Job  Cost  Total  Cost 


1.3 

^2,3 

^3,3 

^4.3 

^5.3 

153 

250 

910 

S.,  , 

S-  , 

S,  o 

Sr  o 

141 

650 

620 

1.1 

2,3 

3.1 

4,2 

5,2 

So  o 

So  o 

S,  o 

c 

143 

450 

510 

1  1  ^ 

2,3 

3,2 

4,3 

5,3 

OPTIMUM  SOLUTION  -  FOUND  TIME  59.7  sec. 


Problem  70 


Due  Date  134 

Premium  50 

Pent  lay  35 

Decision  Jobs 

Schedule  Length 

Job  Cost 

Total  Cost 

^1.3 

^2,3 

^3,3 

^4.3 

^5.3 

141 

250 

495 

^1.3 

^2,3 

^3,2 

^4.2 

S,3 

137 

350 

455 

^1,1 

•^2,3 

^3.2 

^4.2 

^5,3 

133 

500 

450 

OPTIMUM  SOLUTION  -  FOUND  TIME  60.4  sec. 


Problem  71 

Due  Date  105 
Decision  Jobs 


^1,3  ^2,3  ^3.3  ^4,3 
^1.2  ^2.3  S.;  "^4.3 


Premium  50  Penalty  80 

Schedule  Length  Job  Cowt  Total  Cost 

S^  2  115  250  1.050 

S,\  108  350  590 

5,3 


OPTIMITM  SOLUTION  -  FOUND  TIME  48.6  sot. 
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Multiple  Pairs  Interchange  Results 


Problem  72 


Due  Date  50  Premium  10  Penalty  100 

Decision  Jobs  Schedule  Length  Job  Cost  Total  Cost 


^1.3 

‘‘^2,3  ^3,3 

‘‘’4,3  S,3 

77 

250 

2,950 

^1.2 

^2,3  ’"3,2 

^4,2  ^5,3 

70 

400 

2,400 

OPTIMUM  SOLUTION  -  FOUND  TIME  43.3  sec. 


Problem  73 


Due  Date  160 

Premium  10 

Penalty 

20 

Decision  Jobs 

Schedule  Length 

Job  Cost 

Total  Cost 

^1.3  ^2,3 

s,-  . 

^•3 

S.3 

250 

200 

^1,3  ^2,3  ^3,2 

"4,3 

S.3 

300 

190 

OPTIMUM  SOLUTION  -  FOUND  TIME  50.4  sec. 

Problem  74 

Due  Date  66  Premium  10  Penalty  200 

Decision  Jobs  Schedule  Length  Job  Cost  Total  Cosi 

S,  _  S„  -  S.  ,  S,  ,  S,  -  6S  250  850 

1,3  2,3  3,3  h,3  5,3 

OPTIMUM  SOLUTION  -  NOT  FOUND  TIME  22.8  sec. 
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Multiple  Pairs  Interchange  Results 


Problem  75 


Due  Date  70 

Premium  30 

Penalty 

60 

Decision  Jobs 

Schedule  Length 

Job  Cost 

Total  Cost 

^1,3  ^2,3  ^3,3  ^4,3 
^1,2  ^2,3  ^3,2  ^4.2 

^5.3 

^5.3 

250 

400 

670 

400 

OPTIMUM  SOLUTION  -  FOUND  TIME  37.5  sec. 

Problem  76 

Due  Date  470 

Premium  30 

Penalty 

70 

Schedule  Length 

Job  Cost 

Total  Cost 

482 

836 

1,676 

448 

1,689 

1,029 

441 

1,739 

867 

OPTLMUM  SOLUTION  TIME  468 

sec. 

Problem  77 

Due  Date  28C 

Premium  30 

Penalty 

50 

Schedule  Length 

Job  Cost 

Total  Cost 

/ 

836 

1,536 

270 

1,681 

1,381 

276 

1,208 

1,088 

267 

1,304 

914 

TIME  380  sec. 
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APPEN’'/IX  E 
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APPENDIX  E 

Branch  and  Bound  Prograrn  Description 

A.  Main  Program,  Subroutines  and  Tapes; 

One,  Bonnie,  Clyde,  Five 

The  program  will  call  tape  .99.  to  find  the  number  of  tlie 
problem  tape.  Format  (14) 

B.  Problem  Tape; 

1.  First  Line  LIM,  TIME,  BON,  TI.EX,  LAST,  IJND,  JY 
Format  (15,  3F5.0,  315)  where 
LTM  =  Number  of  Decision  Nodes 
TIME  =  Due  Date 
BON  -  Penalty 
TLEX  =  Premium 

LAST  =  Node  Number  of  Artificial  Tlnlsh  Job. 

t 

IJND  -  Ni'ml)cr  of  "other"  Interdependency  Constraints 
JY  “  Number  of  Decision  Nodes  with  More  than  Three  Alternatlv 
"LIM  +  Sx”  Lines 

Job  Number,  Job  Number,  Job  Number,  Cost,  Cost,  Cost 
Format  (315,  3F5.0,  IS! . 

If  there  are  not  exact  1 y  three  alternatives  then  in  the 
l.'ist  position  of  Che  line  put  l!ie  difleri'ace  between  the  number 


ai\d  three. 
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3.  The  next  sot  of  lines  show  precedence  relations  In  the 
reduced  network. 

Predecessors  Numlier,  Successor  Kutnher,  Time  Format 

(2l/».  F/».0) 

4.  A  line  with  zero  In  column  I 

5.  Next  IJK’D  Lltics. 

"Other"  Interdependency  Constraints 
Key,  Job  Number,  Job  Number 
Format  315 
Key:  1  for 

2  for  ^ 

3  for  - 


ONE  MADTRN 
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OMF  fAniPM 


0f)/?2  lii2T.2 


0  n  0 1  n 

0  0  0  20 

n  If'Ff.'S  1  OM  Mt’ATiK  50  00  )  -  KL  1  ST(  200, 2  )  ,OOST(  200  ) 
n  1  MFf.'S  ION  FP.n(  non  ) ,  on  (  5, 050  )  ,TR2  ( 10  ) ,  TB3(  10  ) 

00030 

niMFf.'GinM  10(02) 

0  0  0  It  0 

COffO^’  MpATM,Rn,TBl,TB2,Tn3,TCOST,K,LIM,OOST 

OOOfiO 

PALI-  non.CK 

ccocn 

RFAD  (00,100)  KFILF 

00070 

ion 

F0P31AT  (|4) 

GOOOO 

130 

FOP” AT  (315, 3F5. 0,15) 

OOOOO 

14  0 

FOPMAT  ( 21 4, F 4.0) 

00100 

MCn!!llM  =  l 

0  n  1 1 1) 

Trnr.T=o 

001  20 

PTAP  (KFII.P,  200  )  Lll', TIME, RON, TLhX, LAST,  I.JMP 

00130 

JY  =  >IY  +  Lr' 

OOl'tO 

LF'’  =  3*JY 

0U150 

K'^nF  =  I.FI*+2 

OfMCO 

PO  in  |-l,KrCF 

00170 

101 

KI.!ST(  1  ,2)=n 

0  010  0 

KLir.T(LF‘'+l,l)=l 

00100 

KL  I  5T(  LF‘’+2,  1  )  =LAr.T 

00200 

LF‘'  =  l  Ff'"2 

00210 

290 

fop:', AT  (  !5,3F5.0,3I5) 

00220 

ri=l 

0  0  23  V> 

no  i.n  i=i,LPt‘,3 

0  0  2  4  0 

L=(+2 

0  0  2  00 

RFAP  (^'FI  I.F,130)(KLIST(  J,l),.)  =  l  ,L)  /  (CnST(  J), 

n  Ti  2  0  0 

f:=f:+i 

0  0  2  70 

IF(  lOfr'-l  )  -00  )  230,239,23^ 

0  0  2  0  0 

23S 

no  152  0=1,2 

0  0  20  0 

fo:  =  1  +'■ 

0  0  3  0  0 

IF(COCT(0.::))  149  ,  150,15  0 

00310 

150 

iF(-:OOT('‘f')-^OF^(l  )  )lfU)52,ir.2 

0  0  3  2  0 

151 

KTFf‘Pl=Pn5T(l‘0 ) 

C  G  3  3  0 

KTF*T2  =  KL  I  5T('P',  1 ) 

)  0  >  1,  0 

CPST('’4')="'0;'.T(  1  ) 

n  7  r,  r 

!<!.I5TC‘?’,1)-KL!0T(I,1) 

00500 

C05T(  1  )  =KTF;'P1 

0  0370 

KLI5T(  !  ,1  )=KTr  ’n2 

00 '00 

152 

roMT  ro'f 

('  ( •  3  0 

140 

TrOOT  =  T''r-^  f  ♦rpt;!  (  i  ) 

•MiiiPO 

SOr  ^rpOT'  1  ) 

0  0  4  1 

'>  H  0 

np  ir^  j7|,!. 

0  0  0 

00.31  (  ,1  )  =r05T  (  .!  )  -POO 

,  t'. !  3  r. 

K7'=KI  1  ^'T(.',  1) 

:  04  40 

KO 

op  AT"  =  3  ‘“'PP  -  I  -  2 

'  i !'  !j  '■  r 

1 0 

r  01;  rO 

DP,MT  p  ^Kriip" 

'  1  il  / 

2 

FCP’-'T  ' 

;  0  0 C 

•  •p*  rii(  1  1 

>.n;  T  M  f  1  ••  T  7  _  T  ,  j  /  4  -1 

7  7  r  7 

1  7'^ 

1  ^  ^ 

rr  «n  LI ,  12,  r'7P(  !  1 

r  7  7  '7 

ir(tn  1  n  ,  1  ■ 

0  <■  3  0 

7  0  0 

i  ,  1  )  r.fPAT"(  I  1  > 

'7  ■  f,  '7 

:  -t  -7"i';T!if  l  2  ) 

■  r  r  ■;  0 

1  -  1  ♦  1 

JY 


0nF70 
onFao 
oo5on 
noRon 
OCGIO 
00620 
00630 
00640 
00631, 
00  6  60 
onr.70 
oonsc 
00600 
00700 
00710 
G0720 
00730 
00740 
00750 
00760 
00770 
00730 
00700 
OCSOO 
OORIO 
00820 

00830 
00840 
00850 
00860 
00370 
00880 
G0890 
009C0 
a  00 1C 
00020 
0  003G 
00940 
00950 

nooeo 

00070 

00930 

00900 

oinoo 

01010 

01020 


115  LF.X=I‘l 

IF  (I.JNH)  188,175,188. 

188  M«0 

on  i»i,ijNn 

RF6n  (KFn.F,27n)  KFY,I1,I2 
no  TO  (310,320, 320), KFY 
310  M=M+1 

N(M,l)«MPATH(il) 

M(M,2)»MPATM(I2) 

(50  TO  320 
320  no  328  K=l,3 

IMIKF=(MPATH( ll)/3)*3+K 
IF(rn5T(IMIKF))  328,321,521 
521  r(  IMIKF-MPATIK  |1))322,328,322 
322  M=M+1 

IJ(M,1)=|MJKF 

IJ(M,")=MPATH(I2) 

328  rONTINUF 

GO  TO  (320, 329, 325), KEY 
325  13=12 
12=11 
11=13 
^FY  =  2 
GO  TO  320 

329  CONTIMUF 
IJN’0=r’ 

270  FORMAT  (315) 

175  DO  190  1=1, LEX 
KnMP=KBn(|,2) 

IF  (KLL5T(K0MP,2))  182,181,182 

181  KlIST(KOMP,2)-I 

182  D"  185  J=l .LEX 

Ih  (KRn(J+l,2)-KBD(l,2))  185,184,185 

184  KRP(l,2)=J-*-l 
GO  TO  190 

185  rOMTIf'UF 

K'^n(  1 ,2)=0 

190  CONTINUE 

CALL  JORTM(KTIMF) 

PRINT  1000,K.T|MF 
1000  FORMAT  (lOH  TIMF  USFn,|5) 

CALL  CPM  (KPn,FRn,KLIST,IJ,IJNn, jO) 
CALL  JOBTM(KTIMF) 

PRINT  1000,KT|MF 

CALL  EXIT 

END 


BONNIE  MADTRM 


f^ONMIE  MADTRN 


00010 

00020 

00030 

OOOliO 

00050 

oooec 

00070 

00080 

^0100 

OOllC 

00120 

C0150 

00160 

00170 

OOlSO 

00190 

00230 

00240 

00250 

00260 

00270 

002S0 

Q0290 

00300 

00310 

00320 

00340 

0C350 

00^60 

00370 

00380 

00390 

00400 

OO-MO 

00420 

00430 

00450 

00500 

00510 

00520 

&0530 

00540 

00550 

00560 

00570 

00580 

00590 

00600 

OOGin 

00C20 

90G30 

00640 

00650 

00660 

00680 

00690 


05/22  1417,0 


1000 

1900 


2010 

2015 


2018 


2019 
20  20 
2100 


2109 


SL'RROUTINF.  TPM  (KB^,  FBD,  Kll  ST,  IJ  J  JNP,  1 0  ) 

DIMFr'SlOf'  IO(62),START(5000),SLNGTH(60,10),SL(10) 

niMFf'SION  «PATH(509n)  ,  P.O(  3, 650 )  ,COST(  200 ;  ,TB1  ( 10 )  ,  L I  STBN(  100  ) 

DIMF.NSIOM  TB2(10),LfVF(55C),FBD(600),IJ(50,2) 

DIMFMSION  IFEAS(10),n(200),LISTEM(10),LPATH(10),TB3(10),KDC(60) 
DIflENSION  LOi;T(6C),KB0(600,2),LSTD!  J(  62,2),  KL I  ST(  200,2  ),KZERO(  650) 
C0f9i0ri  *'PATH,Br,TBl,TB2,TB3,Tr0ST,K,LIM,C0ST,TIME,50N,TLEX,LM,M,J'f 
locate ( MPL'M ,  KOUM ) « (  2+ L I M)  * ( MPUM- 1 )  +  KPUM 
RFAn  (80,10)  (KnC(l),l=l,LIM) 

FOP.MAT(30|  2) 

KLNrTH  =  5nOO/(HM+2) 
no  20  l=l,KLNGTH 
KZFRO{n  =  | 

Bn(i,i)=o 
BD(2,1)=0 
BD(3,l)r  0 
LEPP=1 
K^!TP  =  1 
nr:Nf^=i 

i.PLAr=i 

M=1 

LI''F(1)=1 
LM  =  0 
|PM  =  0 

Rr-6?.>'n=9999999, 

KPATM(l)sG 
KF.OP=3*(JY)  +  2 
n{KFOp-l)=l 
KSTOP=KFOp-2 
D(KFOP^=l 
FORMAT  (1514) 

KZTOP=0 
KZPLAf:  =  l 
GO  TC  7000 
LIi  =  LnrATF(M,0) 
iQUiT=3-inc.'.!r:) 

GO  TO  3000 

!F(LL-1)  3700,2100,2015 
DO  2020  N=1,LL 
LT  =  LI.“M 

DO  2019  L=1,LT 
LP1=LISTFM(L) 

LP2=Li5TFM(L+l) 

IF  (T.’.3(LP1)-TB3(LP2))  2019,  2019,2018 
LTEf’P  =  LISTFM(L) 

LISTFM(L)«LISTFM(L+1) 

LISTFM(1.  +  1)=LTFWP 
CONTI MUF 
CONTICIIF 
MPATM(LM+2;=1 
N=L!STFM(1) 

L  =  LPATil(N)  +  LM+2 

MPATH(L)=IN+N 

START(L)=SL(N) 

IF(LPATH(f!)-LSF)  2109,2110,2110 
L=LPATH(N)+1 


I 


0700 

00710 

00720 

00740 

C0750 

00760 

00770 

00860 

00870 

00880 

00890 

00900 

00910 

00920 


f  009?'^ 


00940 

00990 

00960 


00970 


I 


f- 


0O9S0 
00990 
01000 
01010 
01020 
01030 
01050 
OlOGO 
G107Q 
01080 
01090 
01100 
OHIO 
01130 
01140 
01150 
01160 
01180 
01190 
01200 
01210 
01220 
01230 
01^40 
01250 
012  60 
,1270 
01280 
01290 
01300 
01310 
0H20 
0i330 
”1340 
01350 
0i360 
01370 
01460 
01470 
0148U 
01400 
01500 


I 


DO  2115  l»l.,LSF 
J=LM+l+2 

MrATH(J)«LSTnN(l-lH) 

2115  STA:^T(J)=SLNGTH(I'1,N) 

2110  Bn(i,M)=TBi(:n 
Bn(2,M)-TB2(N) 

BP{3,H)»T83(N) 

IF(LL-l)  3700,3700,3000 
3000  K0UNT=LPLAC+1 
DO  3080  M=2,LL 
KMTn=KNTR+l 
KZPLAr^KZPLAOH 

IF  (KZPLAC-KLMnTH)  3050,3050,3060 
3060  PRINT  3354, LEND, KZTOP,KZPLAC 
3334  FORMAT  (18'l  LFND  KZTCP  KZPLAr,3l5) 
KZPLAr=l 

3050  MENn  =  KZFPC(KZPl.Ar) 
LMFNn=LOCATF<MFMn^C) 

MPATH  (  LMFNi.'+l )  =MPATH{LM+1  )  +  l 

MPATH(LM£Nn+2)=l 

KM  =  LiSTF^^{^') 

LIOTFM(N)=MFNn 

IF  (f.PATH(LM+l))  3051,3051,3055 

3051  MPATH(LMEN0+3:=IN+KN 
SlART(LMENn+3:=SL(KN) 

GO  TO  3086 

3055  LSP=LSF+2 
J  =  0 

DO  3085  l=3,LSP 
Kn=LMt:Nn+i 

;F(l-LPATH(Kf!}-2)  3070,3065,3070 
3065  MPATH(KR)*IN+KN 
START(KR)=GL{KM) 

GO  TO  3085 
3070  J=J+1 

MP,ATH(KB)  =  LSTDIJ(J,1) 

START (i<B)=SLNGTH(vJ,  KM) 

3035  rONTIflUE 

3036  Rn(l,MFMn)=TRl(KN) 
BP(2,MFfir)=TR2i:'N) 
Bn(3,MFNn)=TB3(KN) 

3080  CONTI MUF 

IF  (KOl'MT-LFNn)  3120,3120,3155 
3120  KTEf‘P  =  0 

no  3150  LPlJSM  =  KOU'''T,LENn 
Ml KF=LENn+LL-l-KTEMP 
MFL  =  l.FNn-KTEf'P 
Lr'E(M(KF)=L;'/F(MFL) 

3150  KTEMP=KTFMP+1 
3155  no  3200  r'  =  2,LL 
M|KF*KOUNT+N-2 
3200  LIVE(MIKF)«LISTEM(N) 

LENn=LEHn+LL-l 

3700  MPATH ( LM+ 1 ) ^MPATH ( LM+ 1 ) + 1 
IFCLL)  5000,5000,4000 
4000  IF(MPATH(LM+l)"LIM)1000,2i:5,2125 
2125  I F(Bn(3,M)-RESRnn) 2150, 2130, 5000 
2130  |RM=inf’+l 

LISTBM(  inN)-»M 
on  TO  4036 
2150  BFSRNn=Bn(3,M) 


01510 

ipn=i 

01520 

LIRTOMC  IRN)=fi 

01530 

4036 

CALL  JORTfMMTI) 

015»t0 

PRIfiT  3333,LFNn,MHNr,nn(l,M),RO(2,M),BD(3,M),MT| 

01550 

3333 

FORMAT  (22M  LFMR  MFMO  ROUMOS  TI MF, 2 1 5, 3F10 . 1 J 5 ) 

015  CO 

4050 

IF  (l.PLAC-LFMP)  4100,4060^4060 

01570 

4060 

LPLAC=1 

01580 

(?0  TO  4150 

01590 

4100 

LPLAC=LPLAr+l 

01600 

4150 

IFCLFMO-IBM)  6000,6000,4151 

OICIO 

4151 

M=LIVF(LPLAC) 

01620 

LI’  =  LOCATF{M,0) 

91630 

IF  (Bn(3,M)-BFSBNn)  4000,4000,5000 

OlGliO 

5000 

IF  (LFMO-LPLAC)  5001,5001,5005 

0-1650 

5001 

LFLAC=1 

01650 

CO  TO  5011 

01670 

5005 

DO  5010  N=LPLAC,LENn 

01680 

5010 

LIVF(N)=LIVE(N+1) 

01690 

5011 

LENn=LFND-l 

01700 

IF(KZTOP-KLNCTH)  5020,5030,5030 

01710 

5020 

KZTOP=KZTOP+l 

01720 

GO  TO  5935 

01730 

5030 

KZTOP=l 

01740 

5035 

KZERO(KZTCP)=M 

01750 

GO  TO  4150 

01760 

7000 

HIM=MPATM(LM+1) 

01770 

IM=(KnC(MIK+l)“l)*3 

01730 

Kf'=MPA.TH(LM+l)  +  2 

C1790 

MIN=KPr(f|M+l) 

OISOO 

CO  TO  1900 

01810 

8000 

LL  =  0 

01820 

KSTnp=KFOP-2 

01830 

DO  8011  l=l,KSTOP 

01340 

8011 

n(  1  )=0 

01850 

LSF  =  0 

01?60 

IF  (K!'-2)  9200,9260,9210 

01870 

9210 

no  9250  L=3,KM 

01330 

M!KE=LM+L 

OISOO 

l=MPATM(f'i:<F) 

01900 

LSF=LSF+1 

01010 

LSTDI J(LSF,1)=1 

01920 

LSTPI  jUsF,2)=KLIST(  1 ,1) 

01930 

9250 

0(  1  )-l 

ai940 

9260 

L5F=LSF+1 

01950 

LSTOI  J(LSF,;i)=KEOP 

01960 

LSTnM(L5F,2)=KLIST(KF0P,l) 

01070 

DO  9900  K=l, lOUIT 

01980 

IO=IN+K 

01990 

D( IO)=l 

02000 

DO  3220  1=1, LSF 

02010 

iF  (KLIST(IO,l)-Lr>TD|J(|,  2))  8230,8220,8220 

02020 

8220 

CONTINUE 

02030 

8230 

LPATM(K)=I 

02040 

IFdUN'O)  3030,8800,8030 

02050 

"030 

IF  (K,N-2)  8806,8800,8050 

020CO 

,050 

DO  3300  |L=1,IJND 

02070 

|f'|KF=|d(ll.,l) 

02080 

|MFL  =  |..'(  lt.,2) 

02090 

TORT  =  n(  IMIK6)  +  n(  |f»FL)-2. 

02100 

IF  (TOST^  8300,9900.8300 

190 


c 


i 
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n21l0 

3300 

rPf.'TIMLir 

^2120 

8800 

IRO-f-1 

J213U 

IF(IDO)  8805,8810,8805 

.12140 

8805 

DO  8306  IZ»1, IRO 

''•2150 

MIKE  =  l.f'+2+IZ 

:2160 

8806 

SIMCTIK  IZ,K)=START(MIKF) 

.'^2170 

8810 

CALI.  TOME(KRD,KLIST,D,FBD,LSF,LSTDIJ,  1  ,  1 0,  SLNCTH,  SI ) 

h21S0 

CKOST^onSTC 10) 

J2190 

CALL  MOMFY  (CKOST) 

’2200 

IF(TB3(K)-BFSRL'n)  9400,9400,9  500 

r.2210 

9400 

LL=LL+1 

02220 

LISTFfULL)=K 

02230 

no  TO  9900 

02240 

9500 

CONTINUE 

J2250 

9900 

D( I0)=0 

('2260 

no  TO  2010 

02270 

6000 

DO  6020  IK=1,|RN 

02280 

l=LI5TRN( IK) 

0  2290 

MIKF=LOCATE( l,0)+2 

12300 

JPO=LIM+2 

02310 

PRINT  e070,Rn(l, 1 ),Bn(2, I ),BD(3, 1 ) 

12320 

no  6025  J=3,JPa 

12350 

MIKE=MIKF+1 

0  2540 

N|KE=‘^PATU(MIKF) 

02350 

6025 

LOUT(J)=KLIST(Nli<E,l) 

02360 

6020 

PRINT  6050, 1, (LOUT(J),J»3,JPQ) 

)2370 

DO  6030  M=1,!BN 

12380 

6030 

PRINT  60Gn,BFSRNn,LISTBN(M) 

02390 

PRINT  6071, KNTR 

12400 

C071 

FORMAT(5H  KNTP, 15) 

12410 

G050 

FORMAT  (1X,I4,5M  PATH, 11 1 6/ (IIX, 10 1 6) ) 

r'2420 

GOGC 

FORMAT  {1X,F10.1,5X, (4) 

:.2430 

6070 

FORMAT  (1X,7H  BOUNDS , 3F10 . 1 ) 

'2  440 

PRINT  6999, (KOr( 1 ), 1=1, LIM) 

:2450 

6999 

FORMAT  (IIH  NODF  ORDFR,15|4) 

0  2  460 

RETURN 

'2470 

.083+2. 

783 

END 

i; 
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CI.YPE  MADTRN  /  22  I'llD.f) 


cocin 

oco2n 

noo’,0 

00040 
0005C 
00060 
00070 
0  003  0 

poooo 

00100 

00110 

00120 

00130 

00140 

0C150 

0017C 

00130 

CC190 

C0200 

0C21C 

00220 

00230 

00240 

00230 

00370 

00380 
0  0400 
00410 
00420 
0  0  li  3  0 
0  0  4  4  0 
0  04C0 
00470 
0  0  48  0 
00400 
00500 
00510 
00520 
00530 

>.  uooo+.rC'C 


SliRROMTIf'F  Tnf‘r(f'T’p  J(l.|ST,D,Fnn,L5TJ.STni.J,  I  f'SFPT,  1 0, SI  MPTH.SL) 
DP'FNSinM  SIT'PT'KGO,  10),SL(10) 

niMpMSIPN  n(2on),LSTnU'(C2,2).KLIST(  200,2),Kr’t'»(60i\2),LP/'TH{10) 
^l^'FflSIOM  TI.MPTM(  200  ),TRl(ir),FPn(GO0),MPATM(3000)  ,BP(3,f.30) 

D I  MF.HS  I  ON  TR3  ( 10  )  ,Tn2(  10  )  ,COST(  200  ) 

COMMON  MPATH,RP,TRl,TB2,TB3^TrOST^K,LIM,COST,TIME,BON,TLFX,LM,M,JY 

MIKF=3*JY+2 

PO  150  l=l,f'ilKr: 

150  Tl.NPfM(I)=0 
IOO=lNSFnT“l 
IFdRO)  120,137.120 
120  DO  125  I  =1,  IPO 
KO  =  LSTnivI(  I  ,1) 

125  TI.NnTM(KO)=Sl.r'nTH(  1  ,K) 

137  l=IBQ 

LST=^LST+1 

PO  200  .J=INSFPT,LST 
IF  (.I-INSFRT)  130,135,130 
130  1=1+1 

r!0',/=LSTn|0(l  ,1) 

RO  TO  155 
135  NO\.'=IO 

155  l(MnW  =  KLIST(MOW,2) 

ICO  ITEST  =  KBn(KNnVM) 

IF  (P(ITFST))  300,300,250 
250  F  =  TLNrT'!(  ITt‘ST)  +  FPP(f^Mnv.') 

!F  (TLNnTM(NPW)-F)  200,300,500 
2G0  TLNrT'KMO’.')=F 
300  KMnH  =  Knn(|'finV',  2 ) 

IF  (Kf'PW)  irn,  200,100 
200  COMTI'M'F 

TBl(K)=TLNrTI!(f^IKF) 

L5T  =  !.5T-1 

DO  900  l=IMSEnT,LST 
Kn  =  LSTP!.l(  1 ,1) 

900  SIMCTHC  I  ,K)=TLNr;TH(Kn) 

SL(K)=TLNnTH( ID) 

RETURN 

END 
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00010 

00050 

00030 

00040 

00050 

00060 

00070 

00080 

00090 

00100 

00130 

00140 

900+.316 


SURRnilTJNE  MONFY(rKOST)  FIVE  MADTAN 

niMENRION  MPATH(5On0),Bn(3,C50),TBl(10),TB2(10),TB3(:0),COST(200) 
COMMON  MPATH,BO,TRl,TB2,TB3^TCOST,K,LIM,COST,TIMF,EON,TLEX,LM,M,JY 


TB2(K)«BD(2,M)+CKOST 
A1-TB1(K)-TIME 
|F{A1)  9650,9670,9670 


9650  Al»TLEX*Ai 
GO  TO  9680 
9C70  A1»B0N*A1 
9680  T33(K)-TB2(K)+A1+TC0ST 
RETURN  . 

EKO 
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